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Page  2. 

X 

All  theoretical  and  experimental  data  about  interaction  of 
electron  with  atom  necessary  for  conducting  calculation  are  in  detail 
examined.  The  methods  of  the  optimum  organization  of  the  algorithms 
of  calculations  are  described.  Are  given  the  characteristics  of 
programs  for  computers,  and  also  the  results  of  calculations,  obtained 
with  their  aid,  in  the  form,  convenient  for  the  practical  use  during 
the  solution  of  all  possible  problems  of  applied  character. 
Calculations  are  done  for  8  elements  and  substances  with  the  atomic 

number  from  3.3  to  82  in  the  wide  region  of  wave  energy  (0. 4-6.0  MeV) 

✓ 

and  angles  of  incidence  of  electrons  (0-81°),  which  makes  it  possible 

v7  <"  ' 

to  carry  out  reliable  interpolation  and  to  use  the  obtained  results 
for  any  elements  and  complicated  substances.  '  ‘ 

Consecutive  and  systematic  presentation  of  material  makes  book  of 
available  for  wide  circle  of  specialists,  who  work  in  field  of 
application  of  electron  accelerators  for  practical  purposes,  engineers 
and  technologists,  connected  with  construction  of  different  isotopic 
instruments,  and  also  for  students  of  old  policies  of  VUZ  [Institute 
of  Higher  Education]. 
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INTRODUCTION. 

Transport  problem  of  charged/loaded  particles  (electrons  and 
positrons)  in  substance  for  period  of  several  decades  draws  interests 
of  researchers.  So  prolonged  an  interest  in  this  problem  is  explained 
by  a  large  quantity  of  problems,  which  arose  recently  :n  the  new 
branches  of  science  (space  technology,  radiation  chemistry  technology, 
nuclear  engineering,  etc.).  Experimental  research  in  these  regions  is 
conducted  by  broad  front  and  contains  the  enormous  range  of  energies 
of  particles,  beginning  from  several  keV  and  ending  with  tens  of  GeV. 
For  obtaining  the  necessary  information  usually  it  is  necessary  to 
carry  out  several  separate  experiments,  to  use  precision  instruments, 
which  is  natural,  it  requires  large  efforts/forces  and  significant 
financial  expenditures. 


Much  was  done  recently  in  direction  of  theoretical  description  of 
process  of  transfer  of  particles  through  substance.  Now  the  theory, 
based  on  the  solution  of  the  kinetic  equation  of  transfer,  although  it 
is  sufficiently  perfect  from  a  mathematical  point  of  view,  it  can, 
strictly  speaking,  to  be  used  only  for  the  infinite  homogeneous  media. 
Account  in  this  theory  even  of  the  simplest  boundaries  is  complex 
problem.  Theory  unable  to  correctly  consider  fluctuations  in  the 
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Large  difficulties  appear  during  attempt  at  consecutive  account 
in  equation  of  transfer  of  secondary  radiation:  6-electrons, 
bremsstrahlung  or  products  of  nuclear  interactions. 

All  these  difficulties  led  to  creation  of  whole  series  of 
approximation  methods  of  calculation.  However,  the  essential 
simplifications,  utilized  in  them,  lead  to  the  noticeable  divergence 
of  calculations  from  the  results  of  experiment,  which  does  not  always 
succeed  in  eliminating  by  the  selection  of  semi-empirical  constants  or 
relationships/ratios.  But  the  difficulties,  connected  with  the 
consideration  of  complex  boundaries,  still  remain  in  this  case. 

Page  4. 

Thus,  was  explained  that  there  is  class  of  problems,  such,  as 
calculation  of  passage  of  particles  through  protection  of  complex 
form,  contribution  of  secondary  radiation  after  protection,  energy  and 
angular  particle  distributions,  backscattering ,  dependence  of 
probability  of  forming  interstitial-vacancy  pairs  as  functions  of 
depth  of  penetration  of  particles,  obtaining  spatial  distribution  of 
absorbed  in  substance  energy,  which  cannot  be  correctly  solved  by 
existing  theory  or  its  approximate  modif icat ions ,  but  experimental 
methods  they  are  excessively  labor-consuming.  The  alternative  method 
of  the  solution  of  problems  by  the  radiation  transfer  is  the  method  of 
Monte  Carlo  (method  for  random  testing).  This  method  has  already  long 
ago  been  successfully  used  for  the  calculations  of  the  passage  of 
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neutrons  and  7-quanta  through  the  substance. 


Application  of  method  to  calculation  of  passage  of  charged/loaded 
particles  and,  in  particular,  electrons  was  not  so/such  intense.  Is 
explained  this  fact  by  the  large  interaction  cross  section  of  the 
charged/loaded  particles  with  the  substance,  as  a  result  of  which  for 
obtaining  the  data  even  with  10%  accuracy  necessary  to  fulfill  lO’-lO* 
random  tests,  which  virtually  cannot  be  realized  on  the  contemporary 
computers . 


However ^  when  into  diagram  of  Monte  Carlo  method  was  introduced 
series/row  of  improvements,  time  for  conducting  calculation  noticeably 
decreased,  and  advantages  of  method  -  sufficiently  common  algorithm  of 
calculation,  automatic  account  of  boundary  conditions,  even  most 
complicated,  largest  volume  of  information,  obtained  in  one  cycle  of 
calculations,  and  others  -  select  it  in  comparison  with  others  by 
Monte  Carlo  method,  it  can  be  done  sufficiently  high  and  in  many 
instances,  compared  with  attained  in  experiment,  if  necessary  accuracy 
in  initial  data,  that  characterize  elementary  report/event  of 
interaction  is  provided.  For  the  electrons  with  the  energy  to  10  MeV 
this  condition,  as  a  rule  is  satisfied. 


Examination  of  different  aspects  of  application  of  Monte  Carlo 
method  for  solving  of  those  formulated  above  problems  composes  content 
of  this  work.  We  did  not  approach  the  consecutive  presentation  of  the 
essence  of  the  Monte  Carlo  method.  At  present  there  are  several 
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I  excellent  monographs,  dedicated  to  this  question. 

i  Page  5. 

1 


The  target,  by  which  we  were  guided,  first,  to  let  to  reader 
fundamental  information  for  conducting  such  calculations,  including 
initial  data,  tested  by  different  authors  versions  of  networks, 
recommendation  regarding  the  organization  of  programs  and  different 
methods  of  accelerating  the  calculation.  Secondly,  so  to  systematize 
the  results  of  our  calculations  so  that  them  it  would  be  possible  to 
use  in  practice  for  all  possible  evaluations/estimates  during  the 
engineering  developments,  connected  with  the  application  of  electron 
beams,  isotopic  sources,  etc.  Probably,  the  results  of  calculations 
will  be  useful  for  researchers,  who  work  in  different  regions  of 
physics. 

All  calculations  are  carried  out  on  computers  BESM-3,  BESM-4  for 
series  of  elements  of  periodic  system  in  interval  of  energy  of  0.4-6 
MeV,  which  is  of  greatest  interest  for  practice.  The  mean  time, 
necessary  for  conducting  one  cycle  of  calculations,  is  40-60  min. 

Even  with  such  relatively  larger  expenditures  of  time  (during  the  use 
of  the  more  advanced  machines,  for  example  BESM-  6,  this  time  it  can 
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be  reduced  almost  by  an  order)  economic  effectiveness  of  calculations 
in  comparison  with  the  realization  of  the  corresponding  experiments  is 
indisputable . 


Rapid  increase  in  number  of  works,  carried  out  by  Monte  Carlo 
method  in  USSR  and  abroad,  is  best  confirmation  of  this  position. 


Authors  express  appreciation  to  academician  of  the  AN  [Academy  of 
Sciences]  of  KazSSSR  Zh  S.  Takibayev  for  the  constant  attention  with 
the  fulfillment  of  this  work. 
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Page  6. 

The  principal  notations,  utilized  in  the  book,  and  the  value  of 
certain  constants. 

/3-v/ c  -  velocity  of  electron. 

P=(m,/3c)/VfT-/32  -  relativistic  momentum  of  an  electron. 
p-p/v'T-^2  -  relativistic  momentum  of  an  electron  in  ones  m„c. 

E  -  kinetic  energy  of  electron, 
e  -  kinetic  energy  of  electron  in  ones  m0c2. 
e+1  -  electronic  energy  in  ones  m0c2. 
k  -  photon  energy  in  ones  m0c2. 

Z  -  atomic  number  of  element. 

n,  -  number  of  atoms  in  1  cm2  of  substance. 

a«Z/137. 

T(x)  -  gamma  function. 

N*6.02*1023  mole'1  -  Avogadro  number, 
m, *9. 1*10" 2  *  g  -  mass  of  electron. 
e=4.8*10'l#  cgs  esu  -  electron  charge. 
c=3«101#  cm/s  -  speed  of  light. 

r,  -  e2/(m,c2 )*2. 82*10* 1 3  cm  -  a  classical  radius  of  electron. 
h*=6. 625*10" 2  7  erg*s  -  Planck's  constant. 

Try 

>f,c2*0.511  MeV  -  rest  energy  of  electron. 

0*=Z2r,  2/137=5.78*lO'2,Z2  cm2. 

a,  *h2/(4ir2m,e2  )*0. 529*10'  *  cm  -  radius  of  the  first  Bohr  orbit. 


DOC  =  86063201 


Page  7 


PAGE 


5 

v 

! 


Chapter  I. 
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CHARACTERISTICS  OF  ELEMENTARY  EVENT. 
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In  this  chapter  we  will  consider  characteristics  elementary  event 
of  interaction  of  electron  with  atom,  which  lie  at  base  of 
calculations  by  Monte  Carlo.  Are  investigated  the  field  of 
applicability,  agreement  with  the  experiment  of  different  formulas, 
which  describe  elastic  electron-nuclear,  electron-electron 
interaction,  and  also  inelastic  process  of  braking  electron  emission. 


S  1.  Elastic  electron-nuclear  scattering. 


Differential  cross  section  of  scattering  in  solid  angle  of  dSi*=2ir 
sin  0d0  electron  in  nuclear  field  with  pure/clean  Coulomb  potential 
Ze/r  is  described  by  Rutherford's  formula 

d  3  _  rJZ*  1 

Ts  “  (I=sni» ’  (1, 1} 


where  p  -  relativistic  momentum  of  an  electron  in  ones  m,c.  Formula 
(1.1)  corresponds  to  Born  approximation,  and  in  it  the  effects,  caused 
by  electron  spin,  are  not  considered.  In  this  same 
approximation/approach  the  difference  in  scattering  of  electron  and 
positron  is  not  exhibited. 


More  complete  examination  of  electron  scattering  was  carried  out 


m 


"T.  «r.  -r* 


WWW 


3 


,  '  V-7  V 

«C  >  V*  •  a  T  V-  T  •  1  •  1  i 
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by  Mott  [1].  Using  Dirac's  equation,  he  succeeded  in  obtaining  exact 
expression  for  the  section  in  the  form  of  the  weakly  converging  series 


•  <L2) 


where  F-Fo+Fi  and  G-Go+Gi 


i  r / i  —  4— ) 

Fi  -  -f  i  2  +  (/  +  1)  Di ! ,]  ( -  iy  P,  (cos  0), 


G0  -=■  —  i 


T[ct«24]  F»' 


Gl  "T*  2t/2£>/-(H-l)2DHl]  (-  1)'  P,  (cos  0), 


exp[— ir.  (j‘-  i*)‘ 

* \r 

(1*- 4 

i* 

7l 

-«■)'•  i  ~ ! 

■'  j 

e,  p,  /3  -  the  energy,  pulse/momentum  and  the  velocity  of  electron 
before  the  scattering. 


Page  8. 


Thus,  expression,  which  stands  in  curly  braces  of  formula  (1.2), 
describes  relation  of  sections  of  Mott  and  Rutherford,  i.e., 

R  -  “■  | •  The  detailed  calculations  of  this  relation  as  the 
functions  of  the  scattering  angle  by  0,  atomic  number  Z  and  energy  of 
electron  E  were  carried  out  on  computer  by  Doggett  and  Spencer  [2], 
and  their  results  for  the  electrons  were  given  in  Table  1. 
Applications/appendices.  Analogous  calculation  was  carried  out  by 
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Sherman  [3]. 


Special  feature  of  formula  (1.2)  is  such,  that  scattering  cross 
section  to  large  angles,  close  to  180°,  requires  addition  of  very 
large  number  of  terms  of  series/row  and,  consequently,  data  of  table 
in  this  region  of  angles  can  be  imprecise. 

Bartlett  and  Watson  [4]  obtained  from  (1.2) 
approx imation/approach,  valid  for  small  scattering  angles 

A  -  -  1 H — cos  v  (1  —  cos  0y . ,  (1.  3) 

/  du  y  2 

where 


cos  7  -  Re 


r  +  ia,tj  r  (i-iafi) 


(1.  4) 


Page  9. 


r 

The  numerical  values  cos  X  are  given  in  Table  2  of 
applications/appendices . 


Since  formula  (1.5)  is  obtained  during  solution  of  a  precise  section 
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in  the  series/row  according  to  degrees  a  (terms,  which  contain 
a\  n  ^ 3  they  are  throvn/re jected) ,  it  is  applicable  only  for  a<<l, 
i.e.,  for  small  Z.  For  the  rough  estimates  the  formula  can  be  used  up 
to  Z=30.  As  showed  checkout  calculations,  values  R  for  carbon  in  this 
formula  coincide  with  the  accurate  results  of  Doggett  and  Spencer. 

Given  above  formulas  do  not  consider-*Paa-  shadowing  of  coulomb 
field  of  kernel  by  atomic  electrons. 


First  calculation  of  shadowing  was  carried  out  in  Born 
approximation  of  Bethe.  The  tabulated  values  of  the  shielded  section 
of.  Rutherford  are  given  in  work  [6].  According  to  Goudsmit  and 
Saunderson  [7],  the  modified  section  of  Rutherford  taking  into  account 
the  parameter  of  shadowing  rj  takes  the  form 


where 


e^=0! J/4=(hZ>^  )/(8xaaP) ,  61  -  the  so-called  angle  of  shadowing. 


As  can  be  seen  from  formula  (1.6),  effect  of  shadowing  to  become 
essential  only  in  very  small  angles,  moreover  scattering  cross  section 
is  reduced. 


Moliere  [8],  using  WKB  [Wentzel-Kramers-Brouillin]  approach, 


obtained  expression  for  calculating  parameter  of  shadowing: 
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Nigam  et  al  [9]  made  more  precise  Moliere's  theory  for  multiple 
scattering  and  for  the  parameter  rj  they  give  the  following  expression 


ij*  -  1.7  -10-»  ~jr  Z'<‘  jl  +  4a Xo  In  /. 

+  ^P  +  1.448?]j. 


(1.8) 


where  Xo=0.00825(t/l-/3 3//3)Z3  3 . 


Page  10. 

As  can  be  seen  from  Jable  1,  where  are  given  values  Y.*»  and  v..v  for 
different  energies,  the  first  are  systematically  more  than  the 
seconds,  although  only  for  large  Z  and  low  energies  of  electrons 
difference  between  them  is  substantial.  Convenient  in  the 
calculations  by  Monte  Carlo  formula  for  the  section  proposed  Spencer 
[10] .  He  spread  formula  (1.3),  which  coincides  with  a  precise  formula 
of  Mott  (1.2)  for  the  small  scattering  angles,  to  all  angles,  after 
introducing  the  correction  term  H(0)  - 


d°M  _  q 

&  ~  da  l1  +  pi  ^ 


a?  cos  7  (1  —  cos  &)  '•  -f  H  ('i>  ,  (1.  9) 


whence 


H  (0)  =•  R  —■  1 - ~  a'}  cos  7  (1  —  cos  !>)  •  •  (1.  10) 

v  a 


Relation  is  assumed  to  be  known  from  precise  calculations, 

for  example,  of  Doggett  and  Spencer.  The  account  of  shadowing  in 
(1.9)  was  carried  out  by  Spencer,  using  the  parameter  of  Moliere’s 
shadowing  Ti*.  Furthermore,  it  took  into  account  correction  for 
inelastic  electron-electron  scattering;  for  this  Z,  is  substituted  by 


ijjp 

fir 


9 

& 

I® 

$ 
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Z(Z+1).  Then  scattering  cross  section  is  registered  in  the  form 


2?L'L  UJ* — X)  f  1  4-  -4=  a}  cos  7(1  —  cos  rJ  r 

p*?1  (t  cottU-r 2i))*  [  ‘  r  2 


I-21J)1'*  1  H(0)], 


(1.  ID 


where 


—  C08  0  4- 2-»i)' •  , 


(1.  12) 


[(Z  4-l)ln4vjM]-^u  -lnj0t16Z-*^l-l-  3,33 -£)j 


-  the  so-called  factor  of  Fano, 


u  -  constant,  approximately  equal  to  5. 
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48 

10 

20 

52 

03 


257.44 


5,83 


23 

08 
04 
90 
430,90 


462,74  |  011.81  |  785,64  |  1068,79  |  1477.90 
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Convenience  in  expression  (1.12)  consists  in  the  fact  that 
correction  term  H (0)  can  be  with  high  accuracy  approximated  with  the 
help  of  the  same  polynomials  1  -  cos  6  +  2r\ ,  which  enter  into  formula 
(1.11)  s 

*  j 

H(0)=  Vi/  (0)  (l  -coe'H- 2^'.  (1.  13) 


Already  for  n=5  sum  (1.13)  differs  from  the  precise  value  of  H(0)  not 
more  than  on  !$ 


It  must  be  noted  that  account  of  shadowing,  as  it  is  conducted  by 


Spencer,  is  very  approximate.  Actually,  the  sections,  calculated  by 
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formulas  (1.11)  and  (1.12),  in  practice  do  not  differ  from  the 
appropriate  unshielded  sections  of  Doggett,  as  is  evident  from  Table 


Formulas  given  above  can  be  used  also  for  positrons.  For  this  it 


is  nect-isary  in  them  to  everywhere  replace  Z  by  -  Z.  As  far  as  cos  y 


(1.4)  is  concerned,  the  structure  of  expression  is  such,  that  it  does 


not  depend  on  sign  Z. 


Consecutive  examination  of  shadowing  is  carried  out  by  Lin  [11]. 


The  starting  point  of  this  work  is  Mott's  formula  (1A-I09v)  [18]  for 


an  arbitrary  form  of  shadowing.  Some  results  of  calculating  Lin  are 


given  in  Table  3  and  show  that  the  shadowing  leads  to  a  decrease  of 


section  during  the  small-angle  scattering  and  an  increase  in  the 


section  during  the  scattering  to  large  angles.  Recently  these  results 


were  confirmed  by  the  calculation  of  Buhring  [13]  for  mercury.  That 


that  the  shielded  section  is  more  than  the  unshielded  for  large  angles' 


of  scattering,  especially  important  during  calculations  of 


back-scatter  factor. 
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Table  3. 


a  - 

Relation  for  Au  and  Cu. 

a  .ir^np 


400 

KJB  w" 

200 

iou 

50 

Kta  V 

\ 

\ 

2 --79 

2=29 

2=79 

2=29 

l! 

2-29 

2=79 

2  =  29 

10 

0,771 

0.937 

0,630 

0.939 

0,469 

0.869 

0,307 

0,732 

30 

0.935 

1,03 

0,892 

1.00 

0,794 

0.986 

0,631 

0.966 

60 

0,996 

0.981 

0.987 

1,00 

0.967 

1.01 

0,870 

1.02 

90 

0,997 

0.984 

1,01 

0,996 

1,02 

1.02 

1,00 

1.04 

120 

1.01 

0.996 

1,03 

1.01 

1.07 

1.04 

1,11 

1.05 

150 

1.08 

1,03 

1,05 

1.01 

1.11 

1.01 

1.19 

1.06 

Key:  ( 1 ) .  keV. 
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One  should  note  that  according  to  other  theoretical 
evaluations/estimates  (Olsen  [18])  /.  l  for  0-180°.  It  is  at 

**  HHSkp 

present  difficult  to  comprehend  these  contradictions  to  theory. 


Unfortunately,  are  absent  systematic  experimental  data  over 
scattering  cross  sections,  especially  to  large  angles.  Most  complete 
measurements  ^»««« *■  in  the  interval  of  angles  of  30-150°  were 

carried  out  by  Rainwater  and  Rester  [14].  Figure  1  shows  the  part  of 
their  results  for  Al ,  from  which  it  is  evident  that  with  an  accuracy 
to  10%  experimental  sections  will  be  coordinated  with  the  unshielded 
section  of  Mott.  Remains  incomprehensible  the  systematic  divergence 
,  from  1  for  the  energies  E>2  MeV.  More  precise  experiments  of 
Dougal  [15]  and  Kessler  and  Weichert  [16]  better  will  be  coordinated 
with  experimental  data  [16]  for  Hg  with  the  calculations  of  Lin  and 
Sherman  (unshielded  section). 

Thus,  scarce  experiments  confirm  need  for  correct  account  of 


shadowing 
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Fig.  1.  Fig.  2. 

Fig.  1.  Comparison  of  experimental  elastic-scattering  cross  sections 
(point  on  graphs)  with  section  of  Mott  for  different  angles  and 
energies  of  electrons. 

Key:  (1).  MeV. 

Fig.  2.  Comparison  of  values  of  experimental  sections  with 
calculations  of  Lin  (2)  and  Sherman  (1). 
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However,  the  absence  of  systematic  data  over  the  shielded  sections 
makes  it  necessary  to  use  in  the  Monte-Carlo  calculations  the  actually 
unshielded  sections.  Error  introduced  in  this  case  will  be  noticeable 
at  the  large  scattering  angles  and  with  low  energies  of  electrons. 

There  is  one  additional  source  of  error  in  scattering  cross 
section  on  large  angle  -  radiation  effects.  As  Schwinger  [17]  showed, 
the  electron,  which  is  diffused  in  the  Coulomb  field  of  the  nucleus, 
with  probability  1  emits  photons.  Therefore  in  the  energy 
distribution  of  electrons  besides  the  quasi-elastic  peak  is  present  a 
low-energy  "tail",  caused  by  the  emission  of  rigid  photons.  The 
theory  of  this  phenomenon  is  complicated,  and  its  results  at  present 
it  is  possible  to  examine  only  as  qualitative  [18]. 

Let  us  define  radiation  correction  *r  as  relation 

/r  \l  dijd*  dE-j(  %  ju  •  (l.  14) 

where  d-£-dl,  -  double  differential  section; 

Et  and  E 2  -  energy  of  electron  before  and  after  scattering; 

(da/dK)  o  _  section  of  purely  elastic  part  of  scattering  for  example, 
sect  ior.  of  Mott )  . 

According  to  Schwinger,  for  A1  and  E.=l  MeV  v  '  0,97,  and  this 
value  was  recently  confirmed  in  experimental  work  of  Berger,  et  al 
[19].  In  earlier  work  of  Olsen  (see  [18]),  where  a  radiation 
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correction  is  obtained  in  first  Born  approximation,  it  is  shown  that 
with  energies  of  electrons  E>1  MeV  and  at  scattering  angles  ~i80°/r  l. 
Figure  3  gives  graphs  *r  (E)  for  the  different  8. 

If,  actually,  calculations  of  Olsen  are  correct,  then  account  of 


radiation  correction  can  noticeably  increase  scattering  cross  section 
by  large  angles. 
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Fig.  3.  Graph/diagrams  of  dependence  of  radiation  correction  >-.r  on 
energy  of  electron  at  different  scattering  angles. 

Page  15.  It  should  be  noted  that  as  results  of  Olsen  due  to  roughness 
of  approximation/approach  is  not  obtained,  probably,  also  essential 
dependence  *r  (Z).  However,  in  the  complete  scattering  cross  section 
radiation  correction  virtually  is  not  manifested. 

Feeding/conducting  the  result  to  study  of  problem  about 
scattering  cross  sections,  one  should  stress  that  results  of 
calculations  by  Monte  Carlo  contain  certain  systematic  error,  caused 
by  inaccuracy  of  initial  theoretical  sections. 

Further  research  both  theoretical  and  experimental  is  necessary, 
in  order  to  decrease  effect  of  this  error. 

S  2.  Electron  scattering  on  the  electron. 


Let  us  consider  kinematics  of  process  elastic-g  of  electron 
scattering  in  free  scattering  center.  Let  mx  -  mass  of  incident 
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particle  (electron),  and  P  and  Pt  -  its  pulses/momenta  before  and 
after  collision.  Let  the  mass  of  the  scattering  center  m2  -  its 
pulse/momentum  after  scattering  .  To  the  particle  scattering  m2 
rested.  Then  from  the  laws  of  conservation  of  energy  and  momentum  we 
obtain  for  the  energy,  transmitted  to  recoil  particle  E2  (kinetic 
energy  of  this  particle): 


Ei  —  2 m2c? 


pt  ca  co,-  ll 


and  maximum  energy  of  recoil  particle 


Eimt*  '  2mzC'  mj>c4+mi»c<+2 

If  scattering  center  -  atom  with  atomic  weight  of  A,  then  after 
small  simplifications,  connected  with  neglect  of  terms,  which  contain 

relation  *:<£■!  and  j ?,+!»,<?*  ^  i.we  will  obtain 

m»  "  ml 

c  2f(£+2m,c»)  _  2mE(E (-1,022)  JJ  ,,  ,,, 

" ~  A  ^  u'17) 

where  E  -  energy  of  electrons  in  MeV. 


Energy,  transmitted  to  recoil  atom,  is  small  in  comparison  with 
energy,  transmitted  to  electron,  and  by  it  are  usually  overlooked  in 
calculations  according  to  braking  of  electrons.  However,  in  problems 
associated  with  the  radiation  defects,  which  appear  under  the  action 
of  fast  electrons,  the  account  of  energy  of  efficiency  is  essential. 

Page  16. 

The  scattering  angles  of  the  particles  participating  in  the  collision 
in  the  laboratory  system  of  coordinates  is  calculated  with  the  help  of 


If  two  identical  particles  (two  electrons  in  our  case)  collide, 
then  after  collision  it  cannot  be  distinguished,  which  of  electrons 
scattered.  It  is  conditionally  accepted  that  the  electron  with  the 
larger  energy  is  scattered,  the  maximum  energy,  devoted  by  electron 
to  another  electron,  will  be  in  this  case  of  E/2,  but  not  E,  as  this 
follows  from  (1.16)  for  colliding  not  the  identical  particles  of  the 
identical  mass,  for  example  electron  and  the  positron. 


Quantum-mechanical  expression  for  scattering  cross  section  of 

electron  on  electron  is  given  by  Moller  [20]  in  the  form 

«'a  -orr  v  <1  hM  l _ 1 _*!1  ,  1  , 

ilw  "  0  «*(»  |  2)  [w1  w  ( 1  —  ui)  (Ml)*'  (1  u/)‘  1 

(T+lj*  •  (*•  ^) 

where  we  -  transmitted  energy  in  units  m,c*  and  0£w£l/2. 


From  formulas  (1.18)  and  (1.19)  it  follows  that  scattering  angl 
of  electron  (i.e.  electron,  which  has  high  energy  e(i-w))  is  equal 


u 


i  =  arc  cos 


f< l  w)  (M-2)pi 
[  <1  m)  t  i-2  »- 


(1.21) 


but  an  electron  with  smaller  energy  we  moves  at  an  angle  0,  to 


N 

A 


8 


i 


n 
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direction  of  movement  of  falling/incident  electron 


-  arc  cos 


u>(- 
5  un  ,  2 


(1.22) 


In  number  of  cases  it  is  convenient  to  use  section  of  Moller  for 


scattering  in  solid  angle  of  Q-fi+dfi.  Using  kinematic 


relationship/ratio,  is  easy  to  obtain  expression  for  do/dft  in  the  form 


,  2  r  li-  i  u* |_ _ 1 _ _  _ 2  (2M  1) _ 

‘  u"  ^  l(*  -r2)i*8in4#1  (^-+-1)*  h2)i*in*  20, 


(1.  23) 
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Maximal  scattering  angle  e  is  determined  from  (1.21)  with  w=l/2  and 


is  equal  to 


0l(n«*  —  arcr08 


(1.  24) 


Let  us  note  that  formula  (1.20)  is  valid  only  for  electron  scattering 


on  the  free  quiescent  electron  1 . 


FOOTNOTE  1 .  Moller 's  formula  was  generalized  in  work  [21]  in  the 


case,  when  both  electrons  move.  ENDFOOTNCTE. 


The  presence  in  bonding  electron  is  usually  considered  by  the 


introduction  tc  certain  threshold  they  are  of  a  size  energy  w’e,  which 


corresponds  tc  smallest  possible  transmission.  Ir.  principle  of  w'e  it 


is  possible  to  identify  with  the  electron-  binding  energy  in  the  atom. 
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although  for  these  energies  formula  (1.20)  is  yet  not  applicable. 


Total  cross  section  of  electron-electron  scattering  is  easy  to 
find,  integrating  (1.20)  in  limits  of  w  '  and  u>  — 


3  (w  >  J  dw  dw  “  z*ro~  ^  [p-  ~  T^7  -i~ 

w' 

■  [r  i i  a  3(t)  =  0-  (1,25j 


Knowledge  of  total  cross  section  e  -  e  -  scattering  makes  it 
possible  to  calculate  mean  free  path  X  with  respect  to  this  scattering 


/.  J- -  1 


The  process  of  scattering  positrons  on  electrons  was  examined  by 
Baba  [22],  and  for  section  of  this  scattering  expression 

Jj  ,  9-tr  r .  - —  itiiH _  if  I  _  ,2  ,,,2  I  _ 


j  2*r°‘  ••(€  1-2)0)'  l|/  (•+1)*“'  +  2  (t+ll  W  j 

-r  irr, )  - [t  +  t+i  -  5TTTT?  -  (.Vi  “»]|  •  “• 261 


where  OSw^l,  is  obtained. 
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S  3.  Braking  electron  emission. 


*.y 


5** 
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Inelastic  scattering  of  electrons  with  emission  of  bremsstrahlung 
intensely  was  studied  in  recent  decades.  In  one  of  the  most  complete 
surveys/coverage,  dedicated  to  this  problem,  Koch  and  Motz  [12]  give 
more  than  30  formulas  for  the  differential  cross  sections  of 
bremsstrahlung,  obtained  in  different  kind  approximations/approaches. 
The  reason  for  so  large  a  number  of  formulas  limited  on  the 
applicability  is  explained  by  the  impossibility  to  at  present  compute 
the  matrix  elements  of  the  wave  functions  of  electrons  in  shielded 
coulomb  field  of  atom. 

Are  distinguished  two  forms  of  approximate  sections  of 
bremsstrahlung:  nonrelativistic  and  relativistic.  The  first  are 
obtained  by  Sommerfield  [23]  and  are  valid  for  the  electrons  from  the 
energies  into  several  keV  and  for  the  substance  with  small  Z 
(shadowing  in  these  sections  not  taken  into  consideration).  From 
relativistic  region  is  allotted  usually  separately  a  region  of 
super-relativistic  (E>50  MeV).  The  sections  of  bremsstrahlung  for  if 
are  given  in  work  [24]. 

Of  greatest  interest  for  the  practice  are  intermediate  energies 
from  tens  of  keV  to  several  ten  MeV.  The  sections  of  bremsstrahlung 
in  this  case  are  obtained  in  the  first  Born  approximation,  which  sets 
the  specific  limitations  to  the  applicability  of  these  formulas.  They 
are  valid  for  1)  small  Z,  2)  far  from  the  maximum  energy  of  braking 
quantum,  3)  speeds  j3~l. 

However,  in  spite  of  these  limitations,  formulas  in  first  Born 
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approximation  are  extensively  used  even  there,  where  conditions 
indicated  are  not  satisfied.  This  is  explained  by  the  relative 
simplicity  of  the  expressions  for  the  section,  which  make  it  possible 
to  carry  out  the  qualitative  analysis  of  the  obtained  experimental 
results. 


Let  us  consider  only  several  formulas  most  frequently  used  in 
calculations,  Monte  Carlo  (detailed  data  is  presented  in 
survey/coverage  [12]). 


Page  19. 

Initial  in  many  calculations  is  the  expression  for  the  differential 
cross  section  of  the  bremsstrahlung  of  Bethe-Heitler  [25]: 


&(&)’  [l-'fc 

[4  (•.+  «*-«*]  +  [-><•+  I)*-!*]  - 


('#  +  !— Pu  cob  #„)* 
?pp„bii\  ^ain  cob  <J>  [4(  1 4  1)  1)— 

(*4  1  —  j>cob0)  (<„+ 1— ^cobOj) 


4~ 


Mtn*11..—  'ipi' ,,»in  'Inin  co»<P)l 

(Hi  ,pcon0)  (iu-|  1-  p„vo»0u)  I’  ‘ 

where  q^Po-p1-  k,  the  pulse/momentum,  transmitted  to  the  nucleus,  and 
qi*.  p-  p  •  -r  -  A-~ 2p0  A  cop  Vi 2pk  cca&  —  2pp0(cce  &ccb0o4  sin  0  X  «u»  con<P); 

k  and  It  -  energy  and  impulse  of  photon  in  units  =u;J  m0C; 
e „ ,  e  -  initial  and  final  energies  of  electron; 

$0 ,  p-  initial  and  final  momentum  of  an  electron; 

^  & 

dtr  e  -  angles  between  p„,  p  sndAto  respectively; 
azimuth  angle 


d&*  =  sin 6od0od<t>,  dii,  «=  sin  0d0d<P; 


I 

<«s«« 

•i'* 

•»w 

<*,v 

'»*,**• 

¥ 
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F(q,  Z)  -  the  atomic  shape  factor,  whose  numerical  values  are  given  in 
work  [83. 


Conditions  formulated  above  for  applicability  of  Born 
approximation  are  written/recorded  with  the  help  of  following 
inequalities 

(«„+!)  _  j  and  )= _  t  L 

137  K  (S+2)«o  137  v  («+2)t 

Integration  of  formula  (1.27)  can  be  carried  out  only 
numerically.  However,  for  case  of  F(q,  Z)=0  integration  can  be 
carried  out  analytically,  and  for  the  differential  cross  section  of 
photon  emission  with  the  energy  in  interval  of  k-k+dk  is  obtained  the 


express  ion 


MR 


k»U*.+UF  («  +  l 1  (  (e,ri  X)  (t  +  1)+.^  Xo  - 

- ~  'tPi>P\  A. 

U„  UK6  t-lji-/>a  ,  I  1)  (-  null  (1.  2») 

■*'  p*  *•  "•  p*  Pu"  /  ]  I 


where  Z.  -  21n(e,4'1)  LllLLE^-  xo-anK-J-l+p*);  X-21n(;~i  +  p). 


S* 
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Calculations  show  that  effect  of  shielding  is  essential  only  for 
low  photon  energy.  Usually  to  evaluate  the  possibility  of  neglect  of 
shadowing,  investigate  the  dimensionless  parameter  7,  which  has  the 
significance  of  the  ratio  of  the  atomic  radius  in  Thomas-Fermi  model 


to  the  maximum,  impact  parameter 


_ lOOfcZ  ~  ■ 

(«o  t-1)  U  t  l  -  *) 


(1.  29) 


«A 
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From  (1.29)  it  is  evident  that  when  k-*Of  7-*0and  shadowing  is 
substantial,  and  when  fc— y>i  and  screening  effect  is  negligible. 

It  is  considered  that  the  shadowing  can  be  disreyarded/neglected  up  to 
the  values  7~15  which  makes  it  possible  instead  of  precise  formulas  to 
use  maximum  expressions  for  7-“*. 


Schiff  [26]  proposed  a  number  of  simple  analytical  expressions 
for  section  taking  into  account  shadowing,  most  general/common  of 
them  takes  the  form 


(1.  30) 


i  t-  /  i±!  I  ^  -  -i-  Ini 
dk  137 *  \[l  1  Ufl  /  J  L  4  3 

2  <rl  l<f»;(T) _ L  ,  «|l 

f  4  3  ln  Z J)  ' 


In  the  case  of  complete  shadowing  7=0,  and  In  183  and  0^-4 

In  183  -  2/3,  then  formula  (1.30)  passes  in 


-  -5-^77 


(1-31) 


For  j<2  is  used  formula  (1.30),  and  values  <£,(-,)  and  0.(7)  are 
removed/taken  from  graphs  of  figure  4. 
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For  2<7<15  it  is  necessary  to  use  expression 
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moreover  function  C(7)  is  given  on  graph  of  figure 
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*•  u  -  '4 


Fig.  4.  Plotted  functions  ^£7)  and  ).  1  -  calculation  with 
precise  wave  functions;  2  -  calculation  according  to  Thomas -Fermi 
model . 
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Frequently  in  region  of  so-called  intermediate  shadowing  2<7<15 
is  used  more  exact  expression  of  differential  cross  section,  obtained 
by  Bethe  [27] 

S-  TB?  II1  +  1ST /  -  t  St]  w  '  1  - 


—f-.rct,!  !>]+£-,[£  l»«+6*>  + 

+  »"*«»-•»  +  t])’ 


(1.  33) 
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where 


b  - 


:i(io  t-l)  (-o-rl  -k)Z  '• 


111* 


cxnJ 


M 


i _ 1  *  j2 ,  i^.y 

(O)  L  2(io-f-l)  (i  !-l)  J  Mill/ 


In  the  case  of  complete  absence  of  shadowing  (y -*<*>)  both  formula 
(1.33)  and  (1.32)  they  lead  to  one  and  the  same  maximum  expression 


*1  =  tMf,  4-  /i±L\2  _  A  i±i  1  v 

d*  137*  [  '  \£«-(-l  J  3  £„+l  J  A 


x  [la  *l-+»  ■•*!>— 1], 


(1.  34) 


*  «  5  6  7  8  9  IQ 


Fig.  5.  Plotted  function  C(y). 
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Heitler  [28]  deduced  value  <p,  named  total  cross  section  and  defined  as 


the  section  of  energy  loss  to  the  emission 


(1.  35) 


In  the  approximation/approach,  used  by  Schiff, 


-  _  18  Z!rQ* 
’P~  3  137  * 


Segre  [29]  gives  the  following  formula  for 


(1.  36) 


* ‘-ar  [ln  2 


(1.37) 


Sections,  calculated  according  to  given  above  formulas,  compare 

with  experiment  [12'].  It  is  shown  that  only  for  the  wave  energies  of 
/tfeV 

electrons  2  m^rCEo<10  MeV  the  theoretical  sections  with  an  accuracy 
to  10%  will  be  coordinated  with  the  experiment.  In  extremely 
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relativistic  region  of  energies  these  formulas  give  the  overstated 
values,  while  in  the  region  of  energies  <2  MeV  -  those  understated. 

In  extremely  relativistic  region  of  energies  it  is  possible  to 
remove/take  limitation  of  Born  approximation  about  smallness  Z  and  to 
obtain  expressions  for  sections,  valid  for  any  Z. 


We  introduce  only  a  few  formulas,  which  are  analogs  of  formulas 
(1.30)-  1.34) , 

3 -W+OTPP-V -'**>] 


-t&PP-t -«*>])•  <i-38» 


where  f(Z)  = 


1,2021a1  anti  Majiux  Z 
0,925a2  an*  ^ojibniHX  Z. 


Key:  (1).  for  small  Z.  (2).  for  large  Z. 


For  complete  shadowing 


d «  _  r,« 

dk  187  * 


fill 
‘o+i  r 


(1.  39) 


For  p^a-<2  Oo/dk  is  computed  from  formula  (1.38).  Functions  <f\(T) 
and  are  found  from  the  graphs  of  Figure  4. 
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In  the  interval  2<7<15  formula 
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!  -  1  j 

,  .  -u  •  1)  (.,  •  ll 

(vi ) 

-  ■  1  ] 

111 

H 

-'r  -Ci 7)  -f(Z)  !.  Cl.  401 


is  used  and  C(i)  is  found  from  graph  of  figure  5. 


Expressions  for  total  cross  sections  take  form 

0  lafL|In1832'  14  fiZ>  .  11.41) 

and 

2  (*o  -H 1)  — |-  /<Z)].  (1.42) 

According  to  evaluations/estimates  of  Bethe,  these  formulas 
ensure  accuracy,  best,  than  2%  for  energies  E>50  MeV  and  Z~80. 
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Fig.  6.  Comparison  of  experimental  spectra  of  bremsstrahlung ,  which 
emerges  at  angle  8  (small  circles),  with  theoretical  (unbroken  curves) 
with  different  energies  of  electrons. 

Key:  (1).  sterad.  (2).  MeV. 
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Formulas  (1.30) -1.34)  can  be  improved  for  the  purpose  of  the  decrease 
of  the  limitation  of  that  superimposed  by  Born  approximation.  Usually 
this  is  achieved  by  introduction  to  formulas  (1.30)-  1.34)  corrections 
in  the  form  of  the  corresponding  factors.  This  correction  method  of 
sections  assumes  that  the  theory  correctly  transmits  distribution 
according  to  the  energy  of  the  photons  of  bremsstrahlung.  For  the 
high  energies  of  electrons  (E>4  MeV)  this,  probably,  it  is  correct. 

For  low  energies  there  is  a  noticeable  difference  in  the  form  of  the 
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spectrum  from  the  experiment.  Difference  is  especially  substantial 
near  the  maximum  photon  energy  (Fig.  6). 

For  nonrelat ivist ic  energies  of  Elwert  [30]  it  proposed 
correction  in  the  form 


[l  -exp 

IrZ  ’ 

|  (‘  1)1  («.|  1  -)  ‘u 

137  V  2) 

'  *  1  1 

1 -z  <  Hi  1 

ll 

!j  <•„  1)  )  (-  -2)i 

U.  43) 


According  to  evaluations/estimates,  during  use  of  this  correction 
it  is  possible  to  obtain  agreement  with  experiment  in  limits  of  10%, 
if  E, SO . 1  MeV .  The  region  of  energies  0. 1-2.0  MeV  is  most  difficult 
for  the  determination  of  corrections.  Here  there  are  no  theoretically 
obtained  corrections  whatever  and  usually  they  use  empirical.  They 
are  located  as  follows.  The  relation  of  the  total  cross  sections  of 
empirical  (experimental)  and  Born  taking  into  account  shadowing  as 
function  e0  and  Z  is  constructed. 
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Fig.  7.  Empirical  curves  (broken)  of  total  cross  sections  of 
bremsstrahlung .  Small  circles  -  results  of  experiment;  triangles  - 
calculation  in  extremely  relativistic  region;  unbroken  curves  - 
theoretical  dependences. 

Key:  (1).  Without  the  shadowing.  (2).  Empirical. 

Page  26. 

Empirical  values  are  obtained  by  mating  of  experimental 
dependence  <P  KC„  (s„,  Z)  with  theoretical  -  <PB(t 0,  Z)  in  extremely 
relativistic  region  (Fig.  7).  Berger  [31]  led  the  calculations  of 
relations  A  for  different  Z  and  energies  of  electrons.  The 

graph/diagrams  of  dependence  A  (Z,  e # )  are  shown  in  figure  8  and  can 
be  used  for  the  interpolation. 

As  can  be  seen  from  curves,  correction  is  maximum  with  energies 
c,=l  0.5  MeV)  independent  of  atomic  number. 
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Fig.  8.  Graphs  of  the  dependence  of  A(Z,  e.). 
Table  4.  The  corrected  sections  of  bremsstrahlung 

K^er*<ie-  (rK  He^Snpaa- 

cKaH  »Hep-  Mc^8jleHHUe  06J..CTB  upx- 

TPOK.  S0,  »«“««*•  Mep|J)0pMy.  MeHOHMM  H0CTt 

Afj*  Jiu)  a, 


0,10-2,0 
2.0  -15 

15-50 

50—500 

Key:  (1) 


d3**A  f  £ 
d  a — A  d  j  L 
A  da, 
A  da, 
da  — da, 

A  da, 
Ada, 
da— da, 

dii 

da, 


He^npaa- 
jieHHue  ce- 
nuHHH  (ho- 
Map  tpopMy- 
jiu)  a, 

Cl) 

OfjjiacTb  upif- 

MeHfeHltA 

M 

OaAaae- 
MftH  TO'I  ‘ 

KOCTfc 

1,28 

It  s  0,0l6„ 

r  20"., 

1,28 

r  -15 

— 

1,32 

2  t<  15 

1  5 

1,30 

7 

±5"u 

1,28 

1  >15 

— 

1,32 

2<’r<i5 

±  3'\, 

1,30 

T<2 

±3“u 

1,28 

li-15 

— 

1,40 

2<7<15 

±3 

1.38 

7^2 

C  3  X 

Key:  (1).  Kinetic  energy  of  electron  E. ,  MeV.  (2).  Corrected 

sections  da.  (3).  Defective  sections  (number  of  formula)  a,.  (4) 

Field  of  application.  (5).  Expected  accuracy. 
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In  Table  4,  borrowed  from  [12],  it  is  shown,  what  corrective 
multipliers  it  is  necessary  to  use  in  different  energy  ranges. 


Approximations/approaches,  utilized  for  calculations  of  section 
do  not  make  it  possible  to  consider  its  behavior  near  maximum  photon 
energy.  in  this  case,  as  more  correct  calculations  by  Fano,  etc. 
[32]  showed,  sections  were  different  from  zero  (in  the  Born 
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approximation  -  0)  and  satisfactorily  agrees  with  the  experiment, 


Since  the  discussion  deals  with  refinement  of  section  directly  at 
point  k=e0,  in  calculations  by  Monte  Carlo  these  corrections  can  be 
disregarded/neglected. 


In  conclusion  it  is  necessary  to  consider  one  additional  possible 

M  U  Ml  M  A  1  M  M  M  M  *  A  A  1  A  A  X  W  A  A  _  A  1  A  A  ^  A  a  M  «  A  «  BA  ^  A  A  A  ^  «  A  A  X  ^  A  A  A  A  A  A  L, 
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Page  28. 
Chapter  II. 


PROCESSES  OF  REPEATED  INTERACTION  OF  ELECTRONS. 

Sections  of  elementary  reports/events  of  interaction  of  electron 
with  atoms  of  substance  are  great,  and  therefore  even  in  small  route 
segments  in  substance  it  tests/experiences  large  collision  frequency. 
This  fact  makes  it  possible  with  the  high  accuracy  to  describe 
electron  motion  in  the  substance,  using  statistical  models.  In  this 
chapter  different  models  of  multiple  scattering  are  examined,  and  also 
a  question  about  the  average/mean  energy  losses  to  the  ionization  and 
the  bremsstrahlung  are  investigated  the  fluctuations  of  these  losses. 

S  1.  Average/mean  energy  losses  of  electrons  to  the  ionization  and 
the  bremsstrahlung. 


Known  that  average/mean  energy  loss  over  a  unit  of  distance  due 
to  any  process  is  determined  with  the  help  of  expression 


dE 

d* 


"•  ] 


(2.  1) 


where  Q  -  energy,  transmitted  in  single  report/event  of  interaction; 
do/dQ  -  section  of  this  transmission; 
n0  -  number  of  atoms  in  1  cm3  of  medium; 

and  <?■»«—  minimum  and  maximum  energy  transfers,  computed 


from  kinematic  relationships/ratios. 
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In  interesting  us  case  of  loss  energies  to  ionization  will  be 
registered  in  the  form 


\  w£odw' 

0 


(2.  2) 


moreover  E  is  measured  in  MeV,  and  path  in  see  integral  in  formula 
(2.2)  necessary  to  compute  separately  for  two  intervals  of  0<w<w'  and 
w’ <w<l/2 .  In  the  first  the  electron  is  considered  as  connected,  and 
integral,  according  to  Bethe  [33],  is  considered  with  the  help  of  the 
expression 


/  d*  .  ..  -  •  U*  f,  2 («.«.•*)* «4(«  i  2)  w'  2i  | 

\  “»  Tw  dw  ~  2icr»  .T nr  [ln  ~ - 7*  u  >  Ir  ]  •  (  * 


where  I  -  the  so-called  ionization  potential  (in  more  detail  about  it 
it  will  be  said  below). 

In  second  interval  for  do/dw  Moller's  formula  (1.20)  is  used. 
Storing/adding  up  both  integrals  and  taking  into  account  that  w'«l, 
we  will  obtain 


(i  £  ft 

u  /n  (>*-'  ~  2  Ho 


(i  r  1»-  |  1) 


[•  t  1)'  |  ,  i 

TT2T[ln 


ii i 


f  (*) 


where 


1  1  /  1  .  *  2trl  I  n  n 

f~  (e)  -  T  IT  (.  i'll  " 


.  (2.  4) 


(2.  5) 


Frequently  formula  (2.4)  is  called  Bethe-Bloch  formula. 


For  positrons  during  calculation  of  second  integral  is  used 
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formula  (1.26)  of  Bhabha  [22].  Expression  for  -dE/dx  is  obtained 
analogous  (2.4),  but  instead  of  f'(e)  it  is  necessary  to  substitute 
the  expression 


(s)  —  2  In  2 


14 

~r 


(2.  6) 


In  formulas  examined  was  not  considered  polarization  of  medium, 
which  appears  during  interaction  of  electron  with  atoms.  This 
phenomenon  known  even  by  the  name  of  the  "effect  of  density"  leads  to 
the  decrease  of  average/mean  ionization  losses. 


Page  30. 


Formula  for  calculating  the  average/mean  ionization  losses  and 
excitation  taking  into  account  the  effect  of  density  is 
written/recorded  in  the  form 


,.VZ  ,,  (m.c')** 

2*rQ2~-moC2  In  - 2y 


MO-*].  (2.  7) 


In  this  formula  energy  losses  are  expressed  in  MeV«cm:/g. 


Correction  for  effect  of  density  was  computed  in  work  of 

Sternheimer  [34].  The  formulas  obtained  by  him  for  6  take  the  form 

5  —  4,60b*  r  c  t  o  Ui  -  *)“  *0  • '  *  • ;  *! 

<>  —  4,606*  c  x  >  *,  (2.  8) 

6  ^  0  *  •  *o, 

where  x=lgp;  x„,  xlf  m,  a,  c  -  constants,  whose  values  for  different 
substances  are  given  in  Table  5. 


From  Table  5  it  follows  that  for  energy  E0<1  MeV  effect  of 
density  can  be  not  at  all  considered. 
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Ionization  potential  I  is  important  constant,  necessary  for 
calculations  according  to  formula  (2.7).  It  cannot  be  accurately 
calculated  by  contemporary  theory,  with  exception  of  hydrogen  and 
helium. 

According  to  Bloch,  for  heavy  elements  it  is  possible  to  use 
approximation  formula 

I  —  13,5  Z.  (2.  9) 

However,  the  experimentally  obtained  values  most  frequently  are  used. 
The  systematic  measurements  I  are  taken  by  Bakker  and  Segre  [35],  and 
also  Golovin  et  al.  [36].  These  data  are  cited  in  Table  5. 
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Table  5. 


U) 

ivl  a  i  u 
pa  a  a 

/. 

L36] 

[3«J 

! 

O’ 

m 

Li 

i 

1  34 

2 

-0,06 

3,07 

0.374 

3,06 

lie 

80,4 

64 

2 

-0,10- 

2,63 

0,413 

2,82 

C 

76,4 

60 

o 

-0,05 

3,22 

0,531 

2,63 

Al 

!  150 

163 

3 

0,06 

4.21 

0,0906 

3,51 

Vt 

1  243 

— 

3 

0,10 

4,62 

0,127 

3,29 

Cu 

279 

314 

3 

0,20 

4,74 

!  0,119 

3,36 

A* 

422 

— 

3 

0,20 

5,75 

1  0,251 

2,86 

Hn 

472 

3 

0,20 

6,26 

!  0,404 

2,52 

W 

660 

4 

0,30 

6,03 

!  0,0263 

3,91 

Pb 

737 

1  766 

4 

0,40 

6,93 

0,0662 

3,41 

U 

|  663 

1 

4 

0,30 

6,69 

0,0652 

3,37 

ICH.), 

1  62 

I 

•  > 

1  “ 

0,12 

2,94 

!  0,393 

i  i 

2,66 

Key: 

(1)  . 

Mater 

ial . 

(2)  . 

eV. 
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By  another  source  of  energy  losses  fast  electron,  moving  in 
substance,  are  losses  to  bremsstrahlung .  Expression  for  the 
average/mean  losses  over  a  unit  of  distance  in  this  case  can  be 


obtained  with  the  help  of  the  formula 

C 

1  dE 


da 


dk. 


(2.  10) 


x  ue  y  7  i  , 

7dI^~Am^C  \Ldk 

0 

The  integral,  which  stands  in  right  side  (2.10),  according  to  (1.35) 
coincides  (e„+l)  o>,  where  </>  —  the  total  cross  section  of 
bremsstrahlung.  Therefore  average/mean  losses  can  be  easily 
calculated,  if  are  known  formulas  for  <P. 


For  e>>l  they  use  following  formula 


i_  UE 
dx 


1) In  [2(«„ -1-1) (2.11) 


which  is  obtained  during  use  of  formula  (1.3"’).  We  give  in  “fables  3 
and  4  applications/appendices  of  the  value  of  ionizing  losses 
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i  a  i: 


and  total  losses 


i  je 


d  '  1  *  designed  with  the  help  of 


formulas  (2.7)  and  (2.11).  These  values  can  be  used  for  the  rapid 
rough  estimates  of  the  losses  of  energy  AE  in  the  small  route  segments 
Ak4v  xr 


(2.  12) 


Frequently  results  of  experiment  or  calculations  according  to 
passage  of  electrons  through  substance  it  is  possible  to  describe  by 
universal  by  curves  or  formulas,  if  distance,  passed  of  electrons,  is 
expressed  in  units  of  complete  landing  run.  By  complete  landing  run 
is  understood  the  maximum  distance,  passed  by  electron,  if  it  moved 
along  straight  line  and  it  lost  energy  continuously. 


Complete  landing  run  in  g/cm*  is  computed  from  formula 


where 


1  dE 
p  dx 


=  /_JL  ££\ 

„  \  P  ax)v  ~  [  f  dx  h 


(2.  13) 
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The  calculations  of  the  complete  landing  runs  in  the  wide  region  of 
energies  from  10  keV  to  1000  MeV  for  different  substances  are  carried 
out  by  Berger  and  Seltzer  (37].  Since  these  calculations  are 
published  in  the  special  report,  we  give  in  Table  5  of  the  Appendix 
values  of  landing  runs  for  several  elements  in  the  region  of  energies 
of  electrons  0.2-10  MeV.  For  determination  R0(E„,  Z)  (Z>6)  the 
tabulated  data  can  be  interpolated. 
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S  2.  Fluctuations  in  energy  losses. 


In  previous  paragraph  we  considered  average/mean  energy  losses  of 
electrons.  In  actuality  electron  loses  energy  discretely,  and  are 
possible  losses  both  large  and  it  is  smaller  than  the  average.  For 
the  first  time  a  question  about  the  fluctuations  in  energy  losses  to 
the  ionization  and  the  excitation  was  examined  by  Landau  [38]. 


Function  of  energy-loss  distribution  in  transit  through  thin 
layer  of  substance  is  obtained  by  Landau  under  the  assumption  that  a) 
losses  to  bremsstrahlung  are  small,  b)  probable  loss  of  energy  \£„ 
satisfies  inequalities  A£M  £0  and  A  a.,  —  electron-binding 

energy  in  atom,  /)• 


In  Landau’s  theory  is  not  considered  correction  for  effect  of 
density  and  that  it  is  more  substantially,  is  not  taken  into 
consideration  effect  of  multiple  scattering  of  electrons  in  this 
layer.  According  to  Yang  [39],  the  account  of  the  latter  actually  is 
reduced  to  the  replacement  of  thickness  of  the  layer  At  to  efficient 
thickness  .\c  ,+ll;  \<). 


By  most  correct  method  account  of  multiple  scattering  can  be 
carried  out  with  the  help  of  calculation  by  Monte  Carlo. 


Attractiveness  of  Landau's  method  consists  in  the  fact  that 


V  / 
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probability  of  energy  loss  in  the  range  AE-  AE+d(AE)  succeeds  in 
expressing  through  universal  function  wL0),  where  X  -  dimensionless 
parameter. 

w  (A£)  d  (AS)  wL  (/■)  d > ,  (2.14) 


cln-.i^i'rax 
\E-ax  ^ In  ,=77r 


0,123 


(2.  15) 


where  a=0.153  Z/A/32  (MeV^cmVg),  x  -  thickness  of  the  layer  in  g/cmJ . 
The  detailed  tables  of  the  values  of  function  «,.('>  are  given  in  the 
work  of  Bo^rsh-Supan  [40]. 
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For  value  X=-0.05  function  u>L (/)  has  maximum.  Energy  losses, 
which  correspond  to  the  maximum  of  distribution,  is  conventionally 
designated  as  the  probable  energy  losses.  From  (2.15)  with  X=-0.05  we 
obtain  formula  for  the  probable  losses 

AE.  =  ax  Jin  -  p*  +  0,373  ] .  (2.  16) 

Knowing  the  function  of  the  distribution  of  losses,  it  is  easy  to  find 
the  average/mean  energy  loss  in  the  form 

-  ax  Jin  —  S' J .  (2.17) 

It  is  important  to  note  that  the  probable  energy  loss  does  not 
coincide  with  the  average. 


From  asymmetry  of  curve  wL( *)  it  follows  that  for  small  x 

average/mean  losses  are  more  than  probable,  and  this  is  connected  with 
noticeable  probability  of  large  losses  of  energy  Afi  A E,.  However*^ 
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with  an  increase  in  the  thickness  of  absorber  Atf.  44.  increases  more 
rapidly  than  AE  due  to  the  further  logarithmic  dependence  on  x  and 
with  the  the  large  x(ax>E/2)  distribution  curve  becomes  close  to  the 
Gaussian,  while  A~£«AE„. 

In  calculations  by  Monte  Carle  it  is  necessary  to  know  function 
of  energy-loss  distribution,  greater  than  given,  i.e.,  integral 
function.  The  probability  of  energy  loss,  which  exceed  AE,  is 
determined  with  rhe  help  of  the  integral 

aam»x 

*  (A£)  d  (A£)  =  d  (A£)  j'  u>(A£')d(A£').  (2.18) 

A  E 

Landau's  theory  was  improved  in  works  of  Vavilov  [41],  Symon 
[42],  Blunck  and  Leisegang  [43].  The  distributions  of  losses 
according  to  Vavilov  are  tabulated  in  the  report  of  Seltzer  and  Berge 
[44]  and  more  frequently  they  are  used  for  the  heavy  charged/loaded 
particles.  For  Symon  it  was  possible  to  consider  the  case,  when 
energy  losses  were  compared  with  the  wave  energy  of  electron,  i.e., 
when  the  thickness  of  absorber  no  longer  was  small.  The  distribution 
function  in  this  case  is  represented  by  the  set  of  curves  with 
different  degree  of  the  asymmetry:  from  the  gaussian  form  with  the 
large  thickness  to  the  distribution  of  Landau  with  the  small 
thickness.  Very  work  of  Symon  is  published  in  the  form  of  theses; 
therefore  it  is  virtually  unavailable,  but  the  results  of  its 
calculations  in  the  form  of  graphs  and  explanations  to  them  are  given 
in  Rossi's  book  [45]. 
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Energy-loss  distribution  on  Symon  is  given  as  the  function  of 

dimensionless  parameter  - r — -  ,  where  AE  -  energy  losses,  and  A E, 

probable  losses,  computed  from  the  formula 


A —  ax  In 


2m  ,i3  ax 

rr~wf: 


A0  -=-axb. 


(2.  19) 


Family  of  graphs  u>.  (1)  is  shown  in  figure  9a,  while  family  of 
integral  function  'M'>  —  in  figure  9b.  In  order  to  use  these  graphs, 

UX 

it  is  necessary  to  at  first  compute  value  where  tlm —  maximum 

energy,  transmitted  to  the  electron  (it  is  computed  from  formula 
1.16). 


It  is  further  necessary  to  find  parameters  j,  b,  X,  depicted  on 
graphs  of  figure  10  as  functions  G.  In  conclusion  is  computed  value 
£,  and  from  it  and  u>,  (i)  (in  Fig.  9a  in  right-hand  corner  is  shown 
normalizing  factor  F,  which  is  function  X).  For  convenience  in  the 
use  of  the  distribution  of  Symon  we  give  6  applications/appendices  of 

— - — -j.  in  the  Table- 
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Fig.  9a.  Graphs  of  the  energy-loss  distributions  w  <:i'  on  Symon. 
Numerals  in  the  figures  -  value  of  the  parameter  >■,  -  normalizing 

factor . 
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Blunck  and  Leisegang  taught  resonance  interaction  of 
falling/incident  electron  with  inner  electrons  of  atom.  The  function 
of  the  energy-loss  distribution  in  this  case  is  obtained  somewhat 
wider  than  the  corresponding  functions  of  Landau  or  Symon.  The 
distribution  of  Blunck-Leisegang  is  family  of  curves,  depending  on 
parameter  b* 

b2  —  q&E  Z‘  >(ax)~2,  (2.20) 

where  qae20  eV. 

Very  distribution  function  is  sum  of  four  gaussian  curves  and  is 


written/recorded  in  the  form 
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n  (Afc')ti  (A£)—  V  - •  l  '  exp - — —  dr  ,  (2.21) 


V  r  y 


where  X  -  dimensionless  parameter,  determined  by  expression 


--  |  In  -  -1,116, 

ax  1  ax 


(2.  22) 


and  constants  c,  are  given  in  Table  6. 


Let  us  note  that  cases  of  large  energy  losses  in  distribut: 


(2.21)  are  not  examined. 
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Fig.  9b.  Graphs  of  integral  distributions  *  on  Symon. 

Page  36. 

When  b2«3,  theory  of  Blunck  and  Leisegang  brings  to  coinciding 
with  theory  Landau  to  distribution.  The  role  of  resonance  broadening 
becomes  noticeable  only  with  the  the  large  b2,  i.e.,  for  the 
relatively  small  thicknesses. 

By  deficiency  in  distribution  in  question  in  calculations  Monte 
Carlo  is  use  of  considerable  space  of  working  storage  of  computer (s) 
for  storing  family  of  curves. 
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Fig.  10.  Graph/diagrams  of  dependence:  a)  j(G),  b)  b(G),  c)  X(G). 
Curves  1,  2,  3,  4  correspond  to  values  02= 0;  0.4;  0.7;  1. 
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Early  experiments  [46]  agreed  well  with  Landau's  theory,  if  we  i 
latter  consider  "effect  of  density"  and  to  take  into  consideration 
multiple  scattering  and  bremsstrahlung .  In  later  works  [47,  48]  it 
was  noted  that  during  the  use  of  instruments  with-  the  best  energy 
permission/resolution  for  the  agreement  of  experiment  with  the  theory 
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it  is  necessary  to  calculate  the  resonance  broadening  of  distributions 
on  the  basis  of  Blunck  and  Leisegang. 


Energy-loss  distribution  for  positrons  was  for  the  first  time 
examined  in  work  Rohrlich  and  Carlson  [49].  Landau's  method  was  used, 
and  for  the  scattering  cross  section  of  positron  on  the  electron  was 
used  formula  (1.26)  of  Bhabha.  For  the  probable  energy  loss  in  this 
work  the  expression 

=  ax  [in  -  ?2  +  0,373  -  2,85] , 

where 


is  obtained 


,  2netnll  Zx 
0  ilmgc*)1 


(«+2)‘ 


(2.  23) 


Nishi  and  Hiromasa  [50]  carried  out  experimental  check 
energy-loss  distributions  of  positrons  with  energy  0.3  MeV  in 
different  substances,  indicate  that  theory  [49]  is  valid  only  for 
small  thicknesses  and  that  with  increase  in  thickness  experimental 
distribution  very  well  will  be  coordinated  with  distribution  (2.21), 
if  average/mean  energy  losses  are  computed  on  (2.7)  with  f+(e).  Thus, 
distribution  (2.21),  which  more  fully  considers  all  effects,  connected 
with  the  fluctuations  of  energy  loss,  can  be  used  also  for  the 
positrons.  In  the  calculations  by  Monte  Carlo,  where,  as  a  rule  are 
examined  thin  layers,  for  the  positrons  will  be  valid  and  the  theory 
of  Symon. 

To  fluctuations  are  subjected  other  losses  of  energy  -  radiation 
losses.  The  function  of  the  distribution  of  these  losses  they  for  the 
first  time  computed  Bethe  and  Heitler  [25]. 


DOC  =  86063203 


PAGE 


&7 


Table  6. 


1 

o 

3 

4 

c  0,174 

O.OM 

0,019 

0,007 

a  0,000 

3,000 

8,600 

11,000 

1.800 

2,000 

3,000 

5,000 
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They  proceeded  from  simplified  formula  (1.31)  for  the  differential 
cross  section  of  bremsstrahlung .  Expression  for  the  probability  of 
the  losses  of  energy  AEr  takes  the  form 


Eu  1»«- 


lr  kF  fH(.h  d(S£)  I  Eu-£E 
<o(E0,  AE,  l)d(AE)-  -g—  - 


,  (2.24) 


•yt  4  ii  cm*  . 

where  a  1,4  10*--  -j-  In  183  Z~‘H-  -vp  (  —  j*  a  f(xx)  -  -  gamma  function. 


Distribution  cj(E0,  AE,  t)  is  calibrated  to  one.  The  form  of  the 
function  is  such,  that  when  ax>l,  more  probable  are  the  large  losses 
as  a  result  of  a  few  reports/events  of  interaction,  and  when  ax<<l , 
more  probable  become  the  low  losses  of  energy.  This  latter/last  case 
most  frequently  is  realized  in  the  calculations  by  Monte  Carlo.  Until 
now,  we  examined  both  sources  energy  losses  independently. 


It  is  more  right,  as  this  is  shown  by  Blunck  and  Westphal  [51], 
to  seek  function  of  distribution  of  total  losses  due  to  both 
processes.  It  is  necessary,  consequently,  to  find  the  probability  of 
the  given  losses  knowing  the  functions  of  the  distribution  of 

losses  tc  the  bremsstrahlung  and  to  the  ionization.  The  unknown 
distribution  is  located  with  the  help  of  the  convolution  integral  in 
this  form: 
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where 


u>(A£„)  \  ui(A£\  d'„  (2.251 

0 


A£u-  A £n-v. 


In  work  [51]  families  of  these  distributions  are  designed,  using 
as  w\sE„  —  >,)  distribution  (2.21)  and  »*Ci)  in  the  form  (2.24).  We  led 
analogous  calculations,  but  the  energy-loss  distribution  to  the 
ionization  was  undertaken  on  Symon.  The  graphs  of  the  integral 
distributions 


£, 

*(A En)  ^  \  w  (A£„)  d.  (A£.,)  (2.26) 

0 

are  shown  in  figure  11. 


On  figure  it  is  evident  that  for  Pb(E0=5.75  MeV)  role  of 
fluctuations  in  radiation  losses  becomes  essential  (for  comparison  in 
the  same  figure  shown  distribution,  which  considers  only  fluctuations  - 
in  ionizing  losses).  For  Si(Eo=0.95  MeV)  distribution  (2.26) 
coincides  with  the  distribution,  designed  only  for  the  ionizing 
losses . 
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Let  us  note  that  even  in  the  case  of  small  thicknesses  and 
relatively  low  energies  E~2  MeV  average/mean  energy  losses,  calculated 
for  Z>50  during  use  of  distribution  (2.24),  are  more  than  sum  of 
average/mean  energy  losses,  calculated  separately  from  each  other, 
which  indicates  importance  of  account  of  fluctuation  of  total  losses. 
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Checking  distribution  (2.25)  showed  that  it  somewhat  wider  than  the 
experimental  [52,  53].  This  is  connected  with  the  approximate 
character  of  expression  for  w (rj).  It  is  probable  that  distribution  of 
Symon  in  the  combination  with  q(t?)  will  give  better  agreement  with  the 
experiment,  since  it  is  somewhat  narrower  than  the  distribution  of 
Blunck  and  Leisegang.  However,  the  further  experimental  check  of  this 
assumption  is  required. 

§  3.  Multiple  scattering  of  electrons. 

Results,  obtained  in  calculations  by  Monte  Carlo,  are  to  a 
considerable  extent  connected  with  detailed  and  checked  during 
experiment  theories  of  multiple  scattering.  We  will  consider  here 
only  two  distributions  for  the  angles  in  multiple  scattering: 
Goudsmit-Saunderson  and  Mol i ere. 


I -,*)  * 
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Fig.  il.  Graphs  of  integral  function  of  distribution  of  energy 

loss  of  electron.  1  -  total  losses;  2  -  loss  only  to  the  ionization 
(E0=5 . 75  MeV  in  Pb) .  For  Si  and  Eo=0.95  MeV  both  distributions 
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1.  Goudsmit  and  Saunderson  [7]  showed  that  the  quantity  of 
electrons  scattered  in  interval  of  angles  0-0+d0  after  the  passage  of 
path  segment  t,  is  given  by  formula 

x 

f  r-c  (h  fj)sin  OdO  —  sin'JdO  V  (/  —  1  _.)  Pl  (cos  0) 

i  0 


where 


X  exp 


\Gt(t')dt 

o 


(2.  27) 


Gi(t)  =  na  |  [1—  P,  (cos  0)]u<J.  (2.  28) 

o 

Here  (do  ii».  o  /dQ  -  section  of  single  scattering,  p,(cos>»)  -  Legendre's 

polynomial.  The  decrease  of  energy  of  particle  due  to  the  ionization 
losses  and  bremsstrahlung  in  distribution  (2.27)  is  considered  through 
the  implicit  functional  dependence  t=t(E),  moreover  the  continuous 
losses  of  energy  are  assumed.  In  such  a  manner  both  in  the 
distribution  for  energy  losses  does  not  succeed  in  taking  into  account' 
multiple  scattering  and  in  the  distribution  for  the  angles  in  repeated 
scattering  it  is  impossible  to  correctly  consider  fluctuations  in 
energy  losses. 


To  advantages  of  distribution  in  question  should  be  related, 
first,  possibility  of  using  any  expressions  for  section  of  single 
scattering,  secondly,  applicability  of  it  for  any  scattering  angles, 
including  large  (~180°).  Use  in  (2.27)  of  the  Mott  cross  sections 
(1.2)  makes  it  possible  to  consider  the  relativistic  and  spin  effects, 


DOC  =  86063204 


PAGE 


essential  for  the  large  scattering  angles,  which  lead  to  the 
difference  in  the  passage  of  electrons  and  positrons  through  the 
substance  into  consideration. 


Deficiency  in  theory  of  Goudsmit-Saunderson  is  weak  convergence 
of  series  into  (l-ll)  especially  for  small  t  and  large  0,  as  a  result 
of  which  for  obtaining  necessary  accuracy  it  is  necessary  to  add  large 
number  of  terms  of  expansion  (to  100).  The  latter  so  impeded  the  use 
of  distribution  that  it  in  practice  was  not  used  in  the  calculations, 
on  it  served  as  the  criterion  of  the  validity  of  one  or  the  other 
simplified  variations  of  the  theory  of  multiple  scattering. 

Page  41. 

Only  after  Spencer  proposed  the  recursion  relations,  which  make  it 
possible  to  compute  Gt  with  the  high  indices  up  to  /=100,  obtaining 
distribution  (2.27)  became  relatively  simple  procedure. 


Examination  of  Spencer's  method  [10]  let  us  begin  from 
calculation  of  integral 


/=  f  G^Odt'. 


(2.  29) 


Let  us  switch  over  to  the  new  variable  s=(R„-t)/(Re ) ,  where  R0  is  the 
total  path  of  electron  in  the  substance.  Then 

t 


/  -=  f  Gt  (t')dt'  -  H0  j  G,*  (•')<*«'•  (2-  a0> 


According  to  Spencer, 
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G,*(*)-Gi(*>  §4^ 


(2.  31) 


where 


<V0>)  '  7/ 


a(«l  C,)’ 


(2.  32) 


and  G*  (1)  and  G^CL)  coincide  with  the  values  of  function  G,  of  the 


corresponding  order,  calculated  for  the  wave  energy  of  electron. 


Substituting  (2.32)  and  (2.31)  in  (2.30),  it  is  possible  unknown 


integral  (2.29)  to  register  in  the  form 


r„c ,  G<«(  1).  « --C; 

C J  *G,vUr  *(1tC2)  • 


(2.  33) 


Parameters  Cx  and  C2  can  find,  solving  system  of  equations  with  the 


the  known  Gi*(s)  and  G»*(l): 


°»*w-7(rte 


^•W-iSv 


(2.  34) 


Let  us  consider  the  case,  when  t/R,<<l,  i.e.,  s«l ,  being  usually 


realized  in  the  Monte  Carlo  calculations.  Then,  after  being  bounded  to 

In 

the  first  term  of  expansion  in  the  Maclaurin  series,  after 

a  (  l  i  C  a) 


simple  transformations  we  obtain  convenient  approximation/approach  for 


I  -G,(t  ->0)  t. 


(2.  35) 
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Let  us  find  now  expression  for  Gf.  using  as  section  of  single 


scattering  Mott's  formula  (1.11)  with  correction  term,  described  by 


formula  (1.12).  After  passing  to  new  variable  u=cos0  and  after 


.V  V  /■  /.  j-.  .*  -• 


'.■<vVW 
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fulfilling  in  (2.28)  the  integration  for  the  azimuthal  angle,  we  will 
obtain 

a‘  1-  ?-  z(/  i)  (i  2,  />  j 

i  1 1  - 

-  A.jg,  i,/)  ti .  /, 

l)  ■  haQ  (2.  ib) 

The  confronting  in  curly  braces  (2.36)  polynomials  Q(m,  /)  are 
expressed  as  polynomials  (l  u  2‘,)n  with  the  help  of  the 
relat ionship/rat io 

i 

Q[‘n,  l)  j(l  u  |  p,(u)\du,  (2.37) 

-i 

where  *i  -  parameter  of  shadowing. 

On  the  basis  of  general/common  properties  of  polynomials  of 
Legendre,  Spencer  gives  following  recursion  relations  for  Q(m,  l ) . 

The  first  makes  it  possible  to  pass  from  m  to  m+1: 

Q(m  ,  1,1)  (1  2 \)Q(ik,  l)  -J  Q(m,  1) 

•  l  ~\  Q(m,  ID-  2lLiQ(m,  l  1).  (2.38) 

moreover  Q(m,  0  )e0. 


The  second  makes  it  possible  to  pass  from  l  to  /+1  (for  this  m) , 
moreover,  if  m  -  whole,  then  is  used  (2.39),  and  if  m  -  fractional  - 
(2.40) : 

IQ(  2,  l  t  1)  -  (2/  1)  (1  T  2/,)Q(  -  2,  l)  —  (7  i-  1)Q(  2, 

l  -  1)  -  (21  1)  (1  ~  >,)->,  (2.  3D) 
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i  n 


1 1  i  q\ 


1  I  .  (2.  40) 


where 


i  1,  a  -  1  2  t, 
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The  initial  values  Q(-2.1)  and  Q(-(3/2),  1)  necessary  for  using  the 
formulae  above  are  computed  from  the  formulas 


Q(-2,  1)  =  In  (1-f-fj- 


(1  -i  ',) 

Q  -- f-,  l)  —  2(2^)’ J(1  r 

+ ~,rl-  (2. 4i) 

We  already  indicated  weak  convergence  of  series.  Therefore  for 
obtaining  the  necessary  accuracy,  especially  in  the  region  of  large 
angles  6,  it  is  necessary  to  count  to  100  members  of  series/row 
(2.27) . 

Table  7  shows  convergence  of  values  fie,  t),  calculation  is 
carried  out  for  A1 ,  0=1  MeV,  t=0.03  g/cm*. 


Moliere's  theory  [8]  is  simpler  theory  of  multiple  scattering, 
strictly  speaking,  applied  for  small  scattering  angles  (0<1  rad).  The 
latter  made  it  possible  to  create  the  completed  theory  and  to  obtain 
the  distribution  function  for  the  scattering  angles  for  the  scattering 


m 
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angles  in  the  analytical  form. 
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Table  7.  Convergence  of  series  for  fr  c  n  as  the  function  of  a 
number  of  terms  l . 


'V  '<-)• 

I  \l 

'>  =  18 

45 

72' 

[0 

! 

135 

!  1 

17] 

10 

3.2723697 

0,52174434 

i  —0,46401429 

0,39853767 

—0,23270203 

—0.10793879 

20  : 

3,6323317 

0.17457685 

—0,06470117 

—0,05865163 

0.03858396  ! 

0.01550726 

30 

3,4721717 

0,06237644 

0,00305118 

0,01081927 

0,,. 0403899 

0.00173254 

40 

3.466925;; 

C. 06750488 

0,01025154 

0,00349300 

—0.0000-652  | 

-0.00016018 

50  ; 

3,4688676 

0,06871539 

0,01097318 

0,00424696 

0.00026294 

— 0,00000999 

60 

3.4687265 

0,06865563 

0,01105342 

0.00415893 

0,00033319  . 

0,00001262 

70 

3.4686945 

0,0686.3563  l 

0,01108404 

0,00417140 

0,00032769  1 

0.00001057 

80 

3.4686932 

0,06863721 

0,01106576 

0,00116921 

0.00032640  1 

0,0000099.' 

90  1 

3,4686946 

0.06863796  1 

0,01106610 

0,00416969 

0,0003265*  | 

0.00001006 

100  ! 

1 

3,4686246 

0,0686378  •  j 

0,01106618 

1 

0,00416956 

0.00032665  j 

0,00031010 
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Thoroughly  analyzing  theory  of  Moliere,  Bethe  [54]  showed  that  it 
is  possible  to  introduce  corrective  factor  into  distribution  function, 
making  it  possible  to  spread  theory  to  large  scattering  angles  and 


thus  to  make  its  universal  just  as  theory  of  Goudsmit-Saunderson. 
Below  we  still  will  return  to  the  comparison  of  results  of  both 
theories,  and  now  we  will  proceed  to  the  examination  of  the 
fundamentals  of  Moliere's  theory. 


The  presence  of  only  one  parameter,  so-called  angle  of  shadowing 
(1.7)  is  most  important  result  of  theory,  which  makes  it  attractive 
for  experimenters.  Further,  the  distribution  function  does  not  depend 
on  the  form  of  the  section  of  single  scattering;  only  requirement  for 
the  section  so  that  it  would  pass  at  the  large  scattering  angles  in 
the  section  of  Rutherford  (1.1).  Let  us  note  immediately  that  this 
assumption,  as  it  follows  from  the  discussion,  carried  out  in  Chapter 
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I,  in  reality  occurs.  For  the  probability  of  electron  scatterings  in 
interval  of  d-d+dd  after  the  passage  of  layer  t  Moliere  gives  the 
expression 

f*  (0.  t)»d0^  [r*  (»)  - !  -  B-‘  r ■> (0)  -r  b  -  fu'  (U)  - 

...):>rfl>.  (2.12) 


Condition  for  applicability  (2.42)  to  small  angles  is  satisfied 
by  replacement  of  sin0«0.  For  obtaining  the  accuracy  of  order  1%  it 
suffices  to  consider  in  (2.42)  only  three  first  terms. 


F*  («>)  -  Gaussian  type  function,  and  functions  f1  ii>)  and  fl(j>) 

consider  scattering  of  large  angles  and  divergence  from  Born 
approximation.  The  numerical  values  of  these  functions,  refined  by 
Bethe,  are  given  in  the  Table  of  7  applications/appendices, 
dimensionless  variable  «►  is  connected  with  the  scattering  angle  6  and 
the  parameter  of  theory  B  with  the  relationship/ratio 


where 


t  /<.  *  & 


(2.  42) 


4- nutt*Z(Z~  1) 
Zc*  Pl  c- 


(2.  44) 


In  formula  (2.44)  and  in  those  following  it,  the  correction  for 
mentioned  in  Chapter  I  inelastic  scattering  of  electrons  Z*->Z(Z+1)  is 
taken  into  consideration.  The  new  parameter  B  is  obtained  from  the 
solution  of  the  transcendental  equation 
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where 


ln  1.187  /u- 


.  8680  M  Z  1 1  Z  * 

n  ;i-  a  ( i  ).:ii 


(2.  lb) 


it  is  expressed  as  the  angle  of  shadowing  /u-'~=4-r,  and  is  described 
by  formula  (1.7).  Values  B  usually  lie/rest  at  interval  of  5-20. 
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Let  us  note  that  value  gives  exemplary/approximate 

collision  frequency,  experienced/tested  by  electron  with  passage  of 
thickness  of  the  layer  t.  Since  by  its  nature  Moliere’s  theory 
statistical,  then  collision  frequency  must  be  large,  is  usually  more 
than  20.  This  is  one  of  the  limitations  of  theory.  On  the  other 
hand,  collision  frequency  and,  consequently,  the  thickness  of  the 
layer  cannot  be  large  so  that  the  distribution  would  not  be  distorted 
due  to  the  disregard  of  energy  loss  to  the  ionization. 


Bethe  in  work  [54]  showed  that  if  we  use  in  both  theories  of 
Moliere  and  Goudsmit-Sauderson  the  Rutherford  scattering  cross 
section,  then  between  appropriate  distributions  there  is  dependence 

tv  c(®,  0  -  U,"  „y',cxp[-1V  *>  L  •  <2-47) 

This  formula  is  valid  as  long  as  •  U  <£ l,  thus  far  the  width  of  the 
gaussian  part  of  the  distribution  of  Moliere  SI  rad.  But  in  this  case 
l8  /,-  B J  -1.  Without  the  large  error  it  is  possible  to  disregard  in 

(2.47)  and  1/24.  Then  the  improved  distribution  function,  spread  to 
the  region  of  the  angles  of  0>1  rad,  is  located  with  the  help  of  the 


formula 
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fs  (0)  sin  OdO  =  OdO  "  •• 

t.c  V  B  » 

+  +  (2.48) 

Obviously,  correction  factor  'J* 0)  '  although  improves 

distribution  of  Moliere,  it  cannot  with  full  weight  contain  behavior 
of  function  for  large  scattering  angles.  Is  more  than  that, 
scattering  cross  section  to  the  large  angles  is  sensitiv  to  the 
shadowing,  and  the  atom  model  according  to  Thomas-Fermi ,  which  is  used 
in  Moliere's  theory,  it  does  not  transmit  this  dependence.  In  Figure 
12  we  are  congruent/equate  distributions  (2.48)  and  (2.27).  The 
latter  was  calculated,  using  a  section  of  Mott. 

From  figure  it  is  evident  that  in  region  of  large  angles  for  Au 
difference  between  both  distributions  to  become  essential.  One  should 
indicate  one  more  deficiency  in  Moliere's  theory  -  in  it  the 
difference  in  the  angular  distributions  for  the  electrons  and  the 
positrons  is  not  exhibited. 
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Development  of  Moliere's  theory  is  given  in  works  of  Nigam  et  al. 
[9].  First  of  all  it  seemed  that  the  shielded  scattering  cross 
section,  obtained  by  Moliere  during  the  use  for  the  atom  of 
Thomas-Fermi  model,  is  imprecise  in  all  orders  a=  Z/137,  except  zero. 
Correct  expression  for  the  section  can  be  obtained  by  the  method, 
proposed  by  [55]. 
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Nigam  et  al.  use  relativistic  scattering  cross  section  of  Dalitz 
[56]  with  shielded  Coulomb  potential.  Although  this  potential  is 
worse  than  Thomas-Fermi  potential,  the  procedure  of  calculation  of  the 
distribution  function  according  to  Nigam  is  such,  which  ensures 
coincidence  with  a  precise  distribution  of  Goudsmit-Saunderson . 


Angle  of  shadowing  according  to  Nigam  4Vv.  where  is  given 


by  formula  (1.8).  The  difference  between  the  angles  of  shadowing 
according  to  Moliere  and  Nigam  sufficiently  substantially  for  large  2 
and  low  energies  is  shown  in  Table  1.  Since  into  the  formula  .for  the 
angle  of  shadowing  (1.8)  a  enters  to  the  first  degree  (a2  -  into  the 
theory  of  Moliere),  in  the  distribution  function  is  exhibited  the 
difference  between  the  electrons  and  positrons  [57],  not  observed  in 
Moliere's  theory. 


▼ 


rjri  ir»v  :  y,,.nv7T'"ir-'  w?fv»»^7u 
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Fig.  12.  Comparison  of  the  function  of  distribution  f(0)sin@, 
calculated  for  A1  and  Pb  Eo-0.2  MeV:  1  -  according  to  Moliere;  2  - 
according  to  Goudsmit-Saunderson . 
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Function  of  distribution  was  obtained  in  second  Born 
approximation/approach  can  be  represented  in  the  form,  analogous  to 
distribution  of  Moliere 

f(,J,  t)  -  77a  [r0' +  -£-(/' n'  t/ui)  h  if'*''  -r  P')] ,  (2.49) 

where  f<#\  fa'  and  f2'  -  Mol iere-Bethe  functions  (see  Table  7 
application/appendix  ) ;  fu’  1  and  f^  ’  -  corrections,  which  appear  due 
to  second  Born  approximation 


fll)  —  J  f  (u)[—  V  B  u]du 


(2.  50) 


and 
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r  U) 


!  U  In  u-j  du,  (2.51) 


where 


'f  U)  -  exp  1+4* /c-  (  +  +  - -  In  u 


Jo 


>»U 


l  U  I  B 


exp 


U- M 


(2.  52) 


J0  -  Bessel  function,  and  yc  is  found  from  formula  (2.44).  Other 
parameters,  necessary  for  calculation  fL1)  ’  and  ffJ^' ,  are  determined  by 
the  expressions 


+  11  l,34-10-5~(l-p-)\ 

? 


r ---  ^4-^-  i|x(  V  b. 


Xu  ^  2 

h 


(2.  53) 


/,  - 1.26. 


For  the  coefficient  K  of  Nigam  it  gives  the  expression 


*-  «*l>l+B/.=(l  i’TT 


r  —  (+  i  ~a{s)  [0,5772  -  In  (y.c  [  £)Jjj  .  (2.  54) 
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Parameter  B  is  found  from  the  solution  of  the  transcendental  equation 


B-llnB  b,  (2.55) 

where 

and  X,-  is  computed  from  the  formula 
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In  j.?_(  =  in  _  0.8456| —  j  (1  -  V)  -  0.0772.  (2.56) 

Dimensionless  variable  in  theory  of  Nigam,  just  as  in  Moliere's 

theory,  is  value  o Marion  and  Zimmerman  [58]  they  considered 

u-y  b 

the  simplified  variation  of  the  theory  of  Nigam.  Actually,  in  the 
region  of  energies  0.4-6. 0  MeV  examined/considered  by  us  for  all 
elements,  including  Z=82,  £=*1 ,  /.«•  '  /  ,  K31! ,  /«  ~  10  ^hen  terms  f'1'  ’ 

and  f  '-1'  can  be  disregarded/neglected  respectively  in  comparison 

with  fu1  and  f"J-,  (2.49)  it  is  possible  to  register  in  the  form 


t.c“  B 


AO)  _ L/U)  j _ L/-U) 

/  •  B  1  '  2B" 


(2.  57) 


moreover  in  this  case  in  it  is  found  from  the  equation 

B  —  InB  -6 

W  ft  A 


b  —  In 


2730  [Z 


A 


1544. 


(2.  58) 


Table  of  8  applications/appendices  gives  normalized  functions 
f{6)  (they  are  calibrated  to  f  (0))  for  several  values  of  B: 


£  =  4,  5,  6,  7,  8,  9,  10,  12. 


Calculations  [58]  were  carried  out  for  protons,  but  their 
applicability  to  electrons  was  not  limited  by  anything,  since  formula 
of  Dalitz  for  scattering  cross  section  was  derived  for  particles  with 
spin  1/2.  The  experimental  check  of  the  angular  distributions  of 
stray  electrons  confirms  the  theory  of  Nigam,  although  in  the  majority 
of  cases  is  a  completely  satisfactory  agreement  also  with  Moliere's 
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theory.  A  good  agreement  of  experiment  is  obtained  also  with  the 
theory  Gaudsmit-  Saunderson.  This  indicates  that  the  theories  of 
multiple  scattering  examined  at  present  can  ensure  the  reliable  data 
about  the  angular  distributions,  passing  through  the  substance  of 
electrons . 

Page  47. 

And  nevertheless  we  give  preference  in  our  calculations  to  the  theory 
of  Goudsrrit-Saunderson ,  which  does  not  require  any 

approximations/approaches  and  gives  a  "precise"  angular  distribution, 
if  a  precise  section  of  the  elementary  report/event  of  scattering  is 
known. 

Only  deficiency  in  this  theory,  being  more  of  a  psychological 
nature,  is  need  for  numerical  calculation  of  distribution.  However, 
many  problems,  which  were  being  earlier  considered  hopeless  (in  the 
sense  of  the  analytical  representation  of  the  solution),  are  solved 
rapidly  and  efficiently  with  the  help  of  the  owls  of  belt 
computational  methods. 
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Chapter  III. 

ANALYTICAL  AND  NUMERICAL  METHODS  OF  THE  STUDY  OF  THE  TRANSFER  OF 
ELECTRONS  IN  A  SUBSTANCE. 

In  this  chapter  range  of  applicability  of  analytical  method  of 
study  of  transfer  of  electrons,  based  on  solution  of  kinetic  equation, 
briefly  is  examined.  Primary  attention  is  paid  to  the  method  of 
Monte-  Carlo:  its  special  features  in  application  to  electrons  are 
explained.  Are  examined  the  basic  versions  of  the  method,  when 
tracking  is  conducted  by  a)  by  consecutive  collisions;  b)  with  the  use 
of  theories  of  multiple  scattering.  The  latter/last  version  and  its 
modifications,  used  by  different  authors,  in  detail  is  analyzed.  The 
methods  of  accelerating  the  calculation  are  considered. 

S  1.  Kinetic  equation  of  transfer. 

Kinetic  equation  for  electrons  has  the  same  structure,  as 
equation  for  7-quanta,  and  it  is  written/recorded  in  general  case  as 
follows : 

00 

~  +  uyE  f  U/**  J  dE'  J  du'  > 

K 

<  F(E',  u’,  r\  01 *(£'.  «'  -  E,  a). 


(3.  1) 
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Here  v  -  velocity  of  electron,  F(E,  u,  r,  t  )dE  -  particle  flux  at  the 

A 

moment  of  time  t  at  point  r  of  the  spaces,  energy  and  angle  of  which 
lie/rest  respectively  at  intervals  of  E-E+dE  and  u-iS+du,^!  (E,  u'-*E,  "5)  j 

-  the  probability  (per  unit  of  path)  of  the  fact  that  the  particle, 
which  has  energy  and  angle  E',  u',  will  change  them  as  a  result  of 
collision  for  E,  u  1 . 

FOOTNOTE  1 .  Here  by  u  is  understood  the  pair  of  angles  6  and  <p  in  the 
spherical  coordinates.  ENDFOOTNOTE. 

Usually  temporary/time  term  -j-  ~  is  substituted  by  the  term  9F/9R, 

t 

where  R  is  the  path  of  particle  with  this  energy  (speed)  r=-.  \  v(t  )dt'. 

o 
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Introduction  of  a  new  variable  allows  us  to  pass  by  variable  R  in 
(3.1)  to  so-called  approximation/approach  of  continuous  losses  of 
energy,  when  energy  and  path  are  connected  by  a  unique  analytical 
dependence.  Then  from  the  complex  dependence  for  H(E',  u' ,  -»E,  u)  it 
is  possible  to  isolate  the  energy  part  (it  is  expressed  as  function 
R) ,  and  only  angular  dependence  remains  under  the  integral. 

This  approximation/approach,  for  the  first  time  proposed  by  Lewis 
[59],  applicably  only  to  infinite  medium  and  in  it  are  not  considered 
fluctuations  in  energy  losses  to  ionization,  or  secondary  radiation  (6 

-  electrons,  bremsstrahlung) .  The  solution  of  the  equation  of 
transfer  are  expressions  for  the  spatial  momentsof  the  unknown 
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distribution  (spectral,  absorbed  energy,  etc.),  and  sequence  of  an 
infinite  number  of  moments/torques,  as  is  known,  it  is  equivalent  to 
the  unknown  distribution.  In  actuality  it  is  possible  to  calculate 
only  the  finite  number  (usually  not  more  than  10)  of  moments/torques, 
and,  consequently,  the  distribution  constructed  with  their  aid  will  be 
incorrect,  especially  at  the  large  depths  in  the  substance. 

Spencer  [10]  showed  that  knowledge  of  asymptotic  properties  of 
moments/torques  and  some  special  features  of  unknown  distribution  at 
large  depths  makes  it  possible  to  use  methodology,  worked  out  .for 
7-quanta  [60^  and  to  calculate  function  of  distribution  of  absorbed 
energy  *. 

j 

I 

FOOTNOTE  J.  In  the  principle  it  is  possible  to  obtain  the  spectrum  oi 
electrons  at  the  assigned  depth,  but  due  to  the  utilized 
approximations/approaches  hardly  it  is  possible  to  expect  even 
qualitative  agreement  with  the  experiment.  ENDFOOTNOTE. 

The  distribution  obtained  by  Spencer  for  the  point  and  plane  sources 
of  electrons  is  in  a  good  agreement  with  the  experiment  everywhere, 
besides  the  boundaries  of  layer.  After  the  publication  of  the  work  oi 
Spencer  it  was  very  little  done  for  further  development  of  the  methods 
of  solving  equation  (3.1).  ! 

Among  scarce  works,  dedicated  to  this  question,  should  be  noted 
works  of  Adawi  [61],  where  is  done  attempt  to  take  into  account 
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bremsstrahlung ,  and  Crew  [62],  which  solved  two-dimensional  problem 
and  were  obtained  isodose  fields  in  light  media. 

Page  52. 

The  simplified  model  of  a  deep  penetration  of  electrons  was  developed 
by  Yevdokimov  [63].  At  the  base  of  its  model  lies/rests  the  diagram, 
very  close  to  that  utilized  in  the  Monte  Carlo  method,  according  to 
which  electron  path  is  broken  line,  which  consists  of  the  rectilinear 
segments,  the  direction  of  each  of  which  is  located  by  the  averaging 
of  angular  distribution  (usually  in  the  approximation/approach,  of 
small  angles).  For  this  model  it  is  possible  to  register  the  kinetic 
equation,  whose  solution  can  be  obtained  only  at  the  large  depths. 

Thus,  at  present  there  does  not  exist  detailed  mathematical 
apparatus  for  solution  of  equation  (3.1).  Taking  into  account 
specific  difficulties  connected  with  the  role  of  the  boundaries  (for 
the  charged/loaded  particles,  which  have  small  path,  this  is 
substantial),  it  is  hardly  possible  to  expect  that  the  future  theory 
will  be  able  to  compete  with  the  Monte  Carlo  method,  which  gives  the 
sufficiently  exact  solution  of  any  problems,  although  without  the 
refinement,  inherent  in  a  precise  theory.  In  this  case  the  goal 
justifies  means. 

S  2.  Monte  Carlo  method  in  the  problems  of  the  transfer  of  electrons 


It  is  difficult  to  give  any  precision  determination  of  Monte 
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Carlo  metnod.  This  is  caused  by  the  variety  of  problems  decided  with 
its  aid.  However,  sufficiently  fully  its  essence  can  be  determined  as 
follows:  the  Monte  Carlo  method  is  the  numerical  method  of  the 
solution  with  the  help  of  the  random  tests  of  certain  circle  of  the 
physical  and  mathematical  problems,  which  carry  the  probabilistic 
(stochastic)  character,  or  by  which  it  is  possible  to  compare 
probabilistic  analog.  In  detail  the  proof  of  method  and  the  field  of 
its  applicability  is  stated  in  books  [64-66]  and  group  of  the  works, 
published  in  collectors/collect  ions  [67,  68]. 

This  above  determination  assumes  two  forms  of  problems,  solved  by 
Monte  Carlo  method.  First  of  all  this  of  the  problem  of  radiation 
transfer  and  number  of  the  problems  of  nuclear  physics,  which  are  in 
the  essence  probabilistic.  In  this  case  their  solution  is  reduced  to 
tne  selection  from  the  appropriate  random  number  distributions  so, 
that  their  mathematical  expectation  coincides  with  the  unknown  value 
a.  The  mathematical  expectation  M[$]  is  considered  as  arithmetic  mean 

>t 

n  of  the  selections  of  the  random  variables: 

/*=! 
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On  the  basis  of  the  central  limit  theorem  it  is  possible  to  claim  that 
the  distribution  arithmetic  mean  values  asimptotically  tends  for  the 
normal  distribution  and  M[£]-*a  (on  the  probability).  Error  in  this 
case  does  not  exceed  value  /3[D[  $  ]/N]  0 ,  where  /3  -  constant,  determined 
by  authenticity  level  of  evaluation/estimate,  butD[£]  -  dispersion  of 
random  variable  $  [69]. 
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M*\\l  (3.2) 

Another  circle  of  problems  does  not  have  straight/direct 
probabilistic  character  as,  for  example,  resolution  of  Laplace's 
equation  in  electromagnetic  field  theory.  For  this  equation  there 
does  not  exist  the  analytical  solution  under  the  assigned  boundary 
conditions.  On  the  other  hand,  the  same  equation  of  Laplace  describes 
the  diffusion  of  particles  in  the  region,  limited  by  the  absorbing 
barriers.  Therefore  to  Laplace's  equation  in  the  electromagnetic 
theory  it  is  possible  to  compare  as  the  analog  the  equation  for  the 
diffusion  of  particles,  which  already  can  be  solved  by  the  Monte  Carlo 
method.  Many  problems  of  gas-  and  of  hydrodynamics,  applied 
mathematics  are  solved  thus  such. 

Let  us  return  to  our  problem  of  transfer  of  electrons  in 
substance  and  will  consider  simple  problem.  Suppose  that  on  a  plate 
of  the  specific  thickness  which  is  infinite  in  two  measurements,  a 
charged  particle  beam  (electrons,  protons)  with  known  energy  spectrum 
and  the  angular  distribution  falls.  It  is  necessary  to  determine  a 
number  of  particles,  which  passed  through  the  plate,  their  energy  and 
angular  distribution. 

Interaction  of  charged/loaded  particles  with  substance  is 
complicated  stochastic  process.  Thus,  electron  can  experience  elastic 
scattering  on  the  kernel,  emit  the  quantum  of  bremsstrahlung ,  lose 
energy  as  a  result  of  inelastic  interaction  with  atom  shell,  and 
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proton,  besides  elastic  and  inelastic  scatterings,  can  initiate  the 
nuclear  reactions  (p,  n),  (p,  2n)  ...,  etc. 


Each  of  enumerated  reports/events  of  interaction  is  described  by 
specific  probability  laws,  obtained  from  theory  or  experiment.  Other 
words,  for  each  report/event  are  noted  for  only  probability 
scattei ings  to  the  specific  angle,  the  probability  of  this  energy 


loss,  the  probability  of  the  participation  of  particle  in  the  specific 
nuclear  reaction.  During  the  solution  of  our  problem  by  the  Monte 
Carlo  method  complicated  stochastic  process  is  considered  as  the 
sequence  of  elementary  actor  (simple  Markov  chain),  moreover  the 
concrete/specific  values  of  the  parameters  of  the  particles  of  such, 
as  landing  run,  energy,  direction  of  motion  after  scattering  are 
obtained  from  the  appropriate  distributions  with  the  help  of  the 
random  numbers.  The  procedure  of  the  selection  of  values  from  the 
probabilistic  distributions  with  the  help  of  the  random  numbers, 
according  to  the  conventional  terminology,  is  called  drawing. 
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Thus,  step  by  step  with  the  help  of  drawings  is  realized 
complicated  stochastic  process.  The  overall  diagram  of  calculation  by 
Monte  Carlo  in  our  case  must  include  the  following  steps:  1)  the 
drawing  of  wave  energy  and  angle  of  incidence  from  energy-angular 
source  distribution;  2)  the  drawing  of  landing  run  before  first 
interaction  from  the  range  distribution;  3)  the  drawing  of  the  form  of 
process  (elastic  or  inelastic  scattering,  bremsstrahlung ,  nuclear 


v 
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reaction,  etc.).  As  a  result  of  the  previous  drawing  let  it  be 
established  that  the  elastic  scattering  occurred,  then  from  the 
appropriate  distribution  we  generate  energy  of  particle  and  compute 

scattering  angle  (it  is  possible  to  draw  angle  and  to  compute  energy),  j 

! 

Then  procedure  is  repeated  from  point/item  2  for  the  particle,  in 

: 

which  instead  of  the  initial  values  of  energy  and  direction  of  motion  j 
are  used  the  drawn  energy  and  angle  as  long  as:  a)  energy  of  particle 
will  not  become  less  than  the  certain  minimum  value  Em,,,,  with  which 
it  is  possible  to  consider  it  absorbed  in  the  plate;  b)  particle  will 
not  leave  the  plate.  In  the  latter  case  we  simultaneously  fix/record 
energy,  and  direction  of  propagation.  But  if  in  the  process  of 
inelastic  interaction,  for  example,  nuclear  interaction,  appeared  one 
or  several  secondary  particles,  then  they  continue  from  the  collision 
to  collision  analogously  with  initial  particle.  Thus,  the  solution  of 
problem  with  the  help  of  the  Monte  Carlo  method  is  reduced  to  the 
generation  of  the  sequence  of  random  variables  E, ,  r„,  *.  rl,"-  u\'j, 
which  characterize  the  so-called  trajectory  (history)  of  particle. 
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where  E  -  energy  of  particle; 

r  -  set  of  coordinates; 

u  -  pair  of  angles  in  the  spherical  coordinates. 

Page  55. 

For  explaining  general  laws  governing  passage  of  particle  through 
plate  it  is  necessary  to  realize  large  number  of  such  trajectories  and 
to  subject  to  their  statistical  processing.  Let  us  note  that  the 
results,  obtained  during  this  processing,  contain  all  fluctuations, 
which  associate  real  process. 

By  fundamental  point  in  calculation  procedure  Monte  Carlo  is 
drawing  of  values  of  random  parameters  of  particle  r„  u,  from 
appropriate  probabilistic  distributions.  In  our  case  -  these  are 
angular  distribution  in  multiple  scattering  and  the  energy-loss 
distribution  in  the  inelastic  scattering.  Since  the  drawing  is 
repeated  repeatedly,  its  procedure  must  be  sufficient  efficient,  i.e., 
on  it  must  be  spent  little  time.  Let  us  consider  one  of  the 
frequently  used  diagrams  of  the  drawing  (other  diagrams  are  considered 
in  the  work  of  Marchuk  and  Yermakov)  (see  collector/collection  [67]). 

Let  density  of  distribution  of  certain  random  variable  x  be  given 
by  non-negative  function  p(x),  and  density  of  distribution  of  another 
random  variable  y  will  be  g(y).  Further,  let  y=f(x),  f(x)  being  the 
monotonic  function,  which  for  the  certainty  we  will  consider 
increasing.  This  means  that  for  the  arbitrarily  taken  random  number  £ 
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for  all  x<$  inequality  y<f(£)  is  fulfilled.  The  probability  of  the 
fact  that  x<$  is  accepted  to  call  the  distribution  function  and  to 
express  with  the  help  of  the  integral 

f  (*<;)=  J  p(x)dx.  (3.3) 

—  00 


It  is  obvious  that 


F[u  'fit)]-*  J'  gly)dy 


and  consequently, 


i  f  (6) 

I  p(x)dx=  j  g(y)  dy. 
—  00  — » 


(3.  4) 


(3.  5) 


Page  56. 


If  f(x)  monotonically  decreasing  function,  then  (3.5)  will  take  the 


form 


(3.  6) 


Let  the  random  variable  x  be  evenly  distributed  in  the  interval  [0 
,1],  then  p(x)=l  and  (3.5)  it  is  possible  to  rewrite  in  the  form 

t  (6) 

:  =  j  g(y)dy  =  Fiu).  (3.7) 
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In  formula  (3.7)  it  is  possible  to  express  y  through  £  ,  after 
fulfilling  inverse  transformation 

y  —  F~l  (;).  (3.8) 


Thus,  if  it  is  possible  to  fulfill  conversion  (3.8),  i.e.,  if  it 
is  possible  to  register  analytical  expression  for  F'*(£),  we  obtain 
convenient  method  of  expression  of  random  variable  from  arbitrary 
distribution  law  through  another  random  variable  from  uniform 
distribution  [0,  1]  1. 


FOOTNOTE  1 .  The  numbers,  comprised  of  numerals  0,  1,  2,  3,  4,  5,  6, 

7,  8,  9  so,  that  in  their  location  there  is  no  law  whatever,  are 
called  random.  In  order  to  obtain  random  number,  for  example  4-digit 
in  interval  [0,  1],  is  necessary  the  number,  comprised  of  four 
numerals,  to  multiply  by  10' 4 .  ENDFOOTNOTE. 

If  for  the  solution  of  problem  not  many  random  numbers  are  required, 
then  it  is  possible  to  take  them  from  appropriate  tables  [70].  In 
practice  they  use  the  so-called  pseudo-random  numbers.  These  numbers 
in  a  quantity  from  several  hundred  thousand  to  several  million 
(depending  on  program)  are  generated  by  computer  itself  in  the  process 
of  calculation  according  to  specific  recursion  relations  [71].  Very 
fact  of  obtaining  the  sequence  of  such  numbers  according  to  a  defined 
rule  as  if  contradicts  name  itself  "random  numbers".  However,  such 
sequences  are  used  only  then,  when  they  satisfy  to  all  statistical  to 
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criteria  on  checking  of  "chance"  and  uniformity,  and  therefore 
virtually  they  do  not  introduce  systematic  errors  into  the 
calculation.  We  used  the  program  of  the  generation  of  pseudo-random 
numbers  for  computers  BESM,  given  in  monograph  [72]. 

Page  57. 

As  an  example  of  the  use  of  a  method  of  inverse  function  let  us 
consider  the  diagram  of  the  generation  of  a  random  value  y  from  the 
density  of  the  distribution 

g(u)  =  Ay~a,  (3.9) 

where  the  variable  y  is  changed  in  interval  of  Y°-Y~Yi.>  and  constant  A 
is  determined  from  the  condition  for  the  standardization 

!/i 

j  g(i')dy  =  1. 


In  the  case 

v 

*  =  A  \  y~"  dy  =  "+1  -  y~n+l), 

Vo 


in  question  whence 


"  yon+i-i~-  ’  <3- 10) 

thus,  the  drawing  of  random  variable  from  distribution  (3.9)  is 
reduced  to  the  generation  of  random  number  £  in  interval  [0,  1]  and  to 
the  calculation  of  fraction  (3.10). 


m 
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Unfortunately,  by  no  means  in  all  cases  is  it  possible  to  fulfill 


inverse  transformation  (3.8).  Thus,  it  is  impossible  to  register  in 


the  analytical  form  inverse  function  for  the  function  of  the 


distribution  of  Moliere,  Goudsmit-Saunderson ,  Landau  -  fundamental 


distributions,  utilized  in  the  problems  in  the  passage  of  electrons, 


For  these  functions  it  is  possible  to  find  only  numerical  values 
y 


F  (y)  =  \  since  the  method  of  inverse  function  is  comparatively 


simple,  then  it  it  is  possible  to  use,  also,  in  such  a  case,  when 


function  F(y)  is  assigned  numerically. 


Let  function  F(y)  be  assigned  for  sufficiently  large  number  of 


equidistant  values  which  contain  entire  range  of  change  y  (values 


of  functions  of  Goudsmit-Saunderson  we  they  computed  for  80  points  in 


interval  0-2jt).  Then  by  interpolation  it  is  easy  to  find  such  set 


which  corresponds  to  the  sequence  of  the  evenly  distributed  values 


of  function  F(y) ,  determined  with  the  help  of  the  relationship/ratio 


Fiyj 


(3.  11) 


where  m=l ,  2 ,  3  .  . . .  N. 


Page  58. 


These  values  i/„,  will  be  brought  in  into  memory  of  machine.  The 


drawing  of  value  y  is  reduced  now  to  the  generation  of  random  number  £ 


and  the  determination  of  the  number  of  the  cell  m(m=£N),  in  which  the 


is  written  the  corresponding  value  of  y.  If  m=|N  is  not  whole  number, 


then  value  y  is  located  by  interpolation  between  the  values,  carried 


A  .  *.  ."N 
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into  the  next  arranged/located  cell.  The  behavior  of  both  the  integral 
function  of  Goudsmit-Saunderson  and  of  the  function  of  the 
distribution  of  angles  as  a  result  of  single  scattering  is  such,  that 
they  rapidly  approach  one  for  the  angles  of  scattering  close  to  180°. 
In  this  case  it  would  prove  to  be  that  the  equiprobable  values, 
located  in  two  adjacent  nuclei,  could  correspond  to  two  essentially 
distinguished  angles,  and  the  accuracy,  obtained  by  linear 
interpolation,  would  be  low.  Therefore  the  described  method  it  would 
be  possible  to  use  twice:  to  find  the  values  of  angles  for  the  random 
numbers  from  0  to  0.99  and  additionally  100  values  of  angles  f.or  the 
random  numbers  from  0.9900  to  1.0000.  If  the  memory  of  machine  does 
not  make  it  possible  to  retain  this  further  information,  it  is 
possible  to  use  this  procedure.  In  the  memory  of  machine  it  is  brought 
in  the  sequence  of  unevenly  distributed  random  numbers  F('/,\, 
corresponding  N  to  equidistant  values  (for  example  angles),  and  we 
carry  out  between  them  linear  interpolation,  when  F[yt)-  ^  F(yix).  In  - 
this  case  the  accuracy  in  the  evaluation/estimate  of  random  value  . 
will  be  worse  y-. 

S  3.  Networks  of  Monte  Carlo. 

If  we  knew  precise  laws,  which  characterize  reports/events  of 
interaction  of  electrons  with  substance,  then  result  of  solution  of 
problems  by  method  by  Monte  Carlo  method  would  depend  virtually  only 
on  statistical  accuracy  of  calculation.  However,-  the  interaction 
cross  sections  of  electron  with  the  atoms  are  known  approximately  '-n 
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the  theoretical  approaches  they  are  used  its 

approximations/approaches,  in  the  experimental  data  spread  it  is 
determined  by  measuring  errors).  Furthermore,  usually  is  introduced 
an  also  whole  series  of  simplifications.  In  particular,  it  is  assumed 
that  the  electron  simultaneously  interacts  only  with  one  atom.  Here 
the  effects  connected  with  the  interaction  of  the  electron  with 
several  atoms  are  not  considered.  This  simplification  makes  it 
necessary  to  be  bounded  to  the  examination  only  of  isotropic 
materials. 


Page  59. 


At  the  same  time  during  the  research  of  the  transfer  of  narrow 
electron  beams  in  the  thin  single  crystals  the  simplification 
indicated  cannot  be  accepted,  since  the  phenomenon  of  channelization 
discovered  recently  indicates  the  significant  role  of  the  anisotropic 
properties  of  medium.  With  low  energy  of  the  electrons  (-61  keV)  it  is~ 
necessary  to  also  consider  interaction  with  the  electronic  plasma, 
which  leads  to  the  appearance  of  discrete/digital  energy  losses. 


If  we  assume  that  in  problems  examined/considered  by  Us>  all 
simplifications  indicated  above  occur,  then  in  this  case  attempt  to 
realize  at  trajectory  with  method  of  consecutive  collisions  leads  to 
another  very  essential  difficulty  -  to  high  expenditures  of  machine 
time.  Actually,  due  to  the  high  value  of  the  interaction  cross 
section  of  electron  concerning  the  atoms  for  following  one  trajectory 
it  would  be  necessary  to  repeat  the  procedure  of  the  drawings  of 
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<U 

10s-10‘  times.  After  taking  the  time,  necessary  for  calculating  set 
of  parameters  of  particle  (E,  r,  y)  after  one  interaction,  equal  to 
100  us,  we  will  obtain  that  for  calculating  1000  trajectories 
(necessary  minimum  for  guaranteeing  the  statistical  authenticity  of 
result)  it  will  be  required  by  10 *-  10 5  s,  and  this  already  too  long  a 
time.  Let  us  note  that  in  the  problems,  connected  with  the  transfer 
of  7-quanta,  and  in  the  majority  of  neutron  problems,  it  is  necessary 
to  trace  only  20-30  collisions.  For  a  number  of  reasons  for  the 
electrons  it  is  not  possible  to  use  known  procedures  of  the 
acceleration  of  calculation. 

I^Jfiagram  of  the  "amalgamated"  collisions  l. 

FOOTNOTE  In  the  literature  sometimes  is  used  the  term  -  schematic 
of  +he  "condensed"  trajectories.  ENDFOOTNOTE. 

A  way  out  was  found  in  application  of  model  of  "amalgamated" 
collisions,  i.e.,  when  drawings  are  produced  not  after  each  collision 
from  distributions,  which  characterize  elementary  report/event  of 
scattering,  but  after  passage  of  certain  route  segment  At  from 
distributions,  which  characterize  multiple  scattering  of  particles. 
Consequently,  instead  of  the  "real"  particle  trajectory,  schematically 
depicted  in^igure  13  (fine/small  zigzag),  is  examined  hypothetical 
trajectory  (large/coarse  zigzag).  It  is  difficult  a  priori  to 
predict,  how  this  a  replacement  will  affect  the  results  of 
calculation.  For  the  proof  of  this  approach  it  is  possible  to  use  the 
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following  considerations.  Since  elastic-scattering  cross  section  for 
small  angles  is  Rutherford,  the  probability  of  divergence  of  such 
angles  will  be  large  and,  consequently,  angle  in  multiple  scattering  - 
this  the  result  of  numerous  single  small-angle  scatterings. 

Page  60. 

Under  these  conditions  it  is  possible  to  isolate  a  certain  mean 
direction,  along  which  moves  the  electron,  and  thus  to  fix/record 
segment  Af,.  It  is  clear  that  the  less  the  more  precisely  is 
fixed/recorded  the  mean  direction  of  electron  motion  and  the  more 
precisely  it  is  possible  to  examine  its  motion  near  the  boundaries, 
but  in  this  case  the  time,  necessary  for  conducting  the  calculations, 
substantially  grows. 

Furthermore,  value  Af,  cannot  be  done  less  than  certain  minimum 
value,  below  which  become  inapplicable  theories  of  multiple 
scattering.  Thus,  for  example,  theory  of  Moliere  is  applicable,  when 
collision  frequency  in  segment  Ar,  is  more  than  20,  i.e.,  for  the 
thicknesses  is  more  than  several  mg/cm2.  On  the  other  hand,  value 
lci  cannot  be  done  very  large,  so  as  to  energy  losses  to  the 
ionization  and  the  bremsstrahlung  would  be  compared  with  the  wave 
energy.  Let  us  recall  that  energy  losses  not  at  all  are  considered  in 
the  theories  of  multiple  scattering.  Thus,  the  selection  of  length 
to  a  certain  extent  is  arbitrary,  and  with  this  fact,  probably, 
can  be  connected  some  small,  but  systematic  divergences,  which  are 
exhibited  during  the  comparison  of  calculations  with  the  experiment. 
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Fig.  13.  Diagrammatic  representation  of  electron  paths  in  substance 
Key:  (1).  reflection.  (2).  absorption.  (3).  transmission. 


Nevertheless  is  possible  to  select  so  as  to  reduce  to  a  minimum 
these  divergences. 


In  work  of  Berger  [68]  different  methods  of  selection  ±tt.  are 
investigated • 


1.  Assume  that  is  given  energy  vt  at  the  end  of  the  i  segment; 
let  us  require  so  that  relative  change  in  energy  of  particle  in 
segment  tt  would  compose  always  one  and  the  same  number 


k(k<l ) 


El~  J  (~Jr}dt 

_ it 


(3.  12) 
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where  -  (dE  )/  (-BT)  -  average/mean  energy  losses,  computed  from 
formula  (2.8).  Assigning  value  of  k,  it  is  possible  to  compute  that 
corresponding  At^'* 


With  this  method  of  determination  At,  average/mean  cosine  of 
angle  of  deflection  proves  to  be  in  effect  constant  for  entire 
sequence  of  segments  Usually  k  choose  in  the  form  '•  -=2_,,m,  where  m 
-  whole  number  (in  work  [68]  m=32).  This  is  convenient  because  near 
the  boundary  of  layer  for  a  more  precise  evaluation/estimate  of  energy 

I 

loss  with  the  help  of  interpolation  it  would  be  possible  to  decrease  k 

.  .  -  /  fan 

m  n  of  times,  i.e.  to  examine  segments  At,',  determined  for  oL 
Naturally,  with  this  approach  the  consecutive  account  of  fluctuations 
in  energy  losses  is  impossible  and  they  have  in  mind  only  continuous 
degradation  of  energy  by  small  portions  =*£,. 


2.  It  is  possible  to  assume  that  At,  const,  as  this  is  done  in  one- 
of  first  works  with  use  of  method  of  Monte  Carlo  [74],  where  passage 
of  electrons  with  energy  in  5-55  MeV  region  was  examined  through  layer 
of  carbon.  Tracking  ceased,  when  energy  of  electrons  reached  1  MeV. 


Gradual  broadening  of  angular  distribution  toward  the  end  of 
landing  run  is  deficiency  in  this  method,  which  can  serve  as  source  of 
systematic  error. 


3.  It  is  possible,  finally,  to  assume,  as  we  did,  following  work 
of  Perkins  [75]  tht  At,  varies  with  energy  in  process  of  particle 


M.VY'  V-V- .N 
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motion 


Afi+1  =  kEt.  (3.  13) 

It  is  easy  to  see  that  if  -dE/dt«  const,  and  this  is  implemented 
virtually  always,  when  segment  At ,  is  sufficiently  small  so  that 
energy  losses  on  it  would  be  small  (AE<0.1Eo),  formula  (3.12)  is 
reduced  to  (3.13) . 

Page  62. 

Actually 

Et  —  cAf,.H  = 

where  k*  -  this  coefficient  k  in  formula  (3.12),  and  -dE/dt=c.  In 
order  to  use  formula  (3.13)  for  any  energies  of  electrons,  we  consider 
coefficients  of  k  as  energy-  dependent. 

Table  8  gives  values  k=k(E),  which  are  used  in  our  calculations. 
This  method  of  selection  At,  makes  it  possible  to  consider 
fluctuations  in  energy  losses  and  to  examine,  consequently,  is  more 
fully  the  transfer  of  electrons  in  the  substance. 

According  to  Berger’s  diagram  1  after  determination  u,  ,  goes 
drawing  of  scattering  angle  from  distributions  of  Moliere  or 
Goudsmit-Saunderson,  moreover  input  energy  for  obtaining  this 
distribution  is  energy  E,,  which  electron  in  the  beginning  of  segment 
had. 

Ei  i-~  0,1,2,...  (3.14) 


DOC  =  86063205 


PAGE 


FOOTNOTE  1 .  Apparently,  Berger  uses  a  diagram  taking  into  account 
fluctuations,  but  it  is  nowhere  published.  ENDFOOTNOTE. 


*;S 


Transit  ion/ junction  into  the  necessary  reference  system  is  produced 
with  the  help  of  the  formulas  of  spherical  trigonometry.  Azimuth 
angle  is  assumed  to  be  that  distributed  evenly  in  the  interval  0-2ff, 
and  its  random  value  is  realized  with  the  help  of  the  further  drawing, 
Cartesian  (or  other)  coordinates  of  the  salient  points  of  electron 
path  simultaneously  are  calculated.  The  sequence  of  the  stages  of 
calculation  according  to  our  diagram  is  such:  1)  for  this  energy  of 
electron  we  choose  from  Table  8;  2)  with  the  help  of  coefficient 

of  k  we  find  ■**»  —  *£ ti  3)  we  compute  A,  and  for  current  energy  £» 
according  to  formula  (2.19). 


WE*. 


Key:  (1).  MeV. 

Page  63. 

Methodology  of  drawings  from  integral  function  of  distribution  of 
Simon  (Landau)  for  energy  losses  only  to  ionization  (see  Table  6  of 
application/appendix)  is  analogous  this  described  in  S  2  of  the 
chapter.  Let  us  note  that  the  distribution  function  was  approximated 

0.8463 _ 

by  curve  l*E  for  all  —  !LZi£»  io.  After  finding  £*  ,  we 

V  •io  / 

compute  the  loss  of  energy  A£u.  For  all  Et >0,511  MeV  through  formula 
(2.11)  we  find  energy  losses  to  the  bremsstrahlung  in  segment  Af,.  We 
compute  the  total  losses  of  energy  *£„,  =*■*£*,+  s srl  and  we  find  energy 
of  electron  £,-n  at  the  end  of  segment 

Ei-i-i  —  Ei  &Enl’ 

Upon  consideration  of  fluctuations  in  energy  losses  to 
bremsstrahlung  drawing  of  total  losses  of  energy  is  produced  from  the 
tabulated  distributions  and  those  input  into  memory  of  computer  (2.26) 
for  entire  sequence  of  energies  and  the  connected  with  them  stL.  In 
this  case,  for  current  energy  Ei  is  chosen  the  distribution,  designed 
for  the  value  of  the  energies  nearest  to  £,  Reaching/achievement  of 
certain  minimum  value  £■»»,  which  is  assigned  by  researcher  and  is 
fixed/recorded  in  the  memory  of  machine,  is  the  first  criterion  of  the 
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cessation/discontinuation  of  the  tracking  of  this  trajectory.  In  our 
calculations  £,uln  ^0,02^0,05  MeV.  Electron  energy  £  £“‘°  is  considerd 
to  be  absorbed  directly  in  the  vicinity  of  the  point,  where  this 
condition  is  satisfied,  and  the  tracking  of  new  trajectory  begins,  if 
j<N,  where  N  -  number  of  planned  for  the  tracking  trajectories.  If 
E.  £ fiiiai  we  pass  then  to  the  drawing  of  scattering  angle.  The  logic 
circuit  of  the  calculations  of  energy  of  electrons  is  shown  in  Figure 
14. 


4.  Schematic  of  drawing  of  scattering  angle  0,.. ,,  calculated  off 
direction  electron  motion  in  segment  ±t,,  it  is  shown  in  Figure  15.  In 
this  case  on  input  energy  E,  is  found  out  the  nearest  energy,  for 
which  the  integral  distribution  of  Moliere  or  Goudsmit-Saunderson  is 
known,  and  very  methodology  of  drawing  is  analogous  described  above 
for  the  distributions,  assigned  numerically.  For  the 
transition/ junction  into  the  reference  system,  connected  with  some 
specific  direction,  for  example,  by  the  positive  axis  02,  are  used  the 
formulas  of  spherical  trigonometry 

C08<|»|  =  cos  <{•(_!  cos  6,  -f  sin  <{*.-1  sin  cos  yA.  (3.  15) 

Angles  0,,  f  =  1,  2,...ix  are  shown  in^igiire  16.  Azimuth  angle  '/■< 

is  distributed  equiprobably  in  the  interval  0-2ir,  if  we 
disregard/neglect  polar izational  effects.  The  schematic  of  the 
drawing  of  this  angle  is  trivial:  random  number  £  is  chosen  and  it  is 
computed  x=2 
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Fig.  15.  Logic  circuit  of  calculated  angle  in  multiple  scattering. 
Key:  (1).  Random  number  £.  (2).  number  of  angular  intervals.  (3). 

no.  (4).  yes.  (5).  determination  of  angle  v  from  distribution  by 
interpolation.  (6).  Continuation  of  tracking. 

Page  65. 

Calculation  cos  *  we  produced  not  on  standard  program  for  this 
computer  1 ,  but  we  used  the  table  of  sin  y  (y  =  0° 90°),  the  registered 
into  the  memory  machine  with  a  space  of  /  . 

FOOTNOTE  1 .  There  is  an  efficient  diagram  of  Monte  Carlo  [76],  which 
uses  the  methodology  of  exception/elimination,  which  makes  it  possible 
to  obtain  immediately  cos  *  and  sin  Xf  using  on  the  average  2,  6 
random  numbers.  ENDFOOTNOTE. 
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Turning  to  this  table  occurs  in  the  dependence  on  £  through  the 
functions 

sin  (^2 - 2«j  ,  i  <  0,25, 

-  sin  |  2*5  -  ,  0,25  <  5  <  0,5, 

sin  ( ~  -  2«).  0,5  '0,75,  (3.  16) 

sin (2*5-  -yj,  0.75<5^1. 


5.  Cartesian  coordinates  of  consecutive  salient  points  of 
trajectory  are  computed  from  relationships/ratios 

i 

*1  =“  V  At,  coe  (3.17) 

1 

( 

x,  —  V  At,  sin  ‘J'i-i  cos  <ft, 

(3.  18) 

V  At,  sin  <J*i— 1  sin  <p, 

(xo=yo=zo=0  -  coordinate  of  impact  point  in  electron).  Angles 
included  in  (3.18),  are  located  with  the  help  of  the  formulas 


sin  (?,-ri  —  <Pj)  — 
COS(f,  !“?()  = 


Bin  (<i»in  Xi 
sin  ’ 

COB  U,-  COS  'Vi  COS  i,  ,  [ 


(3.  19) 


sin  '(-i-sln  'f,,t  , 


It  is  assumed  that  <p„=0.  Since  the  majority  of  problems  is  solved  in 
the  one-dimensional  case,  only  coordinate  z,  computed  from  (3.17) 
interests  us. 
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Page  66. 

For  spherical  geometry  (Fig.  17)  tracking  electron  occurs  not  on 
Cartesian  coordinates  of  end/lead  of  segments  nt,  but  on 
radius-vector  Ii,,  carried  out  into  end/lead  of  this  segment.  In  this 
case  the  reference  system  of  angles  is  conveniently  connected  with 
radius-vector  J*,.  It  is  not  difficult  to  establish/install  simple 
correlations  for  determination  r,  and  ^ 

=  it,-  ;  Arrl—  2Al,,  1  /?,  C0S<5»,,  (3.  20 

cos  +1  n  =  cos  Cf|-H  cos  a,+i -j-  sin  bj+i  sin  *,  tcos/^i.  ^.3.  21 

Values  cos  a, n  and  sin«.-ii  in  (3.21)  can  be  found  from  the 
expressions: 


*'+i  =  to  -  (3. 22 

ain'^ 
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ir  :»r'»r;v  .  v~  '  w  v*  '  •-  V'  ’V 


v  ' \r ;  v*”v- -v  rwr  'V'  * 


boundary  occurred.  Characteristics  (energy,  angle),  with  which  the 
electron  intersects  will  be  input  into  the  storage  cells  of 

computer  corresponding  to  this  layer,  tracking  is  continued,  until 
electron  is  absorbed  in  substance  (A’  ii  it  will  not  cross 

boundary  z\ vp  of  sample,  which  separates/liberates  it  from  the 

vacuum.  It  is  necessary  to  now  consider  the  method  of  calculating  the 
characteristics  of  electron  outgoing  from  the  boundary.  Distance  of 
boundary  we  find  from  the  formula 


-I  Afi 


rW 


(3.  24) 


where  2. .  1.  and  -  coordinate  of  the  ends/leads  of  segment  \f » ,  1 . 
Strictly  speaking,  we  must  play  values  E\  ;  1*  at  the  end  of  the 
segment  from  the  appropriate  distributions  as  this  done  above. 

However,  if  is  easy  to  compute  on  the  basis  of  the  universal 

distribution  of  Simon,  then  for  drawing  'v,,  it  was  necessary  to  have 
a  distribution  of  angles  in  multiple  scattering,  designed  for  &c£- 
This  would  lead  to  the  further  charging  of  the  memory  of  machine.  For- 
the  approximate  determination  of  angles  and  energies  linear 
interpolation  was  used.  Since  the  checking  to  the  intersection  of 
boundaries  is  produced  after  that  how  is  done  next  space  At.n.  we 
know  set  of  parameters  E,,  ti.  zi  in  the  beginning  of  segment  and 
Ent,  ♦h-w  *»'•■  at  the  end  of  it,  then 

ehi  ^E‘  +  ~st77T  (£l!  l  ~e‘> 

and  (3. 25) 

Page  68. 
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In  contrast  to  this  methodology  3erger  [68]  calculates  energy 
E, , ,  from  average/mean  losses  in  segment  it*'.  Undoubtedly,  both 
approaches  lead  to  the  errors  in  the  angular  and  energy  distributions 
of  the  outgoing  electrons,  but  they,  apparently,  they  are  small,  since 
the  results  of  calculations  will  agree  well  with  the  experiment.  One 
should  note  that  tracking  after  the  intersection  of  boundaries 

(2«P  ‘  z‘.v  )  is  continued  with  sets  Et  1,  ti-n.  <-h*  recorded 

in  the  computer  memory,  which  are  obtained  correctly. 

II.  Diagram  of  of  "Successive"  Collisions. 

Fundamental  results  of  applying  Monte  Carlo  method  to  transfer  of 
electrons  are  obtained  according  to  diagram  presented  above  of 
"amalgamated"  collisions.  Nevertheless  several  works  were  carried  out 
according  to  the  diagram  of  "consecutive"  collisions.  The  target  of 
such  calculations,  which  are  implemented  only  for  the  very  thin 
samples,  to  investigate  the  possible  systematic  errors,  introduced  by 
the  diagram  of  the  "amalgamated"  collisions,  and  to  explain  the  effect 
of  interaction  cross  section  on  the  characteristics,  which  especially 
sensitive  to  its  form,  for  example,  scattering  to  the  large  angles  or 
to  the  energy  losses  in  the  layer  of  substance.  This  program  is 
formulated  in  the  works  of  MacCallum  [77]  and  Paul  and  Tatzber  [78]. 
Somewhat  later,  Hara  [79]  fulfilled  the  analogous  calculations  even 
for  a  specific  problem  -  the  evaluation/estimate  of  heat  release  in 
the  window,  through  which  the  electron  beam  is  withdrawn  from  the 
accelerator.  The  logic  circuit  of  the  tracking  of  trajectory, 


S 
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borrowed  from  work  [793/  it  is  shown  in  Figure  18.  The  drawing  of  the 
landing  run  between  collisions  is  produced  from  the  distribution 


P(x,  E)dx  -  :‘,.oxp(  -  \K.x)dx, 


(3.  26! 


where  coefficient  is  connected  with  total  cross  section  :,,{£)  with 
the  relationship/ratio  i>7) 


st  v.(f» 


In  [80]  the  ef f ective procedure  of  drawing  from  this  distribution  is 
described. 

Page  69. 

Depending  on  energy  range  in  question  total  cross  section  is 
composed  of  sections  of  processes  of  1)  elastic  scattering;  2) 
inelastic  scattering;  3)  bremsstrahlung ,  etc.  In  the  following 
drawing  the  form  of  process  from  the  distribution 


Pi  —  .  pi  “  >  •  • -  Pa  “  1  - -  Pi  -  ...  Pr.  I.  (3.  28) 


is  clarified- 


For  this  random  number  $  is  equal  consecutively/serially  with 
relations  (3.28),  and  if,  for  example,  Px<$£pa,  then  we  have  process  2 
and,  etc.  Expressions  for  plr  p2t  jU^  are  obtained  by  the 

-A 

integration  of  the  corresponding  expressions  for  the  differential 
cross  sections.  During  the  calculation  of  (elastic  collisions)  in 
work  [793  Rutherford  section  is  undertaken. 
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Fig.  18.  Logic  circuit  of  tracking  of  electron  by  "consecutive" 
collisions . 

Key:  (1).  input  data  .  .f^F\ 2).  random  number.  (3).  determination 

of  path  ....  (4).  determination  of  energy  loss  to  excitation.  (5). 

Scattering  on  kernel.  (6).  random  number.  (7).  Scattering  on 
electron.  (8).  random  number.  (9).  no.  (10).  yes.  (11). 
determination  of  scattering  angle.  (12).  determination  of  azimuth 
angle.  (13).  determination  of  energy  of  stray  electron:.  (14). 
calculation  of  energy  loss.  (15).  end/lead.  (16).  determination  of 
coordinate  ....  (17).  ...  (backward)  ...  (inside):.  (forward). 

(18).  backward.  (18a).  forward.  (19).  determination  of  energy  on 


DOC  =  86063206 


PAGE 


/v* 


boundary.  (20).  Inside.  (21).  Printing.  (22).  total  losses  of 
energy. 

Page  70. 

Obviously,  there  are  no  fundamental  difficulties  whatever  for  the 
calculation  u x  ,  if  we  use  the  tabulated  values  of  the  section  of  Mott. 
Inelastic  electron-electron  scattering  was  described  by  Moller's 
formula  (1.20).  For  this  process  m»  it  is  computed  from  formula 
(1.25) . 


After  explaining  form  of  process,  we  find  energy  or  angle  of 
stray  electron  from  integral  distributions 


/  oni4t 


<><">«  y  #*/  .f 


^  I 

~  \  Te  dE!  \  (3‘ 


I  ruin 


where  &  £mi0>  0mix,  Em,x  ~  minimum  and  maximum  angles  and  energy  of 

electrons  after  scattering. 


The  manifestation  of  the  minimum  angles  is  caused  by  divergence 
of  section  at  9-0°  (for  example,  Rutherford  without  shadowing),  and 
minimum  energies  -  by  process  of  exciting  electron  shell.  The  account 
of  energy  loss  to  the  excitation  is  produced  on  the  average/mean 
losses,  using  formula  (2.3).  For  it  is  possible  to  use  the 

expressions  for  the  differential  cross  sections,  given  in  chapter  I. 
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Since  energy  and  angle  are  connected,  then,  after  playing  of 
distributions  (3.29)  or  (3.30)  one  of  these  parameters,  it  is  easy  to 
find  another.  In  the  case  of  elastic  scattering  usually  gaming 
process  is  conducted  from  (3.29),  and  they  disregard  energy  losses 
(scattering  on  the  heavy  nucleus). 

After  computing  scattering  angle,  we  give  it  into  necessary 
reference  system  according  to  formula  (3.15)  and  determine  coordinates 
of  particle.  As  can  be  seen  from  diagram  in  figure  18,  in  it  are  used 
the  same  criteria  of  the  break  of  trajectory  as  in  the  diagram. of  the 
"amalgamated"  collisions.  The  results  of  the  calculations,  carried 
out  by  the  authors  of  the  cited  works,  are  sufficiently  plain.  In  work 
[78]  with  accuracy  of  ±20%  the  agreement  with  the  data  of  Berger  [68] 
for  the  back-scatter  factors  is  obtained.  A  good  agreement  of  these 
calculations  with  the  experiment  according  to  the  coefficient  of 
transmission  Ts  is  noted  also.  More  interesting  information  is 
obtained  in  work  [79].  The  calculated  values  of  absorbed  energy  are 
in  a  good  agreement  with  experiment  [81].  In  Table  9  we  compare  these 
data  with  the  data  of  our  calculations  for  A1 . 

Page  71. 

As  can  be  seen  from  Table  9,  there  is  good  agreement  between 
calculations.  This  can  be  considered  as  confirmation  of  the 
applicability  of  the  model  of  the  amalgamated  collisions  generally  and 
even  in  the  region  of  such  thicknesses  1 ,  where  the  criteria  of  the 
applicability  of  the  theory  of  multiple  scattering  strictly  are  not 
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implemented  (see  latter/last  column  •Table  9). 

U 

FOOTNOTE  1 .  However,  in  this  region  the  calculation  procedure  Jara 
[79]  should  be  considered  as  the  version  of  the  diagram  "catastrophic" 
collisions.  ENDFOOTNOTE. 

Belyayev  and  Krupman  [82]  carried  out  calculation  of  angular 
distribution  of  electrons,  scattered  in  very  thin  layer  of  substance. 
The  drawing  of  landing  runs  in  this  case  was  produced  from  the  Poisson 
distribution.  A  good  agreement  with  the  experiment  is  obtained. 

In  considered  diagram,  the  use  of  different  methods  of 
accelerating  calculation  would  raise  speed  and  efficiency  of 
calculation.  For  example,  the  correlated  tests  (see  below)  would  make 
it  possible  to  consider  the  validity  of  one  or  the  other  methods  of 
the  account  of  the  shadowing  and  other  theoretical 

prerequisites/premises.  Conducting  such  calculations  is  of  undoubted - 
interest . 

III.  Diagram  "catastrophic"  collisions. 

In  this  diagram,  proposed  by  Schneider  and  Cormack  [83],  such 
rare  events,  as  scattering  to  large  angle  o  ,  0„  or  large  relative  loss 
of  energy  w  wm,  are  called  "catastrophic"  and  are  examined  separately. 
Particle  trajectory  consists  of  the  segments,  within  which  there  are 
no  catastrophic  collisions,  and  the  tracking  of  particle  is  produced 
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according  to  the  usual  diagrams:  by  Berger,  ours,  etc.  The  most 
"catastrophic"  collisions  occur  only  at  the  ends/leads  of  the 
segments . 
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Table  9. 


£. 

M  >a 

1 

|  o)  AX,  Mi* 

£ ) 

Hmcjio  co 
yAapettMH 

■  100-MK.U 

100  MKM 

j  200  | 

300 

[79] 

|  Ml h 

pacneTu 

[79] 

PACH6TU 

[TO] 

0,4 

0,063 

0,072 

0,165 

0,186 

0,299 

0,288 

8,07 

0,7 

0,052 

0,048 

0,112 

0,122 

0,192 

0,224 

2,77 

1,0 

0,062 

0.019 

0,102 

0,100 

0,168 

0,170 

1.56 

Key:  (1).  MeV.  (2).  m-  (3).  Collision  frequency  in  100  ^m. 

(4).  our  calculations. 

Page  72. 

The  sequence  of  the  states  of  particle  is  written/recorded  with  the 
help  of  the  matrix/die 

£0;  Eu  £,*;  E,,  £,*;••• 

rui  rl»  r  1*.  r'l  •  •  •  •  r<»  r<  •  •  •  • 

U-oi  till  ^1  5  ti2i  ii*2  •  •  •  *  tij ,  tij  *  •  •  • 

< 

The  set  of  parameters  Eit  rit  u,  and  £,\  f.\  u,*  character  izes  the  state  of 
particle  respectively  to  the  i~th  and  immediately  after  the  i-th 
"catastrophic"  collision. 


Diagram  of  drawings  is  such: 


1.  Value  of  Q  is  determined  from  the  distribution 

P(Q)dQ  -  e  **dQ.  (3.  31) 


2.  Knowing  £,*,  we  find  A’, -i  (energy  to  i  +  1  catastrophic 
collision)  from  formula 
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where  %  -  landing  run  with  respect  to  loss  of  energy  w  u-...,,  and 

I  <i£\ 

r  dTj  -  average/mean  energy  losses,  computed  from  formula  (2.7).  In 
the  case  of  electron-electron  scattering  -•  .  1  ,  where  ;<«,j  is 

obtained  from  (1.25)  by  substitution  instead  of  w’ . 

3.  Knowledge  and  A',‘  makes  it  possible  to  find  value  of 

segment  between  i  and  i+1  "catastrophic"  collisions. 


d£; 
A  m  \ 


dx  i 


,a.  32) 


4,  Rate  of  energy  loss  in  "catastrophic" 
the  help  of  drawing  from  distribution 

w  j: 

P(U.  LVm)  =  \  j”  dW  j  \  Ju;  u'“'' 
Um  I 

where  u-  is  determined  by  formula  (1.20). 

(tu> 


collision 


(3. 


is 


found  with 


Using  played  value  of  w,  we  find  energy  of  stray  electron  t’,* 
and  energy  f  5=elektrona  £,*=(1  —  u>)£* 

E„  =  wEt. 

With  the  help  of  (1.21)  and  (1.22)  it  is  easy  to  find  new 
directions  of  motion  of  both  electrons.  Transition/ junction  into  the 
necessary  reference  system  is  completed  employing  usual  procedure. 

Page  73. 

In  the  extreme  case,  when  a',  ^  y,  diagram  "catastrophic" 
collisions  passes  into  diagram,  where  are  considered  only  continuous 
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losses  of  energy  (/.m  =  _).  In  other  limiting  case,  when  w„,  o  (mean  free 
path  between  the  collisions  it  becomes  equal  to  0),  all  collisions  are 
"catastrophic",  i.e.,  this  diagram  passes  into  the  diagram  of 
consecutive  collisions. 

Thus,  diagram  "catastrophic"  collisions  occupies  intermediate 
position  between  diagrams  examined  above.  This  makes  it  possible  to 
hope  that  during  the  proper  development  it  will  be  possible  to 
optimally  combine  advantages  of  both  maximum  diagrams  in  it.  Even  now 
are  clearly  visible  certain  the  valuable  sides  of  the  diagram:,  the 
more  correct  account  of  fluctuations  in  the  losses  of  energy  (drawing 
it  is  conducted  only  for  the  large  losses)  and  obtaining  the  detailed 
information  about  the  energy  and  angular  distribution  of  6-electrons. 
In  the  future  its  application  for  the  study  of  backscattering  is 
possible. 

Recently  diagram  "catastrophic"  collisions  was  used  by 
Kolchuzhkin  and  Shevtsova  [84],  They  showed  that  the  results  of 
calculations  according  to  this  diagram  will  agree  well  with  the 
experiment,  and  they  investigated  the  methods  of  evaluation/estimate 
and  En,.  To  the  diagram  "catastrophic"  collisions  examined  closely 
borders  on  Meissner's  diagram  [85],  in  which  is  done  the  first  attempt 
to  combine  in  one  calculation  the  analytical  methods,  based  on  the 
solution  of  the  equation  of  transfer,  and  the  Monte  Carlo  method. 

This  combined  methodology  has  already  been  used  in  transport  problem 
of  7-radiation  [86]. 
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In  calculation  of  Meissner  small-angle  scattering  and  small 
energy  transfers  are  examined  with  the  help  of  equation  of  transfer, 
and  large  energy  losses  with  the  help  of  Monte  Carlo  method. 

Probably,  this  methodology  will  make  it  possible  to  solve  many 
difficult  questions  of  theory  of  transfer,  in  particular  the  account 
of  fluctuations  in  energy  losses.  Meissner  calculated  in  the  barrier 
geometry  of  the  characteristics,  which  are  obtained  in  the  routine 
calculation  by  Monte  Carlo:  the  coefficients  of  the  transmission  and 
the  reflections,  energy  spectra  and  distribution  of  absorbed  energy. 
Our  calculations  satisfactorily  will  agree  with  Meissner's  results, 
with  exception  of  the  distribution  of  the  absorbed  energy.  Meissner's 
methodology  is  attractive,  since  in  it  it  is  possible  to  use 
advantages  of  analytical  approach.  However,  in  it  there  are  still 
many  vaguenesses,  including  how  to  make  more  precise  diagram  for  the 
case  of  the  heavy  substances,  when  radiation  losses  become  essential. 

In  conclusion  we  give  Table  10,  in  which  in  chronological  of 
order  are  equal  versions  of  basic  concepts  of  "amalgamated" 
collisions,  used  by  different  authors. 
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Table  10.  Versions  of  the  diagrams  of  the  "amalgamated"  collisions. 
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Key:  (1).  Author.  (2).  Energy  losses  to  ionization.  (3).  Energy 
losses  to  bremsstrahlung.  (4).  Angular  distribution.  (5).  Account 


of  secondary  radiation.  (6).  Geometry.  (7).  Use  of  methods  of 
accelerating  calculation.  (8).  Average/mean.  (9).  Moliere.  (10). 
One  measurement.  (11).  Taking  into  account  fluctuations  according  to 
Heitler.  (12).  Gaussian.  (13).  Continuous.  (14).  Moliere, 
Goudsmit-Saunderson.  (15).  Three  measurements.  (16).  Correlated 
tests.  (17).  Taking  into  account  fluctuations  according  to  Landau. 
(18).  Moliere.  (19).  Continuous.  (20).  Our  diagram  [89].  (21). 

Taking  into  account  fluctuations  according  to  Simon.  (22a). 
Average/mean  losses.  (22b).  Combined  distribution.  (23).  Gautsmit- 
Saunderson.  (24).  Analytically.  (25).  Splitting/fission  +  "tape 
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measure"  correlated  tests.  (26).  Monitor  -lizards- and  others  [90]. 

Page  75. 

S  4.  Methods  of  accelerating  the  calculation.  . 

It  is  known  that  success  of  application  of  Monte  Carlo  method  to 
problems  of  transfer  of  7-quanta  and  neutrons  to  a  great  extent  is 
connected  with  use  of  methods  of  accelerating  calculation.  We  have  in 
mind  both  the  optimum  organization  of  the  program  of  calculations  and 
the  special  procedures  which  ensure  the  efficiency  of  drawings,  in  such 
a  case,  when  in  the  direct  tests  only  the  insignificant  part  of  the 
histories  carries  the  information  interesting  us. 

Analysis  of  these  procedures  was  given  in  surveys/coverage  of 
Kahn  [76],  Kalos  and  Goertzel  [91],  in  monographs  [60,  64,  65],  in 
collectors/collections  [67,  68]  and  in  numerous  articles  of  colleagues- 
of  Novosibirsk  data  processing  center.  Let  us  note  that  many 
procedures  developed  for  7-quanta  and  neutrons,  can  be  transferred 
easily  to  the  electronic  trajectories  in  the  diagram  of  consecutive 
collisions.  Since  this  diagram  thus  far  still  rarely  is  used,  we  will 
consider  the  possible  methods  of  the  acceleration  of  calculation  in 
application  to  the  diagram  of  the  "amalgamated"  collisions.  First  of 
all  let  us  formulate  several  general/common  recommendations  regarding 
the  optimum  organization  of  program  for  computers. 

37c  /»  c./ejLsr7*)37' 

1.  ^  that  the  procedure  of  drawings,  repeated,  must  be  simple. 
This  it  is  possible  to  achieve  by  the  preparation  for  initial 
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distributions  in  the  form  of  numerical  matrices/dies  by  the  method, 
described  in  S2  of  this  chapter.  It  is  generally  necessary  to  attempt 
to  bring  in  into  the  memory  of  machine  all  distributions  and  other 
dependences,  to  which  it  is  necessary  to  be  converted  during  gaming 
procedure. 

2.  Algorithm  and  program  must  be  organized  so  as  to  ensure 
simultaneous  calculation  of  all  points  of  unknown  dependence.  Thus, 
in  the  problem,  in  which  is  investigated  the  dependence  of  the 
coefficient  of  transmission  on  the  thickness  of  the  sample,  it  is 
necessary  to  consider  simultaneously  the  passage  of  the  electrons 
through  several  intermediate  boundaries.  In  this  case  besides 
noticeable  gain  in  the  time  it  is  possible  to  decrease  the  statistical 
fluctuations  due  to  the  correlation  of  the  events,  which  occur  at 
adjacent  points. 

3.  Information  obtained  as  a  result  of  calculation,  must  require 
minimum  further  processing.  If  this  information  will  be  used  in 
further  calculations,  it  is  convenient  to  put  out  it  immediately  on 
the  punch  cards. 

Page  76. 

Large  time  gain  can  be  obtained  also  by  application  of  special 
criteria,  which  make  it  possible  to  break  "uninteresting"  trajectories 
(in  sense  of  interesting  us  information).  For  example,  in  the  problen 
of  backscatter ing  it  is  clear  that  only  the  electron,  which  penetrated 
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at  a  depth  less  than  its  residual  range  R0(E)  (E-energy  of  electron  at 
this  depth),  has  the  capability  to  leave  the  sample  into  the  rear 
hemisphere.  Consequently,  after  introducing  into  the  memory  of 
machine  the  table  of  values  of  the  path-energy,  we  can  cease  the 
tracking  of  electron  in  the  case,  "unfavorable"  for  the 
backscatter ing .  The  sequence  of  applying  this  criterion  is  shown 
schematically  in  figure  19.  Checking  is  produced  immediately  after 
the  drawing  of  the  new  value  of  energy  of  electron  and  checking 
E  -This  criterion  must  be  used  in  the  diagram  of  the  output  of 

electron  into  the  forward  half  sphere.  In  the  problems,  in  which  us 
interests  the  absorbed  energy  in  the  layer,  electronic  trajectory  must 
be  traced  in  practice  to  the  end/lead,  and  the  criterion  of 
cessation/discontinuation  is  applicable  only  to  layer  did  S*  i.  In 

this  case  into  the  memory  are  introduced  several  values  ' 

(according  to  a  number  of  layers,  if  they  different  thickness)  or  one 
value  Ed  -  the  energy,  necessary  for  electron  for  the  passage  of  layer 
d. 


In  Figure  20  it  is  shown  how  to  use  this  criterion.  Efficient 
proves  to  be  the  statistical  method  of  the  break  of  trajectory,  known 
by  the  name  "roulette"  [66].  Its  essence  consists  in  the  fact  that 
with  probability  p  the  electron  is  tracked,  and  with  probability  1  -  p 
tracking  terminates.  For  obtaining  the  unbiased  estimate  the  weight  of 
the  particle,  "which  survived"  after  the  application  of  a  "tape 
measure",  increases  from  1  to  1/p.  "tape  measure"  has  sense  to  use  in 


such  a  case,  when  analytical  criterion  is  not  yet  applicable,  but 
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trajectory  does  not  contain  necessary  information,  for  example 
electron  moves  in  parallel  to  boundaries.  As  far  as  the  selection  of 
probability  p  is  concerned,  its  evaluation/estimate  usually  is 
produced  from  the  condition  of  obtaining  the  minimum  of  expression 
D[$]*t,  where  D[£]  -  dispersion,  t  -  time  of  calculation. 


S' 
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Fig.  19.  Logic  circuit  of  use  of  the  criteria  of 
cessation/discontinuation  of  tracking  of  trajectory  on  leaving  of 
electron  from  sample. 

A 

Key:  (1).  determination.  (j!) .  no.  (3).  It  began  beam  of 

trajectory.  (4).  yes.  (5).  continuation  of  tracking. 

Page  77. 

Frequently  "tape  measure"  uses  together  with  methodology  of 
splitting/fission  (branching  [66].  This  methodology,  used  by  us  for 
the  first  time  for  the  electrons  in  the  problem  of  backward  scattering 
[92],  consists  of  the  following:  electron,  for  the  first  time  turned 
back/ago  (Fig.  21),  is  split  into  n  of  the  trajectories,  each  of  which 
has  a  weight  1/n  and  subsequently  is  outlined  independently.  The 
intersection  of  the  assigned  boundaries  of  layers  i  during  the  motion 
to  the  surface  of  barrier  leads  to  splitting/fission  n,.  moreover  the 
parameters  of  particles  (angle,  energy,  coordinate)  at  the  points  of 
splitting  are  memorized.  The  tracking  of  each  of  the  split 
trajectories  is  conducted  until  the  particle  exits  from  the  barrier  or 
before  the  absorption  in  it.  The  diagram,  used  us  for  the 
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memorization  of  the  parameters  of  particles,  and  also  the  procedure  of 
consecutive  return  to  the  points  of  splitting/fission  up  to  first 
point  of  rotation,  is  generalization  of  the  lexicographical  process, 
developed  in  the  work  of  Golenko  [71],  Since  we  are  only  interested 
in  the  back-scattered  particles,  tracking  of  the  particles,  which 
turned  into  the  depth  of  barrier  after  splitting/fission,  is  limited 
with  the  help  of  the  "tape  measure",  moreover  for  p(  was  selected 
value  p*  —  The  use  of  this  methodology  makes  it  possible  seemingly 
to  draw  out"  to  the  surface  the  electrons,  scattered  in  the  depth  of 
barrier  and  which  have  small  probability  of  exiting  leave  into  the 
rear  hemisphere.  Table  11  gives  the  results  of  the  calculations  of 
back-scatter  factor  for  different  orders  of  splitting. 

Order  of  splitting/fission  -sequence  of  values  rc,(i*l,  2,  3... 6 


in  our  case) . 
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Fig.  20. 


Fig.  21, 


fo  A® 

Logic  circuit  of  the  use  of  a  criterion  of  the 

cessation/discontinuation  of  the  tracking  of  trajectory  in  this  layer 

Key:  (1).  no.  (2).  yes.  (3).  energy  ...  is  absorbed  in  layer. 

(4).  in  roan  new  trajectory.  (5).  random  number.  (6). 
determination  of  scattering  angle. 

Fig.  21.  Schematic  diagram  of  application  of  methods  of 
splitting/fission  and  "tape  measure". 

Key:  (1).  keV.  (2).  tape  measure.  (3).  splitting/fission. 
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If  the  first  rotation  occurred,  for  example,  in  the  second  layer,  a 
number  of  splittings/fissions  is  taken  by  the  equal  to  2  and,  etc. 

The  efficiency  of  this  calculation  procedure  is  evaluated  according  to 
minimum  D[|]*t.  It  is  reached  at  the  order  of  splitting/fission  by 
3:2:1:1:1:1.  Gain  in  comparison  with  the  calculation  without  the 
splittings/fissions  composes  ~35%.  More  careful  research  of 
methodology  carried  out  in  the  work  Ogibin  (see  collection  [67]). 

Methodology  of  correlated  tests,  which  consists  in  use  of' one  and 
the  same  set  of  histories  for  research  of  two  close  processes,  for 
example  transfer  of  electrons  and  positrons,  proves  to  be  simple  and 
efficient.  Since  scattering  cross  sections  for  both  particles  barely 
differ  from  each  other,  the  with  the  limited  number  of  independently 
played  trajectories  difference  between  the  coefficients  of 
transmission  can  and  not  be  showed  due  to  the  large  dispersion  of 
calculation.  Let  we  seek  a  difference  in  the  coefficients  of  the 
transmission  6  of  positrons  and  electrons TU)  ?, ^  -  T N. Then  the 

dispersion  of  difference  can  be  evaluated  with  the  helD  of  the  formula 

)  -  2f- 1  \  M'.J  :  M'l],  (M* 

where  £[?]  and  £h]  -  dispersion  in  evaluation/estimate  t„  3-nJ  7\-  and  p  - 
correlation  coefficient.  When  p» 1,  then  the  dispersion  of  difference 
becomes  smali.  In  our  case  electronic  and  positron  trajectories  are 
actually  strongly  correlated,  since  entire  difference  between  them 
consists  in  the  small  shift  of  the  maximum  of  the  angular  distribution 
of  the  positrons  to  the  side  of  smaller  angles.  The  energy-loss 
distributions  are  identical.  We  used  this  methodology  for  the 
systematic  study  of  difference  in  the  transmission  of  positrons  and 
electrons  [93]. 
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±,’,,  =  0,6  Mm,  Z  =13,  q„  -  0J 
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Key:  (1).  Order  of  splitting/fission.  (2).  Number  of  trajectories. 

(3).  Number  of  particles  in  it  is  reflected,  spectrum.  (4).  (rel. 
unit).  (5).  ...  (without  splitting/fission). 

Page  79. 

Study  of  effect  of  angle  of  incidence  on  characteristics  of 
electrons,  passing  through  plane  layer  of  substance,  is  another 
important  application  of  methodology  of  correlated  tests.  In  this 
case  it  is  necessary  to  carry  out  parallel  calculation  according  to 
our  diagram  for  several  angles  of  incidence,  using  a  set  of  those 
played  one  time  of  parameters  £,A.  i.e.  to  actually  turn  trajectory 
around  point  x,=0,  yo=0,  z„=0.  Similar  methodology  was  used  for 
neutrons  in  the  work  Tyad'kin  et  al.  [94]. 

Experiment  of  application  of  Monte  Carlo  method  showed  that  it  is 
hardly  expedient  to  create  one  program  for  solution  of  many  problems 
of  transfer  of  electrons.  In  our  practice  we  use.  several  programs,  it 
is  sufficient  general/common  in  order  to  obtain  the  maximum  volume  of 


DOC  =  86063206 


PAGE 


information,  but  it  is  at  the  same  time  sufficient  narrow,  in  order  to 
for  certain  circle  of  the  problems  of  detecting/exposing  all  reserves 
for  the  acceleration  of  calculation.  The  characteristics  of  these 
programs  are  given  to  [95]. 


«t» 
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Chapter  IV. 

Results  of  calculations. 

In  the  present  chapter  are  presented  results  of  calculations, 
carried  out  by  Monte  Carlo  method,  given  are  data  for  coefficients  of 
transmission  and  and  reflection,  tables  of  spectra  after  barriers  of 
different  thickness,  distribution  of  absorbed  energy  in  plane  layers. 
These  data  are  suitable  for  the  interpolation  in  the  examined  region 
of  energy  of  electrons  0. 4-6.0  MeV  and  angles  of  incidence  of  0-81°. 
Are  described  the  procedures,  with  the  help  of  which  it  is  possible  to 
compute  the  unknown  values  and  the  spectra  for  the  nonmonoenerget i c 
composition  of  initial  flow. 

t 

[ 

S  1.  Integral  characteristics  of  the  passing  electrons. 

Integral  characteristics  are  important  for  estimation  of 
efficiency  of  protection  from  electrons,  for  calculations  of  output  of 
secondary  particles  from  protection,  etc. 

Coefficients  of  transmission  according  to  number  of  particles 
and  on  energy  Tt  are  functions  of  energy  of  electron,  atomic  number, 
thickness  of  the  layer  and  angle  of  incidence  and  are  determined  from 
relationships/rat ios 
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n  < x .  £.1  1  E\ 


TV  -  — 


V, 

,  /<  (Jr.  A'j 


(4.  li 
(4.  2) 


where  «(*,  £,)  -  number  of  electrons,  which  emerge  from  layer  of 

substance,  with  a  thickness  of  x  g/sm2  in  energy  range  E,~-  E,  r  li’,; 


N 0  and  E0  -  number  of  falling/incident  to  layer  electrons  and 
their  energy. 

Page  81. 


In  order  to  consider  statistical  error  during  calculation  of 
coefficients  Ts  and  7V  (i*.v  and  Rt), we  used  conventional  value  of 
standard  deviation  1 


3  --  q 


/ 

I  iVu  di1  ’ 


(4.  3) 


where  q  -  mathematical  expectation  of  random  variable  q; 


N 0  -  number  of  trajectories,  D[q]=q2-(q) 2 . 


FOOTNOTE  1 .  Less  commonly  for  the  evaluation/estimate  is  used  the 
value  of  the  probable  (middle)  deviation  E=0. 67456a.  ENDFOOTNOTE. 


Dispersion  of  value  q  depends  on  two  forms  of  errors: 
fluctuations  in  number  of  particles,  which  intersect  this  boundary, 
and  fluctuations  of  radiation  characteristics  Q  (energy,  direction). 
Then,  according  to  Berger  [60], 
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/  1  p  % 

q  ~  \  P‘V«"  («)J 

where  p  -  penetration  probability. 


In  the  case,  when  we  are  interested  in  the  error  in  coefficient 
of  transmission  according  to  the  number  of  particles,  Q  —  p  -  Ts ,  Q  1 


|  va 


Standard  deviation  during  the  calculation  of  the  coefficient  of 

transmission  according  to  the  energy  is  computed  from  the  formula 

^  ,  f  1 ~~Tn  £*-(£)*  ,, 

J  T*|  {i  b) 

Dependence  Tvu)  and  TjU)  is  a  smooth  curve,  which  collapses  with 
height  x.  If  we  isolate  on  this  the  curved  straight  portion  (i.e. 

•  d  T 

the  region,  where-—  -const)  and  to  continue  it  before  the 

dx 

intersection  with  axis  x  (see  Fig.  22),  then  the  corresponding  depth 
of  penetration  was  called  extrapolation  length  rNt. 


As  SELIGER  [96]  showed,  it  is  convenient  to  express  depth  of 
penetration  of  electrons  in  units  of  complete  path  R0  J. 


FOOTNOTE 


It  should  be  noted  that  many  authors  do  not  frequently 


give  values  of  R„ ,  as  a  result  of  which  the  uncerta inty/ indeterminacy 
upon  transfer  to  the  true  thicknesses  appears.  ENDFOOTNOTE. 
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In  this  case,  at  least,  for  the  energies  <1  MeV  the  experimental 
curves  of  transmission  for  the  different  energies  virtually  are 
combined.  The  same  conclusion  follows  also  from  our  calculations 
[89].  Subsequently  it  was  proposed  to  express  depth  of  penetration  in 
the  units  of  extrapolation  length.  Since  landing  run  r,  -  important 
characteristic,  you  give  the  empirical  formula  for  the  extrapolation 
length,  obtained  during  the  interpolation  designed  <■  N ,  for  elements  C, 
Al,  Si,  Ti,  Ge,  Sn,  Pb  in  the  region  of  wave  energies  from  0.4  to  6.0 

O 

MeV  and  angles  of  incidence  from  0  to  81  l 


rvAE  Z  .i.  \ _ 2^[  107.2- Z  ,  292,7 -Z 

’  ’  ™  2  [  5,4422-1312,0  4,1632+661,3 

,  2-2,797  _4 

1  63,862+587,6  E  J*  Z‘CM  »  .  (4*  7) 


where  4'o  -  angle  of  incidence; 

E  -  energy  of  electron  in  MeV; 

Z  -  the  atomic  number  of  element. 
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Fig.  22.  Dependence  r.v(*)  for  perpendicular  incidence/drop  in 
electrons.  Unbroken  curves  -  the  results  of  calculation,  point  - 
experiment.  1  -  Cu,  E,=6.0  MeV;  2  -  Al ,  Eo=4.0  MeV;  3  -  Pb,  Eo=1.0 
MeV;  4  -  Sn,  Eo=0.4  MeV;  n  -  scale  index. 
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By  analogy  with  extrapolation  length,  calculated  from  their 
curves  of  transmission  TN(x),  the  extrapolation  length  from  curves  of 
transmission  is  designed  according  to  energy.  Empirical  formula  for 


r*»  takes  the  form 


v  ,i.  \  ^  2—187,0  ,  1400,0— 2  n 

lAE,  +u)  -cob  2  [  2344,0-  7.626 Z  +  3T,822  +  3f08,Q  E  ^ 


2+14.14 

80,02+2207 


fer*']' 


(4.8) 


DOC  =  86063207 


PAGE 


The  values  of  extrapolation  length  obtained  by  us  will  satisfactorily 
agree  with  the  experimental  data  of  Ebert,  etc.  [97],  Kononov  [98], 
Radziyevskiy  and  Osanov  [99],  moreover  it  is  better  with  the  results 
of  the  first  work. 


% 


g! 


! 

i 


i 


From  the  data  of  "table  12  it  is  evident  that  calculated  by  us 
values  rN»  systematically  less  than  designed  by  Baranov,  etc.  [100]. 
Apparently,  this  difference  is  caused  by  the  fact  that  in  the 
latter/last  work  the  fluctuations  in  energy  losses  to  the  ionization 
and  the  bremsstrahlung  are  not  considered.  For  describing  the 
empirical  dependence  tn(x)  and  T*(x)  was  used  the  formula. 

2V*  =  exp|  —  P  »  (4-9) 

which  is  for  the  first  time  proposed  in  the  work  of  Makhov  [101].  The 
same  late  expression  was  used  by  Mar  [88].  We  obtained  the  following 
empirical  relationships/ratios  for  the  coefficients  a  and  /3  by 
working/treatment  according  to  the  method  of  least  squares  of  the 
curves  of  transmission.  For  TN 

1  2  o^o.oi2£  0,1564)  , 

+0,0125  (£-2)  (50  —  Z)  X 

X(k-0,1564;V'-™  Z<50,£>2  *4,  10) 
ci> 

l  +  fer1’  (u- 0,1564)  B  OCTB^HKX 

ciyiaaz, 


Key:  (1).  for.  (2).  inj^rest^/  (-9-K^Cases. 

6 

where  u  =  cos  tyo,  E  B  Mjb,  cos  81u  =  0,1564. 

M  2,59  -  0,0076  (Z  -  6),  Z  >  6. 


(4.  11) 


m 


■->  ■-«  >  V'.'.v,/,'.  /.  v.  ./t;  ■  ,,  v 


i:w»cw 


V 

V 
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Page  84. 
For  Ts 


a  =  0,708  -f 


Z+  24 

0,932+13,7 


(u-  0,1564). 


,  __  Z-32.6 

^  0,5242  +-10,3  • 


(4.  12) 
(4.  13) 


Calculated  according  to  empirical  formulas  (4.9-4.13)  curves  of 
transmission  TN(x )  and  TE(x)  will  be  agreed  on  with  experiment  (Fig 
22)  with  accuracy,  by  best  10%  for  angles  of  ^o<45°,  ,  and  depths,  for 
which  Tn.*(*)>0,1. 
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i  Table  12.  Comparison  of 


the  calculated  and  experimental  values  of 


extrapolation  length  rSjnc^  . 


s 

6  Mia 

4  Mia. 

2  Mia 

z 

a  8,2 

4  £<5> 

5  S— 

M  MS 

a  0-2. 

«s 

s  5 

a  3 

m  m 

n 

• 

0 

o  *  s 
A  A  o 

a 

SftS. 

5  a 
a  a 

<0  A 

B  et 

& 
2  •> 

®  V  r—y 

o  *  00 

x  S 

5  o.® 

<S?  a 

■  3 
a  £ 

•  A 

X  * 

h*1 

01  _ 

§§ 

aZ. 

3,6 

_ 

_ 

2,803 

1,882 

0,884 

1,129 

6 

3,160 

8. 192 

2,974 

1,963 

2,097 

l.oao 

1,002 

0,912 

1,020 

13 

2.967 

3,083 

2,910 

1,837 

2,020 

1,867 

0,960 

0,832 

0,874 

14 

2,942 

3,064 

2,897 

1,821 

2,009 

1,846 

0,964 

0,866 

0,861 

2,780 

2,936 

2.762 

1.706 

1,921 

1,740 

0,906 

0,777 

0.786 

32 

2,581 

2,778 

2.583 

1,680 

1,811 

1,613 

0,816 

0,729 

0,730 

SO 

2,267 

2,448 

2,287 

1,396 

1,613 

1,407 

0,738 

0,626 

0,887 

82 

1.881 

1,978 

1,861 

1,153 

1,281 

1,122 

0,546 

0,492 

0,606 

£ 

1  Mia 

0,7  Afj« 

0,4  Afj« 

Z 

5  a- 

*  is 

•  0.2, 

S  a 

a  a 

ad  A 

** 

•  __ 

*  3 

&! 

i  2 

5  ^ 
s  *s 

A 

wr 

x  3 

H  a 

«  40 

a  « 

<V _ 

S§ 

at 

s|- 

w  S  06 
5  3.8 

5  x 

3  x 

«  •» 

X  Bt 

A  „ 

S§ 

ad 

3.6 

.  i 
0,388* 

0,396 

0,512 

0,260 

0,245 

0,330 

0,127 

0,114 

0,166 

6 

0,464 

0,413 

0,463 

— 

0,266 

0,299 

— 

0,117 

0,141 

13 

0.430 

0,388 

0,396 

— 

0,248 

0,266 

— 

0,107 

0,120 

14 

0,426 

0.382 

0,391 

— 

0,243 

0,252 

— 

0,106 

0,119 

22 

0,398 

0,353 

0,367 

— 

0,223 

0,230 

— 

0,0969 

0,108 

32 

0,384 

0,320 

0,331 

— 

0,202 

0,214 

— 

0,0868 

0,101 

50 

0,300 

0,273 

0,303 

— 

0,172 

0,196 

— 

0,0732 

0,0920 

82 

0,188 

0,218 

0,2741 

— 

0,141 

o.m 

— 

0.06841 

0,0833 

Key:  (1).  (normal  incidence).  (2).  experiment.  (3).  our  data. 

(4).  calculation. 


FOOTNOTE  1 .  The  value  indicated  is  undertaken  by  IS  of  work  [99]. 
ENDFOOTNOTE . 
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For  large  angles  and  large  depths  approximation  (4.9)  gives  error 
20-30%.  Formulas  (4.9-4.13)  describe  the  coefficients  of  transmission 
only  for  ZtS .  For  the  light  substances,  for  example  Plexiglas,  we 
give  the  appropriate  values  of  these  coefficients  in  the  Table  9 
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Appendices . 


S  2.  Integral  characteristics  of  the  back-scattered  electrons. 


Backscatter ing  of  electrons  is  of  interest  for  many  branches  of 
nuclear  physics  and  technology.  In  this  direction  many  experimental 
works  are  carried  out.  The  first  data  on  back  scattering  by  the  Monte 
Carlo  method  are  acquired  in  the  work  of  Berger  [102].  As  in  the  case 
of  passage,  they  are  frequently  interested  in  the  integral 
coefficients  of  reflection  R N  and  rk,  which  are  defined  as  analogous  to 
T.v  and  Tt,. 


We  will  here  speak  only  about  complete  coefficients  of  reflection 
Rs  and  Rs,  i.e.  about  coefficients,  which  correspond  to  reflection  of 
electrons  from  semi-infinite  barrier. 


Curves  rse  monotonically  increase  from  0  to  value,  which 
corresponds  to  complete  reflection  coefficient.  Saturation  is 
achieved  actually  at  depths  x~0.3  Re.  We  did  not  attempt  to  describe 
designed  [92]  data  for  Rm  and  Ht  by  empirical  formula.  These  data  in 
the  form  of  graphs  for  different  angles  of  incidence  and  wave  energy 
of  electrons  are  shown  in  figures  23  and  24,  in  which  points  gave  the 
averaged  experimental  results,  obtained  for  the  perpendicular 
incidence/drop  in  the  electrons. 


f Y-  .Vv  • .  . 
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Fig.  23.  Countesses  of  dependence  rw  (z)  for  different  energies  and 
angles  of  incidence  i- 
Key:  (1).  MeV. 

Page  86. 

The  statistical  accuracy  of  the  results  of  calculation  was  not  less 
than  10%,  to  the  decrease  of  the  fluctuations  of  results  contributed 
the  methodology  of  "splitting/fission"  and  "tape  measures".  It  is 
possible  to  consider  that  the  results  of  calculation  will  completely 
satisfactorily  be  coordinated  with  the  experiment. 

Special  feature  of  electronic  trajectories  is  such,  that  an 
electron,  which  has  turned  back  at  small  depth,  as  a  result  of 
"successful"  selection  can  leave  sample.  The  probability  of 
scattering  to  the  large  angle  is  small  and  therefore  the  results  of 
calculation  will  be  sensitive  both  to  the  form  of  the  section  in  this 
region  of  angles  and,  possibly,  to  the  selected  diagram  of  Monte 
Carlo.  In  Table  13  we  are  congruent/equate  the  results  of  our 
calculations  rn  4  with  the  results  of  calculating  Berger  [102]. 


*- ■*  s*  «  ■  *  *  v.  * 
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From  the  data  of  Table  13  it  is  evident  that  results  of  both 
calculations  satisfactorily  will  be  coordinated  with  each  other, 
although  they  are  carried  out  with  the  help  of  different  diagrams  of 
Monte  Carlo.  Meanwhile  the  attempts  to  spread  these  calculations  to 
the  high-energy  region  lead  to  essential  divergence  R s  from  experiment 
[130],  as  is  evident  from  the  data  of  Table  14. 
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Page  88. 


As  can  be  seen  from.2able  15,  E  monotonically  they  grow  with 
increase  of  Z  and  «fu  only  sometimes  are  small  divergences  from 
monotonicities ,  caused  by  statistical  straggling. 


For  obtaining  reliable  data  about  angular  and  energy 
distributions  of  electrons  reflected  is  required  too  long  a  time,  even 
if  we  use  methods  of  accelerating  calculation. 


S  3.  Energy  and  angular  spectra  of  the  electrons,  which  passed 
through  plane  layer. 

Energy  and  angular  spectra  of  electrons,  which  passed  layer  of 
substance,  contain  much  greater  information  than  integral 
characteristics.  The  general  view  of  the  energy  spectra  of  electrons, 
passing  through  the  consecutive  layers  Ge,  it  is  shown  in  figure  25. 
Their  form  varies  from  the  curve  almost  gaussian  with  the  small 
thicknesses  to  the  wide  of  asymmetric  form  with  the  large  thicknesses. 
Unfortunately,  we  did  not  succeed  in  fitting  the  empirical  formula, 
which  describes  these  spectra  at  any  depths,  with  a  sufficient 
accuracy . 

Page  89. 

The  formula  of  the  form 

i(E)dE  =CE°i  Em-EY  dE,  (4.14) 

where  Em  -  maximum  energy  of  electrons  in  the  spectrum,  and  a,  /3  and  C 
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are  functions  Z,  E, ,  >fo,  it  can  be  used  only  on  the  large  depths,  but 
not  for  any  thicknesses,  as  this  is  indicated  in  work  [103]. 


As  sufficiently  rough  approximation  of  spectra  of  fogs  [88]  it 
proposed  expression 


n  ( E ,  E0,  Z,  .*)  dE  --- 

!D 

Aexp[ — b(Ea  —  E)]dE  npn  E  >En, 
0  E  ,  Em. 


(4.  15) 


Key:  ( 1 ) .  with. 

Parameters  b  and  A  are  determined  from  the  relationships/ratios 

(1.53  -  0,0147 Z)IE0,  A  = 

where  7V  -  coefficient  of  transmission. 


Our  results  are  represented  in  the  form  of  matrix/die  of  values 
for  density  of  distribution  n(E)  (f/Mev),  designed  for  several 
thicknesses,  maximum  from  which  correspond  7\v  (*m)  =  0.1,  and  10  values  of 

E 

energies  on  energy  scale  s- ,  where  Em  (MeV)  wall  energy  in 


spectrum,  and 


£*  ^  4?(2rt-rl),  h— 9,  8,  7 . 1,0. 
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Fig.  25.  Energy  spectra  of  electrons  E0=1.5  MeV  after  passage  of 
layers  Ge  of  different  thickness:  1  -  80.7  mg/cm2;  2  -  106.0  mg/cm2; 
3  -  129.4  mg/cm 2 . 

Key:  ( 1 )  .  rel .  un. 

Page  90. 

Values  n(£„,  x)  for  different  substances,  wave  energies  and  angles  of 

incidence  of  electrons  together  with  wall  energies  E  and  thickness  of 

the  layer  x  <^m2  are  given  in  the  Table  10  applications/appendices. 

All  spectra  are  calibrated  to  one.  For  the  production  of  a  spectrum 

taking  into  account  the  falling/incident  electron  values  n  (£’„)  from 

this  matrix  row  must  be  multiplied  by  T„(x).  It  is  natural  that  the 

accuracy  of  reproduction  of  the  true  spectra  with  the  help  of  the 

tables  is  low  for  the  small  thicknesses  (1-2  table  rows),  since  in 

this  case  in  one  interval  4>_*  almost  entire  spectrum  is  included. 

xo 

However,  the  following  lines  of  matrix/die  describe  the  spectrum  with 
the  satisfactory  accuracy. 
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In  figure  26  calculated  energy  spectrum  is  equal  with  spectrum, 
obtained  under  the  same  conditions  on  beta-ray  spectrometer  by  Andreev 
[104].  Agreement  is  very  good.  A  good  agreement  is  obtained  also 
during  the  comparison  of  calculation  from  experiments  [105]  for  Si 
(E„=0.4  MeV).  The  majorities  of  the  authors,  who  investigated  the 
spectral  composition  of  the  passing  through  plane  layer  electrons, 
usually  give  these  characteristics:  the  probable  losses  of  energy 
\E  -- E0~~  E.  and  average/mean  losses  AE=E0-E,  where  E.  -  the  most 
probable  energy  in  the  spectrum,  and  E  -  medium  energy  of  the  passing 
electrons . 

Table  16  gives  these  character ist ics ,  obtained  from  experiments 
of  Kononov  et  al.  [106]  and  from  our  calculat ions . 

Small  systematic  disagreement  between  calculated  and  experimental 
values  is  caused  by  differences  in  energies,  for  which  are  carried  out 
calculations  (Eo=6.0  MeV)  and  experiment  (E0=6.7  MeV). 

Relation  k  =  ,  according  to  data  of  work  [106],  weakly 

increases  from  k=1.15  (Z=6)  to  k=1.60  with  (Z=83). 
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Fig.  26.  Comparison  of  experimental  and  calculated  spectra  of 
electrons  (E„»0.624  MeV)  after  passage  of  layer  Au  with  thickness  of 
38.35  mg/cm2.  The  instrument  width  of  line  is  to  the  right  shown,  A  E 
-  probable  losses  of  energy. 

Key:  (1).  rel.  un.  (2).  keV. 

Page  91. 

According  to  our  data,  for  the  energy  E„=6.0  MeV  k=1.20±0.10.  In  the 
work  of  Baranov,  etc.  [107]  (E,S3  MeV)  k=1.18±0.24  does  not  depend  on 
Z  and  x.  Our  results  in  this  range  of  energies  do  not  confirm  this 
conclusion,  since,  first,  relation  ^  is  function  x  and  has  a  maximum 
for  x/R„ =0 . 15-0 . 03 ,  and  values  k  in  the  maximum  vary  from  1.20  for 
carbon  to  1.80  for  tin;  secondly,  with  an  increase  in  thickness  k-*l. 
This  disagreement  is  the  corollary  of  disregard  into  [107] 
fluctuations  in  energy  losses. 


The  data,  given  in  Table  10  Appendices,  can  -be  used  also  for 
determining  average/mean  from  energy  spectrum  of  electrons. 


A 
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Calculated  by  Monte  Carlo  method  angular  distributions  completely 
satisfactorily  will  be  agreed  on  with  experiment  [108,  109].  In 
figure  27  is  carried  out  the  comparison  of  calculation  with  the 
experimental  data  of  Lonergan  et  al  [109]  for  Al  and  Au  with  the 
energy  of  the  falling/incident  electrons  E0=4  MeV. 

Angular  distribution  with  small  thicknesses,  as  show  experiment 
and  our  calculation,  is  satisfactorily  described  by  Moliere's  curve, 
if  we  introduce  small  corrections  as  this  recommended  in  work  [108]. 
For  the  large  depths  the  distribution  of  electrons  approaches  the 
diffuse: 

f  (0)  —  cob2  6. 
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Table  16.  Comparison  and  a£„4c« 


Cu(Qe)' 


a\  v/  CM*, 
at  J 
^£p«c<i  ^ 


_ 1 _ _ 

m\  1.78  1.25  0,531,89  1,35  0,54  0,27  1,34  0,90  0,45  0,27 

Afj0<*  3,50  2,45  1,07  3,05  2.85  1,12  0,53  2,82  1,72  0,88  0,56 

„  Mm <£|3, 85  2, 50  1,1013.86  2,75  1,16  0,56  2,80  1,76  .0,88  0,50 

I  i 


Sn 

Pb 

X,  i,CM“ 

1,08 

0,83 

0,51 

0,24 

1.13 

0.87 

0,60 

0,25 

2,40 

1,75 

1,04 

0,48 

2,77 

2. 0B 

1,40 

0,54 

2,26 

1,75 

1.10 

0,46 

2,85 

2,20 

1,50 

0,66 

Key: 


MeV. 


e>p€r\rne.nta,  I  vaJoes  — -f0*'  Ci> 

FOOTNOTE  1 .  Computed  values  are  obtained  for  Ge(Z=32)^  For  Cu 
(E#=2 . 65  MeV)  we  obtained  a  good  agreement  of  probable  loss  with  the 
experimental  data  of  work  [48].  ENDFOOTNOTE. 
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Obtained  in  calculation  angular  distributions  for  x/Ro>0.4  are 
described  well  by  this  distribution  (unbroken  curve  in  Fig.  27) 

f  (6)  sin  6 d()  —  3  cos2  0  sin  (4.  16) 

(function  f(d)  it  is  calibrated  to  one).  In  this  case  the  mean  angle 
of  scattering  proves  to  be  constant  and  equal  to 

~z 

j’  0/  (6)  sin  OdO  = 
o 

=  -§-•  <4-17) 

Actually,  the  results  of  our  calculation  give  value  of  0=35°-  40°  for 
the  mean  angle  independent  of  the  wave  energy  of  electron,  atomic 
number  and  angle  of  incidence.  However,  to  a  similar  result  comes 
Baranov  [107].  Other  by  Berger  proposed  expressions  for  the  function 
distribution 

2  cos  & 

/(CO  -1,445  (0.717  cos  0 

r  cos-'  0)  (4.  18) 

it  leads  to  an  increase  in  mean  angle  of  9  from  35  to  42°,  that, 
obviously,  it  is  not  substantial  with  the  evaluation  of  not  very 
important  in  the  problems  of  transfer  characteristic. 

S  4.  Distribution  of  absorbed  energy. 

Beams  of  accelerated  electrons  began  in  recent  years  to  play 
large  role  in  national  economy.  With  their  aid  ultrapure  metals  and 
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Fig.  27.  Angular  distribution  of  electrons  (E0=4  MeV),  which  passed 
layers:  x*1.275  g/cmJ  A1  (histogram  1  -  calculation,  light  small 
circles  -  experiment);  x=1.460  g/cm2  Au  (histogram  2  -  calculation, 
dark/nonluminous  small  circles  -  experiment).  Continuous  curve  - 
f(0)=3  cos20. 

Page  93. 

Finally,  one  additional  very  important  application  -  radiobiology. 

For  the  practical  applications/appendices  it  is  important  to  know  the 
spatial  distribution  of  absorbed  energy,  the  appearing  as  corollary  of 
two  processes  -  energy  losses  to  the  ionization  and  the  excitation  of 
electron  shell  and  multiple  scattering  of  electrons  in  the  substance. 
The  distribution  curves  of  absorbed  energy  at  the  angles  of  incidence 
of  electrons  of'^£458  have  characteristic  bell-shaped  form  with  the 
maximum,  which  lies  on  depth  x~0.3  R,.  The  appearance  of  this  maximum 
is  caused  by  the  contribution  to  the  energy,  absorbed  in  this  layer 
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both  from  the  passing  electron  stream,  and  from  the  flow,  reflected  of 
the  deeper  layers.  Obtained  in  the  experiment  of  the  spatial 
distributions  of  absorbed  energy  -  difficult  problem.  Difficulties 
are  connected  with  the  inadequacy  of  the  methodology  of  the 
determination  of  absorbed  energy  with  the  help  of  ionization  chambers, 
especially  for  the  substances  with  large  Z.  Therefore  at  present  the 
majority  of  experimental  data  is  obtained  for  light  elements  [110]  and 
for  the  so-called  "tissue-equivalent"  substances.  For  the  heavy 
elements  the  most  precise  measurements  of  the  distributed  absorbed 
energy  were  made  by  Nakai  [111].  Is  recently  developed  new 
methodology  [112],  which  over  the  long  term  will  make  it  possible  to 
obtain  the  distributions  of  absorbed  energy  with  the  satisfactory 
accuracy  for  any  substances.  But  thus  far  calculated  methods  remain 
fundamental  for  obtaining  these  distributions. 

Processing  our  calculation  data  made  it  possible  to  obtain 
empirical  expression  for  distribution  in  depth  ~~  of  absorbed  energy 

•V  9 

^_.4[1  +  a(^)*]r„(4).  Sep.  (4.i9, 

Here  AE(x)  -  energy,  absorbed  in  layer  Ax  at  depth  x  in  the  sample; 

Ts  —  coefficient  of  transmission,  computed  from  formula  (4.^.9). 
Constants  A,  a,  b  are  functions  E,  Z,  u  for  the  elements  14SZ£50  are 
located  with  the  help  of  the  formulas 

A  '  5'17(1  -»>;  6 “[1.10(2-1*)- 0,43]-*; 

.  _.,K  f  1.21  r(l  *01  (Z  -  14) 

38,78  +  392, 14  (1  -u)  ’ 

where  E  -  energy  of  electron  in  MeV,  u=cos  | ^  The  accuracy  of  the 
description  of  calculated  distribution  by  formula  (4.19)  composes  ± 
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20%)  for  angles  of  and  about  30%  for  angles  of  45o^0<63°. 

Page  94. 

The  agreement  of  calculated  distributions  AE(x)/Ax  for  Al  and  Pb  (E„-l 
MeV)  with  experimental  results  of  [111],  as  can  be  seen  from  figure 
28fAvery  good.  For  Plexiglas,  important  simulating  connection  in  the 
radiobiology,  we  computed  the  distribution  of  absorbed  energy  [113]. 
They  are  given  in  figures  29-34.  Transition  to  the  absorbed  dose 
(caused  by  flow<2>=!  cm  J )  is  realized  the  multiplication  of  absorbed 
energy  to  the  constant  factor  1.6*10*"*  rad/Mev. 

Knowing  values  of  absorbed  energy  on  this  depth,  it  is  possible 
to  compute  such  important  characteristics  as  value  of  the  induced 
charge  in  dielectrics  and  quantity  of  surplus  carriers  (electrons  and 
holes)  in  semiconductor. 

Recently  Kobetich  and  Katz  [114]  proposed  an  empirical  method  for- 
determination  of  distribution  AE(x)/Ax. 


- ■  .A -i  -■  -  ■  -  ■  i  i M  n  1 1 1  fi i»  il  i 
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Pb.  Eo=1.0  MeV.  Unbroken  curves  -  calculation,  points  -  experiment. 
Key:  (1).  MeV  cm3.  (2).  g/cm* . 
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Calculations  according  to  this  method  are  connected  with  the  numerical 
differentiation  of  a  certain  function,  which  is  the  product  of  the 
coefficient  of  transmission  Tn  and  polynomial,  which  expresses  the 
connect ion/communicat ion  of  residual  range  with  the  energy.  Although 
obtained  thus  distribution  will  agree  well  with  the  experiment,  it  is 
unclear,  as  it  is  possible  to  spread  this  method  to  the  case  of 
oblique  incidence  in  the  electrons. 

S  5.  Application  of  the  obtained  results. 
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Results,  presented  in  previous  paragraphs,  can  be  used  for 
solving  vide  circle  of  practical  problems.  Let  us  name  some  of  them 
here  only. 


1.  Determination  of  characteristics  of  passing  electron  stream 
for  polyatomic  systems,  in  particular  organic.  Calculations  showed 
that  in  this  case  it  is  possible  to  use  the  systems  described  above  of 
Monte  Carlo  for  a  certain  hypothetical  monoatomic  substance  with  the 
effective  values  of  atomic  number  2.^,  of  atomic  weight  .4.4,11  and 
ionization  potential 
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Fig.  29-34.  Distribution  in  depth  of  absorbed  energy  in  Plexiglas  for 
different  wave  energies  E,*6,  4,  2,  1,  0.7,  0.4  MeV,  respectively  and 
angles  of  incidence:  1  -  vo  =  0°;  2  -  *„*27°;  3  -  i\«*450  ;  4  -  s.=63°; 

-  *,-81°. 

Key:  (1).  MeV  cm.  (2).  cm2. 
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These  values  are  computed  from  the  formula: 

(4.20) 

i 

where  M  =  A,  Z,  I;  Mt—  their  partial  values; 

—  the  concentration  by  weight  of  this  component  in  the  system. 
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From  Table  17,  where  are  given  effective  values  2,  A,  I  for 
series/row  of  organic  matter  and  water,  it  is  evident  that  they  differ 
little  from  each  other  (with  exception  of  polyethylene).  Therefore 
data  of  calculations  for  Plexiglas  can  be  used  also  for  the  series/row 
of  other  substances,  in  particular  for  the  animal  cloth. 

Calculations  show  that  coefficients  of  transmission,  found  thus 
for  Nal  (ZJlN,  32),  also  well  will  agree  with  experiment. 

2.  Determination  of  coefficients  and  Tn  for  arbitrary  energy 

and  angular  distribution  of  electrons  2ir  f(E,  u)  dE  du.  It  is  clear 
that  in  this  case 

■^m»x  0 

T„.E -2r.  \  i'r.v.i  iu,  E,  Z,  x)  f  (E,  u)dEclu,  (4.21) 

where  u  -  co_8 'K;  £<nm  and  wall  energies  in  the  spectrum.  We  in  work 

[115]  have  carried  out  calculations  ^Vtfor  the  spectrum  of  the 
electrons  of  Earth  radiation  belt. 

3.  Determination  of  approximate  number  of  secondary  electrons, 

dislodged/chased  from  barrier  by  thickness  x  g/cma  into  front/leading 
hemisphere:  a)  by  7-radiation  and  b)  by  fast  electrons  during 

perpendicular  incidence/drop. 

a.  Thickness  of  sample  is  assumed  to  be  equal  to  landing  run  of 
photoelectron  with  energy  E=hv-I  (I  -  electron-binding  energy  in 
K-layer),  since  with  this  thickness  flow  of  secondary  electrons 
achieves  saturation. 
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Table  17.  Effective  values  Z,  A,  I. 


f?  1 

Bemecrso  j 

<PopM  y.iai 

IljieKCHrjiac^3) 

C„H,0, 

3,6 

6,67 

46, a 

UoJIHCTHpOJlW) 

C,H, 

3,36 

6,02 

44,5 

rlOJlli3TMA«H(r} 

CH; 

2.67 

4,67 

36,3 

Bcu»  (  t) 

H.O 

3,3 

6,0 

40,0 

Key:  (1).  Substance.  (2).  Formula.  (3).  Plexiglas.  (4). 

Polystyrene.  (5).  Polyethylene.  (6).  Water. 

Page  97. 

The  output  of  secondary  radiation  was  studied  in  work  [116],  and  fo 
713  is  obtained  the  following  expression: 

d  o 

•  -  ~  ~  \  fe  x  {Z  jJ-  Ts  [£*(<?).  2.  u]  *- 
0  1 

7\v  [£*<*),  Z,  x.  u]\dudx,  (4.  22) 

where  A  -  atomic  of  element; 


v  -  coefficient  of  absorption  'th 


x  -thickness  in  g/cm2; 


the  differential  cross  section  of  the  Compton  effect, 
determined  by  Klein-Nishina-Tamm  formula; 


the  differential  cross  section  of  photoelectric  emission 
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£,.('{'),  £+{y>  -  corresponding  energies  of  a  Compton-electron, 
scattered  at  angle  4.  and  the  photoelectron,  dislodged/chased  in  the 
direction  4 '• 

Tv  coefficient  of  transmission. 

In  formula  (4.22)  integrated  on  azimuth  angle  is  already  carried 
out.  The  calculations  of  the  output  of  secondary  electrons  for  the 
thickness,  which  corresponds  to  medium  energy  Compton-and 
photoelectrons  taking  into  account  the  static  weight  of  these 
processes,  are  carried  out  in  work  [1183.  In  this  case  the  thickness 
of  sample  will  be  somewhat  smaller  than  the  landing  run  of 
photoelectron.  Entire  process,  beginning  from  the  generation  of 
secondary  electrons  and  ending  with  their  tracking  in  the  sample,  is 
simulated  according  to  the  Monte  Carlo  method. 

Comparison  of  results  of  calculation  according  to  formula  (4.22) 
with  our  results  is  given  in  Table  18. 

If  we  take  into  account  aforesaid  above  about  different 
thicknesses  of  sample  in  both  calculations,  then  data  of  table  can  be 
considered  as  confirmation  of  method  of  analytical  calculation  of 
output  of  secondary  electrons. 

b.  By  analogy  with  7-quanta  it  is  possible  to  calculate  quantit 
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-r  — i  -x — ; 

- - d - i 

i 

Fig.  35.  Diagrammatic  representation  of  "birth"  of  6-electron  at 
depth  x,  g/cmJ  in  sample  with  thickness  of  d,  g/cml. 
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“min  d 


}  \T,lB.Z,x)i2§iFiTxlBW,Z. 


d  —  x,  u]dudx. 


(4.  23) 


where  “  section  of  Moller,  determined  by  formula  (1.23)  for  the 

medium  energy  of  electrons  in  spectrum  E(x); 


£(,),)  -  energy  of  a  6-electron,  which  moves  at  angle  to  the 

direction  of  motion  of  the  initial  electron  (it  is  determined  from 


formula  1.22); 


.  [  £+1.02  T* 

'  (  £+2,01  J  • 


and  E  -  energy  in  MeV. 


Some  results  of  calculations  according  to  formula  (4.23)  are 
given  in  work  [97] .  It  is  shown  that  the  noticeable  contribution  of 
6-electrons  to  the  coefficient  of  transmission  can  be  expected  only 
with  the  small  thicknesses  of  sample  and  for  the  light  elements. 
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4.  Distributions  for  absorbed  energy  can  be  used  for  research  of 
processes,  in  which  participate  fast  polyvalent  ions. 


Table  18.  Output  of  secondary  electrons  to  one  that  falling  of 
7-quanta . 


£.,  M,a  th 

X,'i  r.vi-  <V> 
><,,.1GS 

2,0 

0,564 

10,6 

1.5 

0,384 

7,1 

,  0.034  | 

|  1,55*)  i 

1.25 

0,282 

5.5 

0,614 
7,1»)  ! 

0.7 

0,095 

2.7 

j  Cu 

E  ,  M  m  0 

2.0 

1 

i  1.5 

1.25 

0.7 

X,  i  C.M-® 

0.700 

0.160 

0.046  | 

0,334 

0.467 

0, 114 

n,X  10-* 

7,8 

7.3 

1.06*  1 

4.6  I 

5,53* 

2,0 

Key :  ( 1 )  . 

MeV. 

(2). 

g/cm 

i 

• 

FOOTNOTE  1 .  Data  are  undertaken  from  work  [116].  ENDFOOTNOTE. 
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In  this  case  appears  a  large  number  of  6-electrons,  as  a  rule 
short-run,  and,  consequently,  the  track  of  this  ion  is  cylindrical 
tube  with  the  very  high  density  of  absorbed  energy. 


Information  about  this  energy  is  used  for  evaluation  of  these 
important  characteristics:  1)  light  output  from  ionic  crystals  during 
recording  of  heavy  ions;  2)  sizes/dimensions  of  track  in  nuclear 
emulsion;  3)  efficiency  of  track  defects;  4)  radiobiological  effect 
from  fission  fragments,  etc. 
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Work  [114]  gives  formula  for  energy  distribution  of  6-electrons, 
formed  on  1  cm  of  path  by  ion  with  effective  charge  (Ze)3^ 


d<o 


_3ajL_h.a3ijL)  •  (1  ..  ^ 

/«0  cm  l  ‘  "'/n  "'m  /  '°m  / 


,  (4.24) 


where  N  - 


number  free  electro?^*  in  1  cm3  of  substance 
(Ze)JlMl  Ze[  1  -  exp(  —  125iZ-‘‘)] 


and 


Zm  „  <•"  •«• 

a  r)  - 


Kinetic  energy  of  electron  E  is  determined  on  known  energy  of  binding 
I  [119]  of  the  relationship/ratio  E=<j-I.  From  the  kinematic 
considerations  it  follows  that  between  the  angle  of  emission  f  of  a 
6-electron  with  respect  to  the  trajectory  of  ion  and  its  energy  there 
is  a  dependence  cost=  (~j  \  Since  that  <|>  -  -jj-.  In  order  to  find 

the  distribution  of  absorbed  energy  along  the  track  of  ion,  it  is 
necessary  to  compute  the  expression 


~  -  j  f  (4.25) 


E  m  tm 

\  (IE  rr. 


d  ’  n 


I  0 


where  dE/dr  [/3{r),  E]  -  the  distribution  of  the  absorbed  energy  is 
determined  with  the  help  of  the  formula  (4.19jor  from  the  experiment, 
dJn/dEdr  it  is  obtained  from  (4.24)  taking  into  account  dependence 
/3=/3(r).  The  latter  is  found  from  relationship/ratio  a  landing 
run-energy  for  the  ion.  Limit  *<»  in  integral  (4.25)  corresponds  to 
maximum  energy  Em  - .  <»m  ...  /. 


5.  In  connection  with  requirements  of  radiation  physics  of 


DOC  =  86063208 


PAGE 


semiconductors  recently  especially  increased  interest  in  research  of 
radiation  defects,  formed  in  Si  and  Ge  by  fast  electrons. 


Pagw  100. 


This  is  explained  the  possibility  of  variation  over  wide  limits  of 
energy  of  electrons  and  introduction  thus  .o  the  sample  of  the  defects 
of  different  forms,  and  also  fact  that  the  electrons  of  not  very  high 
energy  (Ea<20  MeV)  virtually  do  not  induce  activity  in  the  irradiated 
object  and  experiment  is  pure/clean.  For  the  same  reasons,  which 
cause  the  distribution  of  absorbed  energy  in  the  sample,  appears  the 
spatial  distribution  of  defects.  Knowledge  of  this  distribution  is 
important  for  a  whole  series  of  applied  questions  of  radiation  physics 
and  technology. 


In  work  [89]  we  showed  that  it  is  easy  to  obtain  him,  combining 
results  of  calculations  of  spectra  of  electrons  at  given  depth  n(x, 

E) ,  obtained  by  Monte  Carlo  method,  with  experimental  cross-section 
a(E)  of  formation/education  of  defects.  A  quantity  of  defects,  formed 
in  the  sample  at  depth  x>xx  (it  is  examined  the  one-dimensional  case), 
it  is  possible  to  compute  according  to  the  formula 


N  (x  x,) 


1  u  A)  u  £. 


(4.  26) 


where  corresponds  to  the  threshold  energy  of  the 

formation/education  of  defects  (in  Si,  E=0.15  MeV); 


maximum  energy  of  electrons  in  the  spectrum. 
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In  order  to  find  distribution  itself,  it  is  necessary  to 


construct  dependence 


v  t  i  <L'  »x’ 
A  U)  ax~- 


Figure  36  shows  dependence  on  depth  of  probability  of  forming 
defects  77(x)~N(x)  in  Si  with  wave  energy  of  electrons  Eo=0.5  MeV.  In 
the  formula  experimental  section  from  the  work  of  Loferski,  etc.  [120] 
was  used.  The  results  of  calculation  (curve  1)  and  experiment  (curve 
2)  satisfactorily  will  be  coordinated  between  themselves.  Curve  3  is 
designed  Yurkov  [121]  with  the  help  of  the  analytical  method,  similar 
to  Spencer’s  method.  The  difference  between  both  calculations  should 
be  related  due  to  the  impossibility  in  this  stage  of  correct  account 
in  the  equation  of  the  transfer  of  fluctuations  in  the  losses  of 
energy  and  multiple  scattering. 


Page  101. 


Subsequent  calculations  [122]  are  carried  out  for  the  broader  band  of 
energies  of  electrons,  moreover  as  the  probability  of  forming  the 
defects  in  n-silicon  by  monoenerget ic  electrons  formula  [123] 


i](£)  —  0,18  In  AE,  , 


<4.  27) 


where  E  -  energy  of  electrons  in  MeV,  is  used.  Threshold  energy  is 
taken  in  this  case  as  the  equal  to  0.25  MeV.  The  results  of 
calculations  for  different  wave  energies  and  normal  incidence  are 
shown  in  Figure  37.  The  case,  when  silicon  is  covered  with  the  layer 
of  carbon,  is  investigated  also,  i.e.,  defect  formation  occurs  after 
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the  passage  of  the  layer  of  protection  x  g/cm2. 

For  convenience  thickness  of  sample  Si  or  protection  from  carbon 
is  expressed  in  units  of  complete  landing  run 

Ro  si  -0,1052  -r 
r  0.6784E  -  0,Q077£2. 

—  0,0950  + 
r  0,5994 E  -  0,005 5£2, 

where  E  -  energy  of  electron  in  MeV. 

It  is  obvious  that  calculated  distribution  in  silicon  can  be  used 
for  approximate  estimate  of  quantity  of  defects,  which  arose  after 
passage  of  protection  from  Al  with  thickness  of  x  g/cm2.  The 
procedure  of  calculation  of  defects  for  the  monoenerget ic  electrons 
easily  is  generalized  to  the  case  of  the  initial  spectrum  of 
electrons.  Furthermore,  instead  of  the  experimental  section  can  be 
used  theoretical  sections.  From  the  comparison  of  the  results  of  such 
calculations  with  the  experiment  it  is  possible  to  explain  the 
validity  of  one  or  the  other  theoretical  assumptions.  Let  us  note 
that  the  more  correct  calculation  of  distribution  of  defects  in  the 
sample  is  obtained  during  the. direct  calculation  of  a  quantity  of 
defects  in  each  segment  of  the  traced  trajectory. 

S  6.  Consecutive  account  of  secondary  radiation. 

Efficiency  of  diagrams  of  Monte  Carlo  examined  for  solving  of 
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great  circle  of  problems  was  demonstrated  only  in  application  to 
primary  emission  (electrons). 


V. 
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Fig.  36.  Distribution  in  depth  of  defects  tj(x)  in  Si,  formed  by 
electrons  Eo=0.5  MeV.  1  -  calculation  by  Monte  Carlo;  2  -  experiment 
3  -  analytical  calculation. 

Key:  (1).  rel.  un.  (2).  mm. 
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Of  no  less  interest  for  the  practice  are  the  problems,  in  which  is 
required  the  detailed  information  about  the  characteristics  of 
secondary  radiation  and,  first  of  all,  braking  and  X-radiation,  which 
emerge  from  the  target,  bombarded  with  electrons  1 . 


footnote  *.  We  do  not  examine  the  case  of  the  very  high  energies, 
which  initiate  electron-photon  cascade/stage.  Similar  problems  were 
solved  both  by  the  analytical  methods  [124]  and  with  the  help  of  the 
method  of  Monte  Carlo  [1253.  ENDFOOTNOTE. 
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It  is  known  that  the  complete  examination  of  these  problems 
contemporary  theory  is  difficult. 
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Fig.  37.  Distribution  of  probability  of  formation  of  defects  >\Rj  in 
Si  for  different  wave  energies  of  electron.  1  -  defects,  formed  in  Si 
under  layer  C;  2  -  defects  in  Si.  a  ~  Ee=6  MeV;  b  -  E0=4  MeV;  c  - 
E0  =2  MeV;  d  -  E,=l  MeV;  e  -  Eo=0.7  MeV;  f  -  Ee  =  0.4  MeV. 

Key:  (1).  defect/e.  cm. 
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Here  also  the  significant  role  play  boundary  conditions.  In  this  case 
in  the  thick  targets  is  substantial  repeated  scattering  and 
fluctuation  in  the  energy  losses  of  primary  electrons  as  well  as 
multiple  scattering  and  self-absorpt ion  of  secondary  radiation.  At 
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present  only  with  the  help  of  the  Monte  Carlo  method  it  is  possible 
obtain  in  entire  space  information  about  secondary  radiation  l. 


FOOTNOTE  l.  Most  efficient  in  this  case  proves  to  be  the  calculation 
procedure  according  to  the  method  "catastrophic"  collisions. 
ENDFOOTNOTE . 


The  report  of  Ruffle  [126]  is  one  of  the  first  works,  carried  out  in 
this  direction.  In  the  process  of  tracking  the  generation  of 
secondary  electrons,  positrons  and  7-quanta  was  examined.  For  the 
drawing  of  photon  energy  differential  cross  section  (1.31)  (case  of 
complete  shielding  7=0)  was  used.  During  the  conversion  of  7-quanta 
into  the  pair  electron  -  positron,  the  formation  of  Compton  and 
photoelectrons  the  drawings  were  produced  from  the  appropriate 
distributions.  For  the  photoelectrons  is  used  Sauter's  formula  [117] 
Drawing  of  energy  (angle)  of  the  secondary  electron  of  was  produced 
from  the  distribution  of  Moller.  The  tracking  of  electrons  was 
conducted  in  the  approximation/approach  of  the  continuous  energy 
losses,  angular  distribution  in  multiple  scattering  was  received  as 
Gaussian.  In  the  calculation  the  annihilation  of  positrons  and  the 
photoproduction  of  neutrons  were  taken  into  consideration  also. 
Obtaining  heat  distribution,  which  is  isolated  in  the  target, 
irradiated  by  electrons  with  the  energy  in  the  interval  of  50-150  MeV 
is  the  fundamental  result  of  calculation. 


Fundamental  research  of  applicability  of  Monte  Carlo  method  to 
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problem  in  question  was  carried  out  by  Berger  [127].  The  distribution 
of  intensity  calculated  by  it  in  the  spectrum  of  bremsstrahlung  in 
incident  direction  in  the  electron  is  equal  in  figure  38  with  the 
experimental  data  of  Rester  etc.  [128].  Agreement  is  very  good  and, 
which  is  especially  valuable,  even  in  the  region,  close  to  the  maximum 
photon  energy.  The  interesting  calculation  of  braking  electron  beam 
by  energy  to  60  MeV  in  the  thick  target  from  lead  led  bibliophiles 
[129].  The  drawings  of  energy  of  bremsstrahlung  were  conducted  with 
the  use  of  differential  cross  section  of  Shiff  in  the  form  (1.33), 
which  ensures  better  agreement  with  the  experiment  in  the  high-energy 
region  of  7-quanta,  they  are  close  to  the  boundary.  Calculation  was 
divided/marked  off  into  two  stages.  In  the  first  were  computed  the 
characteristics  of  the  generated  7-quanta  and  electrons  in  each  of  25 
layers,  to  which  was  decomposed  the  target.  Subsequently  the  tracking 
of  secondary  radiation,  which  leads  to  the  appearance  of  tertiary 
generation  of  the  particles,  was  conducted:  electron-positron  pairs; 
with  photo-  and  Compton-electrons ,  etc. 

Page  104. 

Electronic  (positron)  trajectories  were  outlined  just  as  into  [126], 
using  an  assumption  about  the  continuous  energy  losses.  The  designed 
heat  release  will  agree  well  with  the  experiment.  It  is  possible  to 
hope  that  these,  until  scarce,  calculations  are  soon  considerably 
expanded,  and  practice  will  obtain  the  reliable  recommendations, 
necessary  during  the  design  of  accelerators ,  conducting  the 
calculations  of  protection,  also,  for  purposes  of  technology.  The 
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problem  of  the  beam  shaping  of  monoenerget ic  highenergy  neutrons  is 
also  of  great  interest.  However,  similar  calculations  very 
labor-consuming  and  for  their  realization  are  required  the 
computer(s),  which  possess  the  large  memory  and  the  speed. 
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Table  1.  Mott's  factor;;*,  for  electron  scattering 


E0.  Mje  Q) 


10  ! 

4  ! 

2  1 

i  ! 

0.7 

0.4 

0 

1,000  | 

1,000 

1,000  ( 

1,000 

Z= 6 

1,000 

1 

1,000  i 

15 

0,999  j 

0,999 

1,000 

1,000 

1,001 

1,002 

30 

0,961  ; 

0,962 

0,963 

0,967 

0,970 

0,977  ; 

45 

0.888 

0,890 

0,894 

0,903 

0,911 

0,929  | 

60 

0,788 

0,790 

0,797 

0,814 

0,828 

0,880  ! 

75 

0,668 

0,669 

0,680 

0,705 

0,728 

0,776  | 

90 

0.532 

0,537 

0,551 

0,586 

0,617 

0,683  1 

105 

0,397 

0,403 

0,421 

0,466 

0,505 

0,590  j 

120 

0,269 

1  0,277 

0,298 

0,352 

0.440 

0.501  l 

135 

0.1591 

1  0,1680 

0,1923 

0,254 

0,308 

0,425  I 

150 

0.0742 

|  0,0839 

0,1106 

0.1787 

0,238 

0  366  i 

165 

0,0206 

1  0,0310 

0,0591 

0,1310 

0,1938 

0,328  ; 

180 

0,0024 

j  0,0129 

i 

0,0416 

;  0,1148  ! 

0,1786 

!  °'3“  , 

0 

1,000 

J.,000 

1,000 

1,000 

Z-13 

1,000 

1,000  J 

15 

1,020 

1,020 

1,020 

1,020 

1,019 

1,019  | 

30 

0.997 

0,998 

0.999 

1.001 

1,003 

1.008 

45 

0,935 

0,936 

0,939 

0.947 

0,954 

0,968 

60 

!  0, 83* 

i  0,840 

0,846 

0,881 

!  0.874 

,  0.903 

75 

0,714 

1  0,718 

0,727 

0,752 

0,773 

j  0.818 

90 

0.575 

0,580 

0,594 

0,628 

0,658 

i  0,722 

105 

0,431 

!  0,438 

0,455 

1  0,499 

!  0,538 

1  0,621  | 
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♦j- 

£...  M  ■- 

& 

1  1" 

4 

i  2 

1 

0.7 

0.  i 

0.2 

0.  > 

0.. 

120 

0.294 

0,302 

0,323 

0,377 

0.424 

0,525 

0,674 

0.809 

■  0.902 

135 

0,1742 

0.1831 

0,207 

0.270 

0,324 

0.441 

0,613 

0,769 

i  0.879 

150 

0,0813 

0,5912 

0,1179 

0,1868 

0,247 

0,376 

,  0,566 

0.738 

!  0,861 

165 

0.0226 

0.0329 

0,0614 

0,1342 

0,1680 

0,334 

;  r, 536 

0.719 

!  0.850 

180 

0,0024 

0.0131 

0,0419 

0,1163 

0,1811 

0,320 

0,526 

0,712 

0,846 

Z 

=  29 

0 

1,000 

1,000 

1,000 

1.000 

1,000 

1.000 

1,000 

1,000 

1.000 

15 

1,069 

1,068 

1,068 

1,066 

1.064 

1.069 

1.050 

1,038 

1.026 

30 

1.094 

1.093 

1,063 

1,092 

1,091 

1,089 

1,082 

1,070 

1.054 

45 

1,066 

1.066 

1,068 

1,072 

1,074 

1.080 

1,085 

1.082 

j  1,072 

60 

0,987 

0,988 

0,993 

1.C04 

1,013 

1,032 

1.057 

1.073 

1  1,076 

75 

0.865 

0.868 

0,876 

0,896 

0,914 

0,952 

1,004 

1.045 

I  1  J67 

90 

0.712 

0,716 

0,729 

0,761 

0,788 

0,847 

i  0,921 

1,003 

!  1.047 

105 

0.543 

0,649 

0,566 

0,610 

0,648 

0,729 

0,848 

0,952 

1,022 

120 

0,375 

0,383 

0.405 

0,460 

0,508 

0,611 

0,763 

0,900 

0,995 

135 

0,224 

0,234 

0.259 

0,325 

0,382 

0,505 

0,687 

0.852 

0,970 

150 

0,1052 

0.1157 

0,1446 

0,218 

0,282 

0,420 

0,626 

0.813 

!  0,950 

165 

0,0289 

0,0402 

0,0710 

0,1497 

0,218 

0,366 

j  0,587 

0.789 

0,937 

180 

0,0026 

0,0141 

0.0456 

0,1261 

0,1964 

0,348 

1  0,573 

0,780 

I  0,932 

2=60 


1,000  1,000  1,000 

1,121  1,116  1,109 

1,281  1,123  1,216 

1,290  1,284  1,279 

1,275  1,277  1,277 

1,184  1,197  1,208 

1,028  1,055  1.078 

0,825  0,869  0,907 

0.602  0,664  0,718 

0,390  0,467  0,635 

0,214  0,305  0,385 

0,0683  0.1980  0,285 

0.0581  0.1608  0,251 


1,000  1 

1,000  I 

[  1.000 

1,096 

1,072  i 

l  1,043 

1,197 

1,161 

1,113 

1,264 

1,232 

1,183 

1.276 

1,282 

1,228 

1,227 

1,246 

1,242 

1,126 

1,188 

1,226 

0,988 

1,101 

1,190 

0,832 

0,999 

1,142 

0,682 

0,899 

1,093 

0,537 

;  0,815 

1,053 

0,475 

1  0,760 

1,026 

0,446 

!  0,741 

1.017 

1.000 

1.018 

1.060 

1.120 

1.173 

1,207 

1.225 
1.230 
1.229 

1.225 
1.221 
1.219 
1.21" 


0 

15 

30 

45 

60 

75 

90 

105 

120 

135 

150 

165 

180 


0 

15 

80 

45 

60 

75 

90 

106 

120 

135 

150 

166 

180 


2-82 


1.000  ! 

1,000 

1,000  I 

1,000  ! 

1,000 

1,000 

l.COO  I 

1.127 

1.125 

1,120 

1,108  ! 

1,098 

1,074 

1.042  1 

1.358 

1.354 

1,344 

1,315 

1,290 

1,230 

1.133 

1,658 

1.653 

1,638 

1,599 

1,664 

1.479 

1,328 

1.918 

i  1.912 

1.897 

1,857 

1,819  ; 

1.728 

1.555 

2,044 

2,040 

2,029 

2,000 

1.971  | 

1.896 

1.741 

1,981 

1.980 

1,979 

1,974 

1,966 

1.936 

1.844 

1,726 

1,731 

1,745 

1,777 

1.801 

1.842 

1,859 

1.324 

1.335 

1,366 

1,444 

j  1,510  1 

1  1,640 

1.799 

0.855 

0,871 

0,924 

1,050 

1,159 

1.385 

1.698 

0.422 

0,446 

0,513 

0,683 

0,830 

1.143 

1.592 

0.1158 

0.1444 

0,222 

1  0.4 22 

1  0,595 

0.969 

1,514 

0.0068 

0,0368 

0,1187 

i  0,328 

1  0,511 

0,908 

1.486 

Z 

[i-  p2 

1,000 

1,000 

1,000 

!  i.ooo  | 

1,000 

1,000 

l.OO'1 

1.103 

1,102 

1,097 

I  1,086  1 

1 ,076  ' 

1,058 

1.038 

1.321 

1.316 

1,304 

i  1,270  j 

1,240 

1.174 

1,075 

1.702 

1,694 

;  1.674 

1,619 

1,569 

1,456 

1.267 

2,119 

2,110 

!  2,085 

2,019 

1.958 

1.815  ; 

3,561 

2.413 

2,410 

2,388 

2.326 

2,268 

2,127 

1.859 

2.482 

1  2,477 

2.465 

2,429 

2,393 

2,295 

-  2.083 

2,264 

2,266 

2,272 

|  2,284 

2,290 

2,280  1 

2,198 

1,797 

1,809 

1  1,840 

1,918 

1.981 

2,096 

2,206 

1.191 

1,213 

1,273 

1.423 

1,551 

1,809  ! 

2,141 

0,598 

0,629 

0.715 

0,932 

1,120 

1.510 

2,048 

0.1652 

I  0,203 

0,307 

1  0,572 

0.8f  . 

1.287 

1.970 

O.fOJl 

0,0496 

0,1600 

;  0,441 

0,685 

1.206 

1,940 

1,000 
1.024  l 
1.040  1 

1,157  ! 

1.336  ; 

1.518 
1.672 
1.786  I 

1.866  I 

1,920  j 
1,955  ! 

1.974  ! 

1.978  I 


1.000  I 
1.034 
1.005  ! 

1,080 
1,27  2  ! 

1,521 
1.770  | 

2.961 
2.177 
2.226 
2.436  i 
2.505  I 
2.523  I 


1.000 

1.024 

o.m 

1.023 

1.112 

1.267 
1.435 
1.614 
1.799 
1.978 
2.13C 
2.233 

2.267 


1,000 
1.029 
1.007 
0,978 
1  .038 
1,193 
1.413 
3 ,879 
1.979 
2.287 
2.561 
2.753 
2.826 
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Table  3. 

Average/mean  energy  losses 

to  the  ionization  over  a  unit  of 

distance 

-  —  (MeV*smVg). 

«  f  i 

Mie  O) 

3,33;  4,3116 


,8692 

.5032 

.1792 

,2235 

,9339 

.9479 

.8260 

.7449 

.6494 

,6210 

.5255 


2,2264 
2,1962 
1,9256 
1,6913 
1,7279 
1,5111 
1,5230 
1.4341 
9 


0,8  j  1,0 


2,0326  j  1,9583  1.9286  ;  1,9448  I  2,0 

2,0055  !  1,9325  1,9034  !  1,9202  |  1,9 


1,7633  I  1,7024 


,7 


1,5223 
3 


1,2925  !  1,2825 


1,2115  I  1,1816  |  1. 


1.1280  |  1,1018  1,0965  1.1394  i  1,1921 

1,0197  ;  0,9980  j  0,9945  1,0376  !  1,0877 


4,0 

5,0 

[  2,1026 

2,0310 

2,0767 

1,8065  ! 

1  1.8492 

1.2366 

1.1298 


1,7115 

1,5209 

1,7486 

1,5553 

1.5356 

1,5704 

1.2739  i  1,3057 
1.1650  i  1,1950 


Values  of  complete  landing  runs  R0 ,  g/cm: 

Afj«  ( 


2,0 

1 

3,0 

4,0  1 

5.0 

'  6.0 

0.8384 

1.4410 

1,9278 

2.4012 

2,8629 

0,9507 

1,4595  i 

1,9520 

OJ 

o 

V 

C4 

2.8979 

1,0759  j 

1,6456 

2,1944 

2,7255 

3.2415 

1.2046  ; 

1,8274 

2,4203 

2,9877 

3,5333 

1 , 1762 

1,7823  ! 

2,3582 

2.9064 

3  436o 

,3253 
,3091 
,37E6 
,4347 
1 ,4846 
.4980 
,5761 
.6784 


1,9944 

1.9643 

2,6605 

2,1395 

2,1960 

2,2081 

2,3160 

2,4200 


2.6224 
2,5767 
2.6943 
2,7921 
2,8457 
2,8530 
2,9846 
3, 0688 


3,2158 

3.152; 

3,2374 

3.40C6 

3,4443 

3.4439 

3.5948 

3,6442 


3,7794 
3,9974 
3 , 8458 
3,9715 
4.0000 
3,9894 
4,1561 
4,1604 
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Table  7.  Values  of  Moliere-Bethe  functions, 


2 

1,9218 
1,7214 
1,4094 
1,0646 
0,7336 
0,4738 
0.2817 
0,1546 
0.0783 
0,0386 
0,01681 
0.00630 
0,00232 
0.00079 
0,000260 
7,3  10~D 
1.9  •  10“5 
4.7  10 

1.1  10~6 


0,8466 
0,703 8 
0,3437 
-0,0777 
-0,3981 
-0,6286 
-0,4770 
-0,3183 
-0,1396 
-0,0006 
0,0782 
0,1054 
0,1008 
0.08262 
0,08247 
0,04660 
0,03288 
0,02402 
0,01781 
0,01386 
10,838*10“3 
8.14010'3 
3,831*  lO"3 
2,527*  10~3 
1,739*  10"3 
0.908*  10~3 
0.521*10“3 
0,3208*  10  ~3 
0,2084*  10-3 


2,4929 
2,0694 
1.0488 
-0,0044 
-0,6068 
-0,6369 
-0,3088 
0,0525 
0.2423 
0,2386 
0,1316 
0,0196 
0.0467 
0,0649 
—0,0646 
-0,03568 
-0,01923 
—0,00847 
-0,00284 
0,00006 
'1,0741*10“* 
1.2294-10-3 
0,8326*  10~* 
0,5368*  10~3 
0,3495*  10“3 
0,1584*10~3 
0,0783*  10“3 
0,0417*  10“3 
0,0287*10  3 
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Table  8.  The  angular  distribution  of  multiple  scattering  in  the 
theory  of  Nigem,  etc.  [58]. 


4  1 

5  ! 

8  i 

7  1 

8  1 

9  i 

10  i 

12 

0,0 

1,00000 

1,00000 

l 

1,00000 

1,00000  1 

1,00000 

1,00000 

1,00000 

1.00000 

0,2 

0,94070 

0.94546  ' 

0,94850 

0,95058 

0,95203 

0,95321 

0,95409  j 

0.95537 

0,4 

0,78389 

0,79992 

0,81017  i 

0,81721 

0,82232 

0,82616 

0,82910 

0.83351 

0.6 

0,58102 

0,60800 

0,62535  i 

0,63731 

0.64601 

0,65259 

0,65772 

0.66520 

0,8 

0,38726 

0,41889 

0,43939 

0,45363  | 

0,46302 

0,47192 

0,47811 

0,48716 

1,0 

0.238C0 

0,26632 

0,28491 

0,29793 

0,30752 

0,31486 

0,32063 

0,32913 

1,2 

0,14139 

0,16116 

0,17437 

0,18377 

0,19077 

0,19616 

0,20045 

0,20681 

1.4 

0.086499 

0,096809 

0,10393 

0,10911 

0,11304 

0,11612 

0,11859 

0,12231 

1.6 

0.066661 

0,059859 

0,062258 

0,064103 

0,065557 

0,066729 

0,067692 

0,069178 

1.8 

0,038991 

0,038401 

0,038158 

0,038067 

0,038046 

0,038059 

0,038088 

0,038165 

2,0 

0,026848 

0,025061 

0.023871 

0,023029 

0,022403 

0.021920 

0,021538 

0,020972 

2,2 

0,017934 

0,016280 

0,015071 

0.014158 

0,013446 

0,012878 

0,012413 

0,011702 

2,4 

0,011635 

0.010476 

0,009557 

0,008826 

0,008237 

0,007754 

0,007353 

0,006726 

2,6 

0.0C7994 

0,007163 

0,006465 

0,005894 

0,005425 

0,005038 

0,004709  j 

0,004193 

2,8 

0,006491 

0.004893 

0,004382 

0.003960 

0,003811 

0.003323 

0.003077  i 

0,002690 

3.0 

0,003969 

0,003492 

0,003096 

0,002772 

0.002508 

0,002291 

0,002106  j 

0.001816 

3,2 

0,002997 

0,002592 

0,002271 

0.002016 

0,001810 

0,001644 

0,001502  j 

0,001283 

3.4 

0.002322 

0,001977 

0,001714 

0,001510 

0,001348 

0,001219 

0.001103 

0,000940 

3.6 

0,001824 

0,001535 

0,001320 

0,001156 

0,001027 

0,000926 

0,000839 

0,000708 

3.8 

,  0.001445 

0,001206 

0,001031 

0,000900 

0,000797 

!  0,000718 

0,000648 

4.0 

!  0.001153 

0,000957 

0,030816 

0,000710 

0.C00628 

0,000565 

0,000510 

_ 

4.2 

!  0,000928 

0,000788 

0,000654 

0,000568 

0, 009502 

1  0,000452 

0.000407 

4,4 

.  0.000754 

J 

0.000623 

0,000530 

0,000480 

0,000106 

0,000366 

0,000329 

4.6 

0,000618 

0,000511 

0,000434 

0,000377 

0,000333 

0,030300 

0,000269 

4.8 

j  0,000512 

0,000423 

0,0,0359 

0,000312 

0,000275 

i  0,000249 

0,000223 

_ 

5,0 

[  0,000427 

0,000353 

0,000300 

0,000260 

0,000230 

|  0,000208 

0,000186 

- 

6.2 

i  0,000360 

0,000297 

0,000252 

0,000219 

0,000194 

j  0,000176 

0,000156 

- 

5,4 

j  0,000305 

0,000252 

0,000214 

0,000186 

0,000164 

1  0,003150 

0.000133 

5,6 

!  0,000261 

0,000215 

0,000183 

0,000159 

0,000140 

|  0,000129 

5,8 
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Table  9.  Coefficients  of  transmission  ts  and  T.  for  Plexiglas  of 
*^.=0°  [d  (g/cml)  the  value  of  space  in  the  depth]. 

8,o  d  =0,146  |  £<,*=4,0  d= 0,098  |  £o^2,0  d -0,0475  jso-l.O  d-0,0213  £0  =  0,7  d-0.0136  £o>=0,4  d=0. 00625 
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